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PRESENTACIÓN 


La resolución de problemas con el microordenador reporta varias 
ventajas, unas de carácter general o formativo y otras de carácter 
específico. 


Entre las primeras cabe destacar la adquisición de un hábito de 
rigor lógico en el enfoque y planteamiento de cualquier problema 
relacionado con determinada área de conocimiento o con la vida 
real. Este hábito se afianza según se va utilizando el lenguaje de 
programación con el que nos comunicamos con el microordena- 
dor, en cuanto que dicho lenguaje lteva a una estructuración del 
pensamiento, ante la necesidad de descomponer el problema en 
una serie de pequeños pasos perfectamente relacionados. 


Entre las ventajas de carácter específico hay que señalar el hecho 
de que la actividad intelectual que se desarrolla en la resolución 
de un problema concreto con el microordenador exige un discer- 
nimiento más profundo y riguroso de los conceptos vinculados al 
problema. En este sentido, la utilización del microordenador 
ofrece un método de estudio y un estilo de analizar situaciones 
que se caracterizan por ser más exhaustivos. Por ejemplo, la reso- 
lución de un problema de electricidad con el microordenador 
lleva a estudiar previamente esta parte de la física con más in- 
terés y rigor que si la resolución se hiciera por métodos tradicio- 
nales. 


Con este libro, PROGRAMAS DE APLICACIONES EN BASIC, se 
pretende ofrecer una panorámica de las posibilidades que brinda 
el microordenador en el estudio de las distintas asignaturas y en 
la resolución de problemas. 


El análisis de los programas de aplicaciones incluidos en el libro 
suscitarán en el lector el interés por elaborar otros programas 
para resolver situaciones que se le vayan presentando. Esto con- 
tribuirá a afianzar el rigor lógico y a obtener resultados prácticos 
en el estudio de las distintas áreas de conocimiento. 


LOS AUTORES 


En este libro se han utilizado las instrucciones de gráficos 
del estándar MSX. 


A A A A A A 


1. Introducción 


1. Objetivos de las aplicaciones 


Con las distintas aplicaciones incluidas en este libro pretendemos ofrecer una 
muestra de la gran variedad de problemas que se pueden resolver con un mi- 
croordenador. 


Por otra parte, muchas aplicaciones, sobre todo las más sencillas, pueden servir de 
ejercicios de programación. En este caso, sugerimos comprobar los programas 
que se elaboren ejecutándolos en el microordenador, y a continuación compa- 
rarlos con los que transcribimos en el libro. 


2. Algunas pautas para resolver problemas con el microordenador 


La solución de cualquier problema en informática constituye un proceso de inves- 
tigación y, como tal, no se puede establecer un conjunto de reglas fijas que apli- 
cándolas conduzcan a una solución correcta. Sin embargo, ciertas pautas o con- 
sejos pueden ser útiles en determinadas ocasiones. 


Resumiremos a continuación algunas de estas ideas que facilitan el trabajo. 


e Definición del problema 


Se trata de simplificar y aislar el problema en que estamos interesados, elimi- 
nando toda la información no necesaria, para utilizar así el menor número posible 
de conceptos. 


e Análisis del problema ' 

Una vez definido el problema debemos elegir una estrategia para llegar a una so- 
lución. Como se va a emplear un microordenador es necesario imponer una serie 
de condiciones, que vienen dadas por las instrucciones que éste es capaz de «en- 
tender» y realizar. Este proceso nos llevará así a un procedimiento, es decir, a un 
número finito de pasos (algoritmo), que permiten resolver el problema. 


e Organigramas o diagramas de flujo 


Los distintos pasos que supone el análisis de un problema se pueden indicar me- 
diante esquemas gráficos. El conjunto de estos pequeños pasos vinculados entre 
sí y expresados gráficamente constituye lo que se llama organigrama o diagrama 
de flujo. 


En un diagrama de flujo, las diferentes etapas a ejecutar se representan en recua- 
dros en los que se incluyen las instrucciones específicas. El orden de ejecución de 
dichas etapas se indica mediante el sentido de las flechas que unen los recuadros. 


Como ejemplo, consideremos el diagrama de flujo correspondiente a un pro- 
grama que lee un número N y calcula el producto de los N primeros números na- 
turales si N no es par y en caso contrario no hace nada. 


Los símbolos que se han utilizado son: 


comienzo o final (END) 
del programa 


lectura 
(INPUT) 


escritura 
(PRINT) 


asignación 
(LET) 


toma de decisión 
(IF... THEN ...) 


repetición de las instrucciones 
intermedias 
(FOR ... TO ... STEP ...; NEXT ...) 


Como se pudo observar, a cada símbolo le corresponde una instrucción en BASIC 
(excepto a COMIENZO), por lo que es bastante fácil pasar del organigrama al co- 


rrespondiente programa. 


10 10 
20 20 
30 

30 

40 
E 
HE 
a 
70 70 


INPUT N 


IF N=2*INT(N/2) THEN 80 


LET P=1 


FOR |I=1 TO N 


LET P=Px] 
NEXT | 


PRINT P 


END 


La instrucción 20 averigua si N es par o no. Para ello se utiliza la función INT 
(parte entera). Por ejemplo, si mediante INPUT se introduce N=25, se tendría: 


IF N=2 * INT(25/2) THEN GOTO 80 
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El microordenador hace estos cálculos: 


2 * INT(12.5)=2 * 12=24 


y a continuación compara 25 con 24 (¿N=24?). Como no se cumple la igualdad 
(lo que supone que N no es par), ejecuta la instrucción siguiente a la 20. 


Elaboración de un programa 


Para concretar el procedimiento, pueden tenerse en cuenta los siguientes pasos: 


a) 


b) 


c) 


Traducir todas las informacionés que manejamos a variables y conjuntos de 
variables. 


Traducir todas las manipulaciones de datos a estructuras. Es decir, subdividir 
el programa en posibles bucles y subrutinas. 

La utilización de subrutinas facilita una concepción y presentación estructu- 
rada de los programas. Cada programa puede considerarse, entonces, como 
una asociación de bloques de instrucciones, constituyendo uno de ellos el 
programa principal y los restantes, bloques anexos (las subrutinas o subpro- 
gramas). 


Tratar de desarrollar el programa en un conjunto de bloques lo más com- 
pletos posible en sí mismos. (En este sentido es útil usar REM para identificar 
las distintas partes.) 


y 


Corrección de errores en un programa 


Los programas construidos raramente funcionan la primera vez que se ejecutan. 


El propio intérprete BASIC envía mensajes para señalar ciertos tipos de error (por 
ejemplo, de sintaxis), pero no detecta los errores lógicos cometidos. 


Para solventarlos es necesario estudiar la ejecución del programa estableciendo 
etapas que se van analizando. 


Esto puede conseguirse intercalando en el interior del programa instrucciones 
STOP (que detienen en un lugar determinado su ejecución), para posteriormente 
continuar hasta la siguiente etapa (siguiente STOP) utilizando la instrucción 
CONT. 
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En este proceso de comprobación podemos visualizar los valores que contienen 
las variables utilizando instrucciones PRINT en «modo directo». Así, 


PRINT A(20) 


permite conocer el contenido de la posición 20 de la variable A. 


Debe tenerse en cuenta que RUN borra automáticamente el contenido previo de 
todas las variables en memoria cuando se inicia la ejecución del programa. 


La instrucción CLEAR borra también el contenido de las variables, aunque no 
inicia la ejecución del programa. 


Para evitar este efecto puede utilizarse la instrucción GOTO es «modo directo» 
que permite transferir la ejecución del programa a la línea que se desee. 


e Verificación del programa 
a) Para detectar errores en un programa, podemos ejecutarlo con datos para los 


cuales se conocen los resultados y compararlos a continuación con las solu- 
ciones obtenidas. 


b) Si un programa contiene distintas bifurcaciones, debemos utilizar datos que 
permitan comprobar todas las opciones. 


c) La técnica de construcción de un programa en diversos módulos facilita su ve- 
rificación, ya que usando instrucciones de entrada y salida adicionales (que 
luego se eliminarán) permiten verificar los procesos intermedios. 


d) Una técnica de depuración útil (sobre todo en programas cortos) es seguir los 
pasos del programa con lápiz y papel, simulando el trabajo del microordenador. 


3. Ejercicios resueltos 


1. Expresar en un organigrama el proceso que hay que seguir para 
sumar los N primeros números naturales. 
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Solución 


2. Obtener un organigrama en el que se indique el cálculo de las raíces 
reales de una ecuación de segundo grado. 


Solución 


Las letras A, B, C representan los coeficientes de la ecuación A x? + 
Bx+C€=0 


LAS RAICES 
SON 
COMPLEJAS 


LAS RAÍCES 
SON: R1,R2 
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Expresar en un organigrama el proceso que hay que seguir para ave- 
riguar si un número es primo o no. 


Solución 


El símbolo Elx] representa la parte entera de X. 


Para averiguar si N es primo tenemos que comprobar que no es divi- 
sible por los números que van de 2 a VN (en realidad sólo por los 
que son primos): 


DESDE 1=2AvYyN 


SIGUIENTE I 


ES PRIMO 


4. Realizar un organigrama para obtener el número menor de una lista 
de 50 números. 


Solución 


Las variables que se utilizan son: 


— N (guarda el número leído) 
— M (guarda el número menor hasta el momento) 
— C (cuenta los números leídos) 


Un organigrama sin contador: 


Solución 


5. Encontrar un organigrama para calcular la raíz cuadrada entera de un 
número positivo. 


Solución 


La raíz cuadrada entera de un número R es un número x que elevado 
al cuadrado más el resto r es igual a R. 


VR=x=x?+r=R 


X ES LA RAÍZ 
EXACTA DER 


NI BS EA 
RAÍZ ENTERA 
DER 
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2. Dibujos a partir de la 
circunferencia 


1. Resumen de las instrucciones específicas que se van a utilizar 
en este capítulo 


Las instrucciones específicas que se van a utilizar en este capítulo se resumen a 
continuación: 


SCREEN 2 Pasa a modo gráfico. 
PSET (X, Y) Imprime un punto en el pixel (X, Y). 


LINE - STEP (X, Y) Dibuja un segmento desde el último pixel 
impreso al pixel que se encuentra X posi- 
ciones a la derecha e Y posiciones hacia 
abajo. 


CIRCLE (X, Y), R Dibuja una circunferencia de centro el 
pixel (X, Y) y de radio R. (En los microor- 
denadores IBM, MSX, etc., en lugar de cir- 
cunferencia, dibuja una elipse.) 


n GOTO n Mantiene la imagen en la pantalla. (En los 
microordenadores MSX, presionar [ CTRL] 
y [STOP |, para salir del bucle sin fin.) 


2. Dibujo de una circunferencia «punto a punto» 


Una circunferencia se puede dibujar directamente con la instrucción CIRCLE. Pero 
puede interesar dibujarla punto a punto, para lo cual es necesario determinar pre- 
viamente las coordenadas de uno de los puntos de dicha circunferencia. 


Suponiendo que el radio de la circunferencia es R y que el origen de coordenadas 
coincide con su centro, las coordenadas (X, Y) de un punto cualquiera de la 
misma, se obtienen partiendo de la definición de seno y coseno. 
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sen a => Y=Rsena 


cos QU Ho X=Rcosa 


Si el centro C de la circunferencia se sitúa en el punto (A, B), las coordenadas (X, 
Y) se obtendrán sumando a R cos a el número A y a R sen a, el número B. 


(0,0) 


Nota: En muchos microordenadores, como, por ejemplo, los MSX, las coorde- 
nadas cartesianas no coinciden con las coordenadas de la pantalla, pues 
los valores del eje vertical se ordenan de arriba abajo. 


' 
191-y4-- - —- (x,191-y) 


cartesianas pantalla 
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La relación entre ambas coordenadas es la siguiente: 
Al punto (X, Y) le corresponde el pixel (X, 191 — Y). 


No obstante, al ser la circunferencia simétrica respecto del diámetro horizontal, no 
es necesario tener en cuenta esta relación. 


Si decidimos dibujar la circunferencia con N puntos, tendremos que fijar previa- 
mente el valor del ángulo a. necesario para pasar de un punto al siguiente. Este 
valor se calcula teniendo en cuenta que el ángulo completo correspondiente a 
toda la circunferencia es 360”, que equivale a 2 rt radianes; luego si dibujamos N 


untos, tendremos: 
p , 3 2 


a.=2 * PI/N 


Es decir, el ángulo correspondiente al punto 1 valdrá 2 + PI/N x 1, el correspon- 
diente al punto 2, 2 * PI/N * 2, ..., y el correspondiente al punto N, 2 * PI/N * N = 
=2*Pl 

Luego, las coordenadas de los sucesivos puntos de la circunferencia se podrán de- 
terminar mediante las siguiente fórmulas: 


X=A+R*COS (2 + PI/N + 1) 
Y =B +R SIN (2 + PI/N * 1) 


siendo | = 1, 2, ...., N. 


Para elaborar el programa que dibuje una circunferencia con una serie de puntos, 
hay que tener en cuenta las dos últimas fórmulas, fijando previamente los si- 
guientes parámetros: 
A, B: coordenadas del centro de la circunferencia 
R: radio de la circunferencia 
N: número de puntos 


y tomar Pl = 3.1416, con lo que 2 * Pl = 6.2832. 
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Además, hay que definir la variable |, que tomará todos los valores enteros desde 
1 a N, ambos inclusive. Así se obtiene el siguiente programa: 


INPUT “CENTRO”; A, B 

INPUT “RADIO”; R 

INPUT "NUMERO DE PUNTOS”; N 
SCREEN 2 

FOR | = 1TON 


X = A +R* COS (6.2832/N + 1) 
Y = B +R * SIN (6.2832/N + 1) 
PSET (X, Y) 

NEXT | 

GOTO 200 


Si lo ejecutamos para A = 127, B = 95, R= 50 yN = 100, obtenemos: 


NOTA: 

En esta fotografía y en las restantes de este capítulo no aparecen verdaderas cir- 
cunferencias sino elipses. Ello es debido a que en los microordenadores MSX la 
separación de dos puntos consecutivos de la pantalla en posición horizontal es 
1.4 veces mayor que la separación de los que se encuentran en posición vertical. 
Por esta misma razón, tampoco aparecen verdaderos polígonos regulares en las 
fotografías del apartado 4, Polígonos regulares. 


(Ver Dibujo de «verdaderas» circunferencias y de «verdaderos» arcos de circunfe- 
rencia, páginas 282 y 283 del libro MSX. Programación, Gráficos, Colores y 
Música, de los mismos autores). 
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3. División de un círculo en N partes iguales 


Para dividir una circunferencia en N partes iguales, la dibujaremos previamente 
con una instrucción CIRCLE y, a continuación, uniremos su centro con cada uno 
de los N puntos; así se obtendrán N sectores iguales. 


El programa que realiza esta división es el siguiente: 


INPUT “CENTRO”; A, B 
INPUT “RADIO”; R 

INPUT “PARTES”; N 
SCREEN 2 

CIRCLE (A, B), R 

FOR | =1TON 

PSET (A, B) 

X =R * COS (6.2832/N * |) 
Y =R * SIN (6.2832/N + 1) 
LINE — STEP (X, Y) 

NEXT | 

GOTO 200 


Si lo ejecutamos para A = 127, B = 95,R= 50yN = 7, obtenemos: 
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Si lo ejecutamos para A = 127, B = 95,R=95yN = 21, tenemos: 


4. Polígonos regulares 


La división de la circunferencia en N partes iguales sugiere la obtención de polí- 
gonos regulares de N lados inscritos en la misma. 


El siguiente programa dibuja polígonos regulares pero no la circunferencia en la 
que están inscritos. Esta circunferencia tiene su centro en el pixel (127, 95) y su 
radio es 80. 


INPUT “NUMERO DE LADOS”; N 
SCREEN 2 

A = 127 + 80 

B=95 

PSET (A, B) 

FOR | =1TON 

X = 127 + 80 * COS (6.2832/N + 1) 
Y = 95 + 80 + SIN (6.2832/N + 1) 


LINE — STEP (X — A, Y — B) 
A =X 

B=Y 

NEXT | 

GOTO 200 
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Con la instrucción 10 se introduce el número de lados del polígono. 


Las instrucciones 30, 40 y 50 imprimen un punto en el pixel (127 + 80, 95), 
donde comienza a dibujarse el polígono. 

El bucle 60 — 120 dibuja los N lados. 

Las instrucciones 70 y 80 calculan las coordenadas del siguiente punto de la cir- 
cunferencia que hay que unir al A, B. 

La instrucción 90 LINE - STEP (X — A, Y — B) une los pixels (A, B) y (X, Y). 


Las instrucciones 100 y 110 almacenan en A y B las coordenadas X, Y, calculadas 
por las instrucciones 70 y 80, que serán el pixel de partida para dibujar el si- 
guiente lado. 


En resumen, (A, B) representan las coordenadas del origen de un lado y (X, Y), el 
punto final del lado. 


Si lo ejecutamos, tenemos 
TRIÁNGULO (N = 3) 


CUADRADO (N = 4) 


PENTÁGONO (N = 5) 


HEXÁGONO (N = 6) 


HEPTÁGONO (N = 7) 
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OCTÓGONO (N = 8) 


DODECÁGONO (N = 12) 


POLÍGONO REGULAR DE 30 LADOS 


Nota: Como se puede observar, al aumentar el número de lados, el perímetro 
del polígono regular se aproxima cada vez más a la circunferencia. 
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5. Flor de circunferencias 


La flor de circunferencias es un conjunto de N circunferencias, todas ellas de radio 
R, cuyos centros determinan a su vez otra circunferencia de radio RP, llamado ra- 
dio principal. 

En el programa siguiente, que genera la flor de circunferencias, la circunferencia 
de radio RP tiene su centro en el pixel (127, 95). 


INPUT “NUMERO DE CIRCUNFERENCIAS”; N 
INPUT “RADIO DE LAS CIRCUNFERENCIAS”; R 
INPUT “RADIO PRINCIPAL”; RP 

SCREEN 2 j 

FOR | = 1 TON 

X= 127 + RP * COS (6.2832/N + 1) 

Y = 95 + RP * SIN (6.2832/N * 1) 


CIRCLE (X, Y), R 
NEXT | 
GOTO 200 


Cuando se ejecuta este programa, al introducir los datos hay que evitar que el di- 
bujo se salga de la pantalla, para lo cual debe cumplirse R + RP < 95. 


Flor para N = 6,R = 20 y RP = 67 


30 


Flor para N = 12, R = 20 y RP = 67. 


20 y RP = 67. 


40, R = 


Flor para N 


ón NY 
PEN y 
AnS 


Flor para N = 8, R = 30 y RP = 35. 
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Flor para N = 16, R = 40 y RP = 47. 


Flor para N = 30, R = 25 y RP = 25. 


Flor para N = 20, R = 40 y RP = 30. 


de 


6. Rayos desde el centro de una circunferencia 


El programa siguiente traza N rayos de longitud 70, que parten del pixel (127, 95). 


INPUT “NUMERO DE RAYOS”; N 
SCREEN 2 

FOR|=1TON 

PSET (127, 95) 

X = 70 * COS (6.2832/N + |) 


Y = 70 + SIN (6.2832/N + |) 
LINE — STEP (X, Y) 

NEXT | 

GOTO 200 


Para N = 20 se obtiene: 


Para N = 50, tenemos: 


33 


7. Rayos que parten de un punto situado a la izquierda de una 
circunferencia 


Localizando el centro emisor de los rayos del programa anterior en el pixel (27, 
95) (para que quede suficiente espacio a su derecha) y aumentando X en 100 
unidades en la instrucción LINE - STEP, se obtiene un conjunto de rayos que par- 
ten de un punto situado a la izquierda del centro de la circunferencia de radio 70, 
y centro el pixel (127, 95). 


INPUT “NUMERO DE RAYOS”; N 
SCREEN 2 

FOR |=1TON 

PSET (27, 95) 

X = 70 * COS (6.2832/N + l) 


Y = 70 * SIN (6.2832/N + |) 
70 LINE - STEP (100 + X, Y) 
80 NEXT! 
200 GOTO 


Para N = 40, tenemos: 


Para N = 100, obtenemos: 
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A E E A | 


3. Representación gráfica 
de funciones 


1. Planteamiento del problema 


e Para dibujar la gráfica de una función debemos tener en cuenta dos cuestiones 
previas: 


— ¿Dónde situar el centro de los ejes? 
— ¿Qué escala hay que elegir en cada eje? 


Con el fin de contestar a estas preguntas estudiaremos una función sencilla: F 
(x) = x?, y a partir de ella intentaremos sacar conclusiones generales. 


Se supone que la pantalla del microordenador, en el modo gráfico, tiene 256 
líneas verticales (de O a 255) y 192 líneas horizontales (de O a 191); luego en 
la pantalla se pueden localizar 256 x 192 = 49152 puntos. 


A cada punto de la pantalla se le llamará pixel. 


También se supone que en la pantalla en modo gráfico cada pixel está referido 
a un sistema de coordenadas, cuyo origen está localizado en el vértice superior 
izquierda (columna 0, fila 0).. A los demás vértices les corresponden las coorde- 
nadas indicadas en la figura. 


0123 A Socios 255 
0.0) (255.0) 
2 
3 
191 
(0,191) (255,191) 


El centro de la pantalla corresponde entonces al pixel (127, 95). 
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e Las instrucciones específicas que se van a utilizar se resumen en esta tabla: 


Instrucción Función 


SCREEN 2 Pasa a modo gráfico. 
PSET (X, Y) Imprime un punto en el pixel (X, Y). 
LINE (A, B) — (X, Y) Dibuja un segmento desde el pixel (A, B) al 


pixel (X, Y). 

n GOTO n Mantiene la imagen en la pantalla. 
(En los microordenadores MSX, presionar las 
teclas y para salir del bucle 


sin fin.) 


2. Centrado de ejes 


La gráfica de la función F (x) = x? es una parábola cuyo vértice es el punto (0, 0), 
y ocupa el primero y el segundo cuadrante. Luego los ejes deben estar centrados 
en el pixel (127, 191) que se encuentra en el centro de la línea inferior de la pan- 
talla. 


La subrutina que dibuja los ejes es: 


REM CENTRADO DE EJES 
SCREEN 2 

LINE (0, 191) — (255, 191) 
LINE (127, 191) — (127, 0) 
GOTO 2000 


La instrucción 1020 dibuja el eje de abscisas, y la instrucción 1030 dibuja el eje de 
ordenadas. 


Al ejecutar esta subrutina se obtiene lo que muestra la pantalla... 
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Si queremos que los ejes queden centrados en el pixel (C1, C2), basta generalizar 
la subrutina anterior cambiando 127 por C1 y 191 por C2 en las instrucciones 
1030 y 1010, respectivamente. 


REM CENTRADO DE EJES 
SCREEN 2 

LINE (0, C2) — 255, C2) 
LINE (C1, 191) — (C1, 0) 


GOTO 2000 


Esta subrutina dibuja los ejes coordenados para cualquier centro (C1, C2), pero 
aquí solamente consideraremos los siguietes: 


e Centro (0, 191): determina el primer cuadrante 


Para conseguir los ejes con este centro hay 
que añadir estas instrucciones: 


20 READ C1, C2 


9000 REM DATOS 
9010 DATAO, 191 


e Centro (127, 191): determina el primero y el segundo cuadrante 


Para conseguir los ejes con este centro hay 
que añadir estas instrucciones: 


20 READ C1, C2 


9000 REM DATOS 
9010 DATA 127, 191 
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e Centro (0, 95): determina el primero y el cuarto cuadrante: 


Para conseguir los ejes con este centro hay 
que añadir estas instrucciones: 


20 READ C1, C2 


9000 REM DATOS 
9010 DATAO, 95 


e Centro (127, 95): determina los cuatro cuadrantes 


Para conseguir los ejes con este centro hay 
que añadir estas instrucciones: 


20 READ C1, C2 


9000 REM DATOS 
9010 DATA 127, 95 


3. Localización de la gráfica en los límites de la pantalla 


e Si una vez elegido el centro de los ejes, se empieza a dibujar puntos de la fun- 
ción F (X) = X?, puede ocurrir que algunos de ellos queden fuera de los l- 
mites de la pantalla. Ocurre, entonces, que la ejecución se detiene, apare- 
ciendo un mensaje de error. Es necesario, por tanto, evitar esta situación 
eligiendo adecuadamente los valores externos de X a la izquierda y derecha de 
la pantalla, de tal manera que sus correspondientes valores F (X) no sobrepa- 
sen los valores 191 y O (límite inferior y superior de la pantalla). 


po 20) 


191 
(127,191) 
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Sr el valor máximo que puede tomar F (X) es 191, ¿cuáles serán los valores ex- 
tremos de X que hacen que F (X) no supere dicho valor? 


Si F(X) = X? = 191 entonces X = + V191 = + 13.8203. 
Luego los valores de X estarán comprendidos entre —13 y 13: 
=13<X< 13 


Los valores extremos de X a la izquierda y a la derecha del origen de coorde- 
nadas los indicaremos por X1 y X2, respectivamente. En el ejemplo que ve- 
nimos estudiando, X1 = -13 y X2 = 13. Fijados estos extremos, ya se puede 
dibujar la gráfica de la función F (X) = X*, mediante el siguiente programa: 


DEF FN F (X) =X 7 2 
READ C1, C2 

GOSUB 1000 

READ X1, X2 

REM *xxxx*xx 

REM DIBUJO 

FOR X = X1 TO X2 

Y = FNF (X) 

PSET (C1 + X, C2 — Y) 
NEXT X , 
GOTO 200 

REM xxx x*xx*x**x 
REM CENTRADO DE EJES 
SCREEN 2 


LINE (0, C2) — (255, C2) 
LINE (C1, 191) — (C1, 0) 
RETURN 

REM xxxx 

REM DATOS 

DATA 127, 191 

DATA -13, 13 
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Este programa funciona así: 
— La instrucción 10 define la función FN F (X) =X ? 2 
— La instrucción 20 lee los valores del centro C1 y C2, y la instrucción 30 llama 


a la subrutina de «centrado de ejes» para que dibuje los ejes de coordenadas 
centrados en C1, C2. 


— La instrucción 40 lee los valores extremos del eje de abscisas X1 y X2. 


— El bloque 100 - 140 dibuja la función. 
Al ejecutar el programa se obtiene la siguiente gráfica: 


Se observa que la gráfica obtenida está formada por puntos bastante distantes, y 
que se ha desperdiciado gran parte de la pantalla, a izquierda y a derecha del eje 
de ordenadas. Esto se puede evitar eligiendo adecuadamente las escalas en 
ambos ejes. 


4. Elección de escalas en los ejes 


e Para determinar la escala en cada eje elegiremos las siguientes variables: 


X1: extremo inferior del eje de abscisas 

X2: extremo superior del eje de abscisas 

Y1: extremo inferior del eje de ordenadas 

Y2: extremo superior del eje de ordenadas 

DX: valor de la división en el eje X 

DY: valor de la división en el eje Y 

X3: abscisa del punto de la pantalla que se va a imprimir 

191 - Y3: ordenada del punto de la pantalla que se va a imprimir. 
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El programa que incluye la utilización de escalas es el siguiente: 


DEF FN F (X) =X 1 2 
READ C1, C2 

GOSUB 1000 

READ X1, X2 

READ Y1, Y2 

DX = (X2 — X1)/255 

DY = (Y2 —- Y1)/191 

REM xxxxxx 

REM DIBUJO 

FOR X = X1 TO X2 STEP DX 
Y = FNF (X) 

X3 = INT ((X — X1) / DX + 0.5) 
Y3 = INT ((Y — Y1) / DY + 0.5) 
PSET (X3, 191 — Y3) 

NEXT X 

GOTO 200 

REM *xxxx*xx*xx* 

REM CENTRADO DE EJES 
SCREEN 2 

LINE (0, C2) — (255, C2) 
LINE (C1, 191) — (C1, 0) 
RETURN 

REM *xx*xxx*x 

REM DATOS 

DATA 127, 191 

DATA -13, 13 

DATA O, 169 


Veamos cómo funciona este programa: 

— La instrucción 10 define la función F (X) = X ? 2. 

— La instrucción 20 lee las coordenadas del centro (C1, C2). 
— La 30 llama a la subrutina que dibuja los ejes. 


— La 40 lee los valores extremos del eje de abscisas y la 50 los valores extremos 
del eje de ordenadas. 
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Como se ha visto, estos valores extremos no son arbitrarios sino que se calculan 
teniendo en cuenta los límites de la pantalla. 


En este ejemplo, X1 y X2 se determinaron anteriormente: 
X1 = -13 y X2 = 13 

Los extremos de Y se calculan considerando que 
F (-13) = (-13) =:169 y F (13) = 13? = 169 

y que se cumple F (X) = X? > 0 

Luego, 
Y1 =0e Y2 = 169 


— Las instrucciones 60 y 70 calculan el valor de la división en los ejes de abscisas 
y ordenadas, respectivamente. 


_ AZ2=XI_ 133218) . 26... 
DX = E = 0.10196078 
MW =YT . 1690 
19 191 


DY = = 0.88481675 


Luego en el eje X el salto de 1 pixel va a corresponder a un salto real de 
0.10196078 (se ensancha casi 10 veces la gráfica), y el eje Y va a correspon- 
der a un salto real de 0.88481675. 


— La línea 110 FOR X = X1 TO X2 STEP DX abre el bucle que dibuja la fun- 
ción. En total se van a imprimir 255 pixels, pues X varía desde X1 = —13 
- hasta X2 = 13 pero con paso DX = 0.10196078. 


— La instrucción 120 Y = FN F (X) calcula los valores de la función F (X) = X? 
correspondientes a los valores comprendidos entre X1 y X2. 


— Las instrucciones 130 y 140 calculan los valores X3 e Y3 en función de X e Y, 
respectivamente. Aquí se tienen en cuenta los extremos X1, X2 y las escalas 
elegidas DX y DY en los ejes coordenados. 


Los valores de la variable X3 estarán comprendidos entre 0 y 255, y se calcu- 
lan mediante la expresión: 

XX! 
DX 
Los valores de Y3 estarán comprendidos entre O y 169, y se calculan mediante 

la expresión: 

Y —- Y] 
DY 

— La instrucción 150 PSET (X3, 191 — Y3) imprime el pixel correspondiente. 


— La línea 160 NEXT X obliga a repetir el bucle. 


X3 = INT ( 


+ 0.5) 


Y3 = INT ( + 0.5) 
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Si se ejecuta el programa se obtiene lo que muestra la pantalla: 


Ahora se aprecia mejor la gráfica de la función F (X) = X*, aunque tiene el in- 
conveniente de estar bastante deformada por haber elegido escalas diferentes en 
el eje de abscisas y en el de ordenadas. 


169 


169 pixels 


-13 0 13 
| 255 pixels | 


La gráfica está ensanchada considerablemente, y este ensanchamiento viene dado 
por la relación entre ambas escalas DX y DY. Dicha relación es de 0.10196078 a 
0.88481675. 


0.10196078 / 0.88481675 = 0.11523378 
que impone un ensanchamiento de casi 10 veces. 
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5. Gráficas de algunas funciones 


e F(X) = X? 
Como F (X) puede tomar valores negativos y positivos interesa centrar los ejes 
en el pixel (127, 95). 


0 


(127,95) 


191 


Si hacemos F (X) = X? = 95, se obtiene X = V 95 = 4,562903. 
Elegiremos, entonces: 


Xl =-4 X2=4 
Y1 = -64 Y2 = 64 


En el programa general hay que cambiar las instrucciones siguientes: 


10 DEFFNF (X) =X 7 3 
9010 DATA 127, 95 
- 9020 DATA -4, 4 
9030 DATA -64, 64 


Al ejecutar el programa se obtiene la siguiente pantalla: 


e F(X) = sen X 
La función F (X) =sen X es periódica, de período 2 mn. Luego X variará entre 
X1 =0yX2 = 2 1 = 6.2832. Como, además, el seno varía entre —1 y 1, 
-=1<senX< 1 
los valores extremos de Y serán Y1 = —1 e Y2 = 1. Es evidente que el centro 
habrá que localizarlo en el pixel (0, 95). 
Las instrucciones que hay que cambiar en el programa general son éstas: 


DEF FN F (X) = SIN (X) 
DATA O, 95 

DATA 0, 6.2832 

DATA -1, 1 
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Al ejecutar el programa se obtiene la siguiente gráfica: 


e F(X) = cos X 


La función F (X) = cos X también es periódica, de período 2 n, luego X va- 
riará entre X1 = 0 y X2 = 2n = 6.2832. Como además se cumple —1 < sen 
X < 1, los valores extremos de Y serán los siguientes: Y1 = —1 e Y2 = 1. 


El centro se situará en el pixel (0, 95). 


Las instrucciones que hay que cambiar en el programa general son éstas: 


10 DEF FN F (X) = COS (X) 
9010 "DATA O, 95 
9020 DATA O, 6.2832 
9030 DATA -1, 1 
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Al ejecutar el programa se obtiene la siguiente gráfica: 


e F(X) = TAN X 


La función F (X) = TAN X también es periódica, pero de período y como es 
discontinua en 


e E 
2 2 
elegiremos como valores extremos de X los siguientes: 
XI === +0.1=-1.4708 
X2=  =0.1= 1.4708 
Y como 


F (X1) = TAN (= 3 + 0.1) = -9.9666444 


F(X2) = TAN (LT _ 0,1) = 9.9666444 


2 
hacemos 
Y1 = -10 
1.10 


El centro estará en el pixel (127, 95). 
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Las instrucciones que hay que cambiar en el programa general son éstas: 


DEF FN F (X) = TAN (X) 
DATA 127, 95 
DATA -1.4708, 1.4708 
DATA —10, 10 


La gráfica que se obtiene es la siguiente: 


e F(X)=VX 


La función F (X) = V X da el valor positivo de la raíz cuadrada de X, y sola- 
mente definida para X > 0. Se tiene entonces que si 


X1=0 
X2 = 255 


la raíz cuadrada de estos valores es 
1255 = 15.968719 
Luego los valores extremos de Y son: 


Y1=0 
Y2 = 15.97 


El centro estará en el pixel (0, 191). 
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Las instrucciones que hay que cambiar en el programa general son: 


DEF FN (X) = SQR (X) 
DATA O, 191 

DATA O, 255 

DATA O, 15.97 


La gráfica que se obtiene es la siguiente: 


F (X) = e* 


La función exponencial F (X) = e* siempre es positiva y al ser e” = 1, en el eje 
de ordenadas debemos tomar un valor pequeño para que se pueda apreciar la 
gráfica de la función para valores de X negativos. 


Si elegimos como valor extremo superior de Y 
Y2 =20 
el correspondiente valor extremo de X a la derecha se obtiene así 
* = 20 => In e* = In 20 >x = 2.9957323 
Esto nos lleva a elegir como valores extremos de X los siguientes: 


Xl =-3 
X2=3 


El centro será el pixel (127, 191). 
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Las instrucciones que hay que cambiar en el programa general son: 


DEF FN F (X) = EXP (X) 
DATA 127, 191 

DATA -3, 3 

DATA 0, 20 


El resultado que se obtiene al ejecutar el programa es el que muestra la foto- 
grafía: 


e F(X)=InX 


La función logaritmo natural o neperiano F (X) = In X sólo está definida para 
valores de X positivos y al ser In 1 = O el valor máximo de X no debe ser muy 
grande para que se pueda apreciar la parte negativa de la función. 


Por ejemplo, si elegimos como valores extremos de X 


X1 =0.1 
X2 = 20 
los correspondientes valores extremos de Y se obtienen así: 


In 0.1 = -2.3025851 
In 20 = 2.9957323 


El centro será el pixel (0, 95). 
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Las instrucciones que hay que cambiar en el programa general son: 


DEF FN (X) = LOG (X) 
DATA O, 95 

DATA 0.1, 20 

DATA -2.3, 3 


La gráfica que se obtiene es la siguiente: 


e F(X) “+ 


La función F (X) = L- es discontinua para X = 0, luego hay que evitar valores 
próximos al cero; además, X no debe ser ni muy pequeña ni muy grande, pues 
de lo contrario F (X) = - se haría demasiado grande o demasiado pequeña, 
respectivamente. Así se obtendrá una buena representación. 
Teniendo en cuenta las anteriores consideraciones, tomamos: 

X1 =-5 X2=5 

YI =-5 1123 
y como centro, el pixel (127, 95) pues X e Y van a ser ambos positivos. 
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Al dispararse la función para valores de X muy próximos al cero (discontinuidad 
en el cero) hay que añadir una instrucción como la siguiente: 


115 1F — 0.2 < = X AND X < = 0.2 THEN GOTO 160 
que evita la división por cero. 
Las instrucciones que hay que cambiar o añadir al programa general son: 


DEF FN (X) = 1/X 

IF 0.2 < = X AND X < = 0.2 THEN GOTO 160 
DATA 127, 95 

DATA -5, 5 

DATA -5, 5 


Al ejecutar el programa se obtiene esta gráfica: 


e F(X)= sa X 
La función F (X) = na 


es continua, pues cuando X se aproxima a cero (0), se tiene que F (X) = 


A EA se aproxima a 1. = pa = (.99833417, y cuando X se hace 
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grande (en valor absoluto), m2 


sen = 0.043352523. 


Teniendo en cuenta estas propiedades, elegimos como valores extremos, por 
ejemplo, los siguientes: 


Xl =-15 X2 =15 
11==] Y2=1 


El centro será el pixel (127, 95), pues X e Y pueden ser positivos y negativos. 
Las instrucciones que hay que cambiar en el programa general son: 


tiende a cero, pues el sen X se mantiene 


entre —1 y 1. Así: 
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DEF FN F (X) = SIN (X) / X 
DATA 127, 95 

DATA -15, 15 

DATA —1, 1 


La gráfica que se obtiene es ésta: 


ea) = EL 


Esta función es discontinua para X igual a cero, luego debemos añadir al pro- 
grama general una instrucción como ésta: 


115 IF-1<=X AND X < = 1 THEN GOTO 160 
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dado que si X tomara un valor entre —1 y 1 la función se saldría de la pantalla. 


Además, si X se hace grande (en valor absoluto) Sta 


a 2. = -0.050645861. 


Luego tomaremos por ejemplo los siguientes valores extremos: 


Xl =-15 X2=15 
Yl = 1 Y2 =1 


Como los valores de X e Y pueden ser positivos y negativos, el centro estará en el 
pixel (127, 95). 


Las instrucciones que hay que añadir o cambiar en el programa general son: 


tiende a cero, al mante- 


nerse cos X entre —1 y 1. Así: 


DEF FN F (X) = COS (X) / X 

IF—1 < = X AND X < = 1 THEN GOTO 160 
DATA 127, 95 

DATA -15, 15 

DATA -1, 1 


La gráfica que se obtiene es la que muestra la fotografía: 
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e F(X) = sen X + cos X 


La función F (X) = sen X + cos X es periódica, de período 2n, y al ser la suma 
del seno y del coseno, F (X) variará entre =1 -1 = -2y1+1=2. 


Los valores extremos que debemos tomar son: 


X1=0 X2 = 2n= 6.2832 
11 =-2 Y2=2 


Como X sólo toma valores positivos, el centro lo tomamos en el pixel (0, 95). 
Las instrucciones que hay que cambiar en el programa general son: 


10 DEF FN (X) = SIN (X) + COS (X) 
9010 DATA O, 95 
9020 DATA O, 6.2832 
9030 DATA -2, 2 


La gráfica que se obtiene es la siguiente: 
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4. Rectas y parábolas. Resolución 
gráfica y algebraica de ecuaciones 
de primer y segundo grado 


1. Representación gráfica de funciones afines, lineales y cons- 
tantes, y de las relaciones constantes (familias de rectas) 


Presentación del problema 


Transcribimos las fórmulas de la función afín, lineal y constante, así como de la 


relación constante. 


Función afín 
y=ax+b 


Función 
constante pS 
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Función 
lineal 
y = ax 


Relación x = K (K, constante) 


constante 


Recordar que a es el coeficiente angular o pendiente de la recta afín y lineal, y es 
el parámetro que da la inclinación de la recta. El parámetro b es el término inde- 
pendiente, y determina el punto del eje de las ordenadas por el que pasa la recta. 
Si b = 0, la función afín se convierte en /ínea/, y en este caso la recta pasa por el 
origen de coordenadas. 


Si a = 0, la función afín se convierte en constante, y en este caso la recta es para- 
lela al eje de abscisas. 


Sia>0 yb =0, la recta coincide con el eje de abscisas. 


La exprexión x = k (k, constante) no es una función; es simplemente una relación. 
La gráfica correspondiente es una recta paralela al eje de ordenadas que corta el eje 
de abscisas en el punto k del mismo. Si k = 0, la recta coincide con el eje de or- 
denadas. 


El programa 


El siguiente programa dibuja gráficas de funciones afines, lineales y constantes, y de 
relaciones constantes. 


También permite obtener en la pantalla familias de rectas, variando el coeficiente a y 
el término independiente b. El número máximo de rectas que pueden aparecen si- 
multáneamente en la pantalla es 4. 


Se ha puesto este tope porque con un número mayor, la pantalla puede aparecer 
confusa. 


10 RECTAS 
20 INPUT “PARA DIBUJAR RECTAS VERTICALES TECLEAR V Y RE- 
TURN. PARA DIBUJAR OTRAS RECTAS, PULSAR OTRA TECLA Y 
RETURN”; A$ , 
30 IFA$ = "V” THEN GOTO 1000 
40 REM DIBUJA RECTAS AFINES, LINEALES Y CONSTANTES 
50 DEF FNI (X) = A* X +B 
60 READ R 
710 K=10 
80 SCREEN 2 
90 OPEN “GRP:” AS 4 1 
100 GOSUB 2000 


La instrucción 60 READ R (y su correspondiente DATA) lee el dato R que da el 
número de rectas que se van a representar. 
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La instrucción 100 transfiere el control a la subrutina 2000 que dibuja los ejes de 
coordenadas. 


110 FORI=1TOR 

120 READA, B 

130 IF K = 10 THEN PRESET(K, 6 * (I-1)) : PRINT 4 1, “Y="; A; 
“X+"; B: K = 140: GOTO 150 

140 IF K <> 10 THEN PRESET(K, 6 * (1-2)) : PRINT 4 1, “Y="; A; 
"X+";B:K=10 

150 GOSUB 3000 

160 FORT = 1TO 1000: NEXT T 

170 NEXT | 


Las líneas 130 y 140 escriben en la parte superior de la pantalla de alta resolución la 
ecuación de la recta que se va a representar. 


La instrucción 150 transfiere el control a la subrutina 3000 para dibujar la gráfica co- 
rrespondiente. 


La línea 160 introduce un bucle de retardo de unos dos segundos antes de dibu- 
jarse la siguiente recta. 


180 PRESET(10, 183) : PRINT 4 1, “Para seguir, CTRL y STOP.Y RUN” 
190 GOTO 190 


La línea 180 escribe en la parte inferior de la pantalla de alta resolución el texto que 
indica cómo salir del bucle sin fin de la línea 190. 


200 DATA 
210 DATA 


En 200 DATA se pondrá el dato correspondiente al número R de rectas que se 
deseen representar y en 210 DATA, el mismo número de pares de datos corres- 
pondientes a A y B. 


1000 REM RECTAS VERTICALES 

1010 INPUT “INTRODUCE ABSCISA X”; X 

1020 SCREEN 2 

1030 OPEN "GRP:” AS+1 

1040 GOSUB 2000 

1050 PRESET(120,0) : PRINT 4 1, "X="; X 

1060 LINE (129 + X/0.1, 20) — (129 + X/0.1, 180) 

1070 PRESET(10, 180) : PRINT + 1, “Para seguir, CTRL y STOP.Y 
RUN” 

1080 GOTO 1080 
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Este bloque de instrucciones dibuja rectas verticales, una vez introducido el valor 
de la abscisa por la cual pasará la recta. 


En la línea 1060, la variable X está dividida por 0.1 debido a que dos divisiones 
consecutivas del eje de abscisas están separadas 10 puntos. 


2000 REM SUBRUTINA EJES 

2010 REM EJE DE ABSCISAS 

2020 LINE(O, 100) — (255, 100) 

2030 FOR! = 7 TO 247 STEP 10 

2040 PRESET(I, 97) : PRINT 4 1, “|” 
2050 NEXT | 

2060 PRESET(138, 106) : PRINT + 1, “1” 
2070 PRESET(247, 106) : PRINT + 1, "X” 
2080 REM EJE DE ORDENADAS 

2090 LINE (129, 20) — (129, 180) 

2100 FOR !|= 24 TO 170 STEP 10 

2110 PRESET(127, 1) : PRINT4 1, ” 
2120 NEXT | 

2130 PRESET(120, 20) : PRINT + 1, “Y” 
2140 PRESET(120, 86) : PRINT + 1, “1” 
2150 RETURN 


Esta subrutina dibuja los ejes coordenados con divisiones, señalando con X el eje 
de abscisas y con Y el de ordenadas, y poniendo un 1 en la primera división del 
eje X y del eje Y. 


3000 REM DIBUJA LA GRAFICA 

3010 IFA = 0 THEN X1 = -12 : X2 = 12 : GOTO 3060 
3020 X1= (-7-—BJ/A : X2 = (7 — BJ/A 

3030 IFA<OTHENC=X1:X1 = X2: X2 = C 

3040 IF X1 < -12 THEN X1 = -12 

3050 IF X2 > 12 THEN X2 = 12 


Con este bloque de instrucciones se determinan los valores extremos X1 y X2 
para evitar que los puntos de la gráfica queden fuera de la pantalla. 


3060 FOR X = X1 TO X2 STEP 0.05 
3070 Y = FNF(X) 

3080 X3 = INT(X/0.1 + 0.5) 

3090 Y3 = INT(Y/O.1 + 0.5) 

3100 PSET(129 + X3, 100 — Y3) 
3110 NEXT X 

3120 RETURN 
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Las instrucciones 3060-3100 dibujan la gráfica. En las líneas 3080 y 3090, las va- 
riables X e Y están divididas por 0,1 que establece la misma escala para los dos 
ejes. Se ha elegido este factor de escala debido a que dos divisiones consecutivas 
del eje de abscisas (y también del de ordenadas) están separadas 10 puntos. 


La instrucción 3100 dibuja los puntos, refiriéndolos al pixel 129, 100 (centro de 
coordenadas). 


Ejecución del programa 

Antes de iniciar la ejecución del programa hay que poner en las instrucciones 200 
y 210 los datos que se deseen. 

A continuación presentamos algunos casos: 


e Para representar gráficamente las funciones y = x, y = 2x + 3e y = —4x — 6.1, 
las instrucciones 200 y 210 hay que escribirlas así: 


200 DATA 3 
210 DATA 1,0, 2,3, -4, -6.1 


Tecleando RUN y | RETURN | y pulsando a continuación cualquier tecla que no 
sea V se obtiene este resultado. 


e Para obtener las gráficas de las funciones y = x, y =Xx+ 1, y =x+3ey=x- 
— 6.5, las instrucciones 200 y 210 hay que escribirlas así: 


200 DATA 4 
210 DATA 1,0, 1,1, 1,3, 1, 6.5 


Tecleando RUN y | RETURN | y pulsando a continuación cualquier tecla que no 
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sea V se obtiene una familia de rectas con la misma pendiente (rectas paralelas 
de pendiente 1). 


PARA SEGUIR,CTRL y STOP.Y RUN 


e Para representar gráficamente las funciones y = x + 3, y = 2x + 3, y = -2.5x + 
+ 3 e y = -x + 3, las instrucciones 200 y 210 hay que escribirlas así: 


200 DATA 4 
210. DATA. 1,5, 2,3, 203,9, 1,3 


Tecleando RUN y | RETURN | y pulsando a continuación cualquier tecla que no 
sea V se obtiene una familia de rectas que cortan al eje Y por el punto 3. 


RUN 


PARA SEGUIR,CTRAL yv STOP.Y 


e Si después de teclear RUN y | RETURN | se pulsa la tecla V, el programa conti- 
núa la ejecución en la instrucción, 1000 REM RECTAS VERTICALES detenién- 


dose en la línea 1010 INPUT “INTRODUCE ABSCISA X”; X. 
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Si se introduce el valor 5 para X en la pantalla se obtiene este resultado. 


2. Resolución gráfica y algebraica de un sistema de primer grado 
con dos incógnitas 

Planteamiento del problema 

e Resolución del sistema 


El sistema de ecuaciones 


P1X+0Q1Y =RÍ 
P2X + Q2 Y =R2 


Se puede resolver por distintos métodos. Aquí elegiremos el de reducción: 


(Xx Q2) P1Q2X +0182Y = R1Q2 
(x-Q1) =P2Q1X -—.Q492Y = -—R2Q1 
P1Q2X —P2Q1X =  R1Q2 — R2Q1 
(P1Q2  —P2Q1)X = R1Q2 -— R2Q1 
X = R1Q2 — R2Q1 
P1Q2 - P2Q1 
Análogamente se obtiene el valor de Y: 
y = PIR2 — P2R1 
P1Q2 - P2Q1 
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Si hacemos 


E = R1Q2 — R2Q1 
F =P1R2 — P2R1 
= P1Q2 — P2Q1 
y sustituimos en las igualdades recuadradas, se tiene 
E E 
X= 5 Y =>5 
e Discusión del sistema 
1. SID X= 0, los cocientes 5 + existen. En consecuencia, la solución 


del sistema es única; o sea, el sistema es compatible determinado. 


Si se dan estas condiciones, las dos rectas correspondientes a las dos ecua- 
ciones se cortan, siendo las coordenadas del punto de intersección la solu- 
ción del sistema. 


Rectas secantes 


2. SD=0,E=0yF =0, los cocientes S y - son indeterminados. En 


consecuencia, el sistema tiene infinitas soluciones; o sea, el sistema es 
compatible indeterminado. Esto hace que las dos rectas sean coinci- 
dentes. 


Rectas coincidentes 


63 


3. SD=0,E=00F=0, el cociente = o el cociente + no existe. En 


consecuencia, el sistema no tiene solución; o sea, el sistema es incompati- 
ble. Se tiene, entonces, que las dos rectas son paralelas. 


bd 


Rectas paralelas 


El programa 


El siguiente programa resuelve gráficamente cualquier sistema de dos ecuaciones 
de primer grado con dos incógnitas, siempre que las rectas que correspondan a 
dichas ecuaciones caigan en los límites de la pantalla. 


Una vez resuelto el sistema gráficamente, el programa analiza los valores de E, F y 
D para discutir el sistema, y si éste resulta ser compatible determinado, procede a 
Calcular la solución. Si, en cambio, el sistema es compatible indeterminado o in- 
compatible, da el mensaje correspondiente. 


10 REM SISTEMA DE ECUACIONES 

20 DEF FNF(X) = A * X + B 

25 DEF FNR(X) = INT(100 * X + 0.5)/100 

30 SCREEN 2 

40 OPEN “GRP:” AS 4 1 

50 GOSUB 2000 

60 REM RESOLUCION GRAFICA 

70 FORI|=1T02 

80 READP, Q,R 

85 IF | = 1 THEN PRESET(80, 0) : PRINT 4 1, P; “X+"; Q; “Y="; 

R: GOTO 90 

87 PRESET(80, 10) : PRINT 4 1, P: “X+"; Q; “Y=";R 
90 IFQ = O THEN X = R/P: GOTO 113 
100 A= —P/Q : B =R/Q : GOSUB 3000 
110 GOTO 115 
113 GOSUB 4000 
115 FORT = 1 TO 1000 : NEXT T 
120 NEXT | 
130- DATA. cs coop snep sas 
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— La línea 50 transfiere el control a la subrutina 2000 que dibuja los ejes coor- 
denados. 


— Las líneas 70-120 dibujan las dos rectas. (La línea 100 transfiere el control pre- 
cisamente a la subrutina 3000 que dibuja la gráfica de cada ecuación.) 


— La línea 90 contempla el caso particular en que el coeficiente de Y es cero, en 
cuyo caso la gráfica de la ecuación es una recta vertical, o sea paralela al eje 
de ordenadas (la gráfica la dibuja la subrutina 4000). 


140 RESTORE 

150 REM RESOLUCION ALGEBRAICA 

160 FORI|=1T0 2 

170 READP, Q,R 

180 IFI= 1 THENP1=P:Q1=Q:R1=R: GOTO 200 

190 P2=P:Q02=0Q:R2=R 

200 NEXT | 

210 E=R1*Q2-R2*Q1:F=P1*R2-P2x*R1:D= 
= P1 * Q2 - P2* Q1 

220 IFD=0AND(E = 0 AND F = 0) THEN GOTO 370 

230 IFD=0AND(E<> 0 OR F <> 0) THEN GOTO 380 

240 X=E/D:Y =F/D 


— Con las instrucciones 140-200 se vuelven a leer los coeficientes de las dos 
ecuaciones. 


— La 210 calcula el valor de E, F y D y las líneas 220 y 230 discuten el sistema. 


— La línea 240 calcula X e Y (soluciones), en el caso de que el sistema sea com- 
patible determinado. 


250 IFX<-12 OR X > 12 OR Y > 7 OR Y < —7 THEN GOTO 360 
255 IF Q1=00R Q2 = 0 THEN GOTO 305 

260 H= 100 — Y/0.1 

270 IFH <= 100 THEN $ = 2 : GOTO 280 

275 IFH> 100 THEN S = -3 

280 FOR |= HTO 100 STEP S 

290 PSET (129 + X/0.1, I) 

300 NEXT | 


— En la línea 250 se averigua si las rectas se cortan en los límites de la pantalla. 
En caso afirmativo, el bloque de instrucciones 260-300 traza una línea pun- 
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teada desde el punto de corte hasta el eje de abscisas; así se visualiza el valor 
de la incógnita X. 


Si una de las rectas es vertical (hecho que sucede si Q1 = 0 o Q2 = 0) no se 
dibuja la línea punteada, transfiriéndose el control a la línea 305. 


305 IFP1 = 0 OR P2 = 0 THEN GOTO 360 
310 H= 129 + X/0.1 

320 IFH<= 129 THEN S = 3 : GOTO 340 
325 IFH> 129 THEN S = -3 

330 FOR !|= HTO 129 STEP S 

340 PSET (I, 100 — Y/O.1) 

350 NEXT | 


Estas instrucciones trazan una línea punteada desde el punto de corte hasta el eje 
de ordenadas; así se visualiza el valor de la incógnita Y. 


Si una de las rectas es horizontal (hecho que sucede si P1 = 0 0 P2 = 0 no se 
dibuja la línea punteada, transfiriéndose el control a la línea 360. 


360 PRESET(O, 183) : PRINT +, “SIST. COMPATIBLE DETERMI- 
NADO”: GOTO 390 

370 PRESET (0, 183) : PRINT 4 1, ”S. INDETERM. (RECT. COINCI- 
DENTES)” : GOTO 390 

380 PRESET(O, 183) : PRINT 4 1, “S. INCOMP. (RECT. PARALELAS)” 

390 GOTO 390 

400 END 


2000 REM SUBRUTINA EJES 

2010 REM EJE DE ABSCISAS 

2020 LINE (0, 100) — (255, 100) 

2030 FOR |= 7 TO 247 STEP 10 

2040 PRESET (1, 97) : PRINT 4 1, ”'” 
2050 NEXT | 

2060 PRESET(138, 106) : PRINT + 1, “1” 
2070 PRESET(247, 106) : PRINT + 1, “X” 
2080 REM EJE DE ORDENADAS 

2090 LINE (129, 20) — (129, 180) 

2100 FORI= 24 T0 170 STEP 10 

2110 PRESET (127, |) : PRINT4 1,” ” 
2120 NEXT | 

2130 PRESET(120, 20) : PRINT + 1, “Y” 
2140 PRESET(120, 86) : PRINT + 1, 1” 
2150 RETURN 
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Esta subrutina dibuja los ejes coordenados con divisiones, señalando con X el eje 
de abscisas y con Y el de ordenadas. 


3000 REM DIBUJA LA GRAFICA 

3010 IFA = O THEN X1 = -12 : X2 = 12 : GOTO 3060 
3020 X1 = (-7 — BJ/A: X2 = (7 — BJ/A 

3030 IFA<OTHENC=X1: X1 = X2: X2 = C 
3040 IF X1 < -12 THEN X1 = -12 

3050 IF X2 > 12 THEN X2 = 12 

3060 FOR X = X1 TO X2 STEP 0.1 

3070 Y = FNF(X) 

3080 X3 = INT (X/0.1 + 0.5) 

3090 Y3 = INT (Y/0.1 + 0.5) 

3100 PSET(129 + X3, 100 + Y3) 

3110 NEXT X 

3120 RETURN 


4000 REM RECTAS VERTICALES 
4010 LINE(129 + X/0.1, 20) — (129 + X/0.1, 180) 
4020 RETURN 


Ejecuciones del programa 


e Para resolver el sistema 
—0.5x + y = 3 
2x-y=0 
la instrucción 130 habrá que describirla así: 
130 DATA -0.5, 1,3,2, -1,0 


Al ejecutar el programa con esta instrucción se obtiene este resultado. 
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e Para resolver el sistema 


0.8x- y= 5 
1.6x — 2y = 10 


la instrucción 130 habrá que escribirla así: 
130 DATA 0.8, -1, 5, 1.6, -2, 10 
Al ejecutar el programa se obtiene este resultado. 


e Para resolver el sistema 
E. 
x=y=5 
la instrución 130 habrá que escribirla así: 
130 DATA 1, -1,1,1,-1,5 
Al ejecutar el programa se obtiene este resultado. 
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e Para resolver el sistema 
| 2x + 0y = 8 
0x - 3y=6 
la instrucción 130 hay que escribirla así: 
130 DATA 2, 0, 48, 0, 3, 6 
Al ejecutar el programa se obtiene este resultado 


3. Representación gráfica de la función cuadrática 
f(x) = (mx + n)* + p. (Familia de parábolas) 


Presentación del problema 


La función cuadrática f(x) = ax? + bx + c se puede expresar en la forma 
f(x) = (mx + n)? + p 
cuya gráfica es una parábola. 


Los parámetros m, n y p determinan la forma de la parábola (si es más abierta o 
cerrada) y su posición. 


El punto característico de la parábola es su vértice V, cuyas coordenadas son 
(-— p). 
m 
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La recta vertical que pasa por el vértice de la parábola es su eje de simetría y su 

> n 

ecuación es X= = —p 

Sin = 0 y p = 0 la ecuación de la parábola es f(x) = (mx)? = m? — x? cuyo 

vértice coincide con el origen de coordenadas: | - Al ) (0,0). En consecuencia, 
S. 


el eje de simetría de esta parábola es el eje de ordenada 


leje de 

Í simetría 
on 

IG — — y) 


f(x) = m?x? 


Xx 


vértice 4](0,0) 
eje de simetría 


(x = 0) 


La recta vertical que pasa por el vértice de la parábola es su eje de simetría, y su 


., En TE 
ecuación es X = m* 


El programa 


El siguiente programa dibuja gráficos de funciones cuadráticas. También permite 
obtener en la pantalla familia de parábolas, variando los parámetros m, n y p. El 
número máximo de parábolas que pueden aparecer en la pantalla es 4. Se ha 
puesto este tope porque con un número mayor la pantalla puede aparecer con- 
fusa. 


10 REM PARABOLAS 

20 DEF FNF(X) = (Mx X + N) A 2+P 
30 READ G 
50 SCREEN 2 

60 OPEN “GRP:” AS 4 1 

70 GOSUB 2000 


La instrucción 30 READ G (y su correspondiente DATA) lee el dato G que da el 
número de parábolas que se van a representar. 


La instrucción 70 transfiere el control a la subrutina 2000 que dibuja los ejes 
coordenados. 
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80 FOR!I=1TOG 

90 READ M, N, P 

100 IF | — 1 THEN PRESET (10,0) : PRINT 4 1, "Y = (%; M; “X+"; N;”) 
> 2+”;P: GOTO 120 

110 1F1= 2 THEN PRESET (10,0) : PRINT 4 1, “Y = (; M; “X+"; N;”) 
> 2+";-P: GOTO 120 

115 PRESET (10,20) : PRINT 4 1, "Y =(";M;“X+“;N;“A2+”;P 

120 GOSUB 3000 

130 FORT = 1 TO 1000 : NEXT T 

140 NEXT | 


Las líneas 100 y 110 escriben en la parte superior de la pantalla de alta resolución 
la ecuación de la parábola que se va a representar. 


La instrucción 120 transfiere el control a la subrutina 3000 para dibujar la gráfica 
correspondiente. 


La línea 130 introduce un bucle de retardo de unos dos segundos antes de dibu- 
jarse la siguiente parábola: 


150 PRESET (10, 183) : PRINT + 1, “Para salir, CTRL y STOP” 
160 GOTO 160 


La línea 150 escribe en la parte inferior de la pantalla de alta resolución el texto 
que indica cómo salir del bucle sin fin de la línea 160. 


170 DATA 
180 DATA 
190 END 


En 170 DATA se pondrá el dato correspondiente al número de parábolas que se 
desean representar (G) y en 180 DATA, el mismo número de ternas de datos co- 
rrespondientes a M, N y P. 


2000 REM SUBRUTINA EJES 

2010 REM EJE DE ABSCISAS 

2020 LINE (0, 110) — (255, 110) 

2030 FOR | =7TO 247 STEP 10 

2040 PRESET(I, 107) : PRINT 4 1, “|” 
2050 NEXT | 

2060 PRESET(138, 116) : PRINT 4 1, "1" 
2070 PRESET(247, 116): PRINT 4 1, "X" 
2080 REM EJE DE ORDENADAS 

2090 LINE (129, 30) — (129, 182) 

2100 FOR|= 34TO 180 STEP 10 

2110 PRESET (127, 1) : PRINT 4 1,“ "” 
2120 NEXT | 

2130 PRESET(120, 30) : PRINT 4 1, “Y” 
2140 PRESET(120, 96) : PRINT + 1, “1" 
2150 RETURN 
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Esta subrutina dibuja los ejes coordenados con divisiones, señalando con X el eje 
de abscisas y con Y el eje de coordenadas, y poniendo un 1 en la primera división 
del eje X y del eje Y. 


3000 REM DIBUJA LA GRAFICA 

3010 E =-N/M : IFE<-12 OR E > 12 THEN GOTO 3160 

3015 IF. 7 — P<O THEN PRESET (10,30) : PRINT 4 1, “Si p > 7 no se 
puede repres.” : GOTO 3160 

3020 X1 = ((-1) x SQR (7 — P) — N)/M : X2 = (SQR(7 — P) — N)/M 

3030 IF X1 >X2 THEN C = X1: X1 = X2: X2 = C 

3040 IF X1 <-12 THEN X1 = -12 

3050 IF X2 > 12 THEN X2 = 12 

3070 FOR X = X1 TO X2 STEP 0.05 

3080 Y = FNF(X) 

3090 X3 = INT(X/0.1 + 0.5) : Y3 = INT (Y/0.1 + 0.5) 

3095 IF 110 — Y 3 > 180 THEN GOTO 3110 

3100 PSET(129 + X3, 110 — Y3) 

3110 NEXT X 

3120 REM EJE DE SIMETRIA 

3130 FOR J= 30T0 182 

3135 D=J/6 : IF D = INT(D) THEN GOTO 3150 

3140 PSET (129 + E/0.1, J) 

3150 NEXT J 

3160 RETURN 


La línea 3010 calcula la abscisa del vértice de la parábola y, además, averigua si 
esta abscisa queda fuera de la pantalla, a la izquierda (E < -12) o a la derecha 
(E > 12). Si esto ocurre se decide no dibujar la gráfica, transfiriendo el control a 
la línea 3160. La línea 3015 evita que la ejecución se bloquee en la línea 3020 
cuando 7 — P < O, ya que la raíz cuadrada de un número negativo no existe en 
el conjunto de los números reales. 


El bloque de instrucciones 3020-3050 determina el intervalo [X1, X2]. 


El bloque 3070-3110 representa la parábola. 

Las instrucciones 3130-3150 dibujan el eje de simetría. La línea 3095 hace que r no 
se representen en la pantalla los A UIOS de la parábola cuyas ordenadas sean me- 
nores que —7, respecto al origen de coordenadas elegido: (129, 110). 


Ejecución del programa 
e Si en las instrucciones 170 y 180 se ponen los datos que figuran a continua- 
ción: 


170 DATA 2 
180 DATA 1, 0,0, 1, 0, 5 


Ya 


se obtiene este resultado: 


e Las ecuaciones siguientes corresponden a 3 parábolas pertenecientes a la 
misma familia (el parámetro que varía es n). 


y = (0.8x + 0)? — 2 = (0.8x)? - 2 = 0.64x? — 2 
y = (0.8x + 1) - 2 
y = (0.8x — 4.3) - 2 


Para dibujar estas parábolas, las instrucciones 170 y 180 deberán escribirse así: 


170 DATA 3 
180 DATA 0.8, 0, 2, 0.8, 1, -2, 0.8, 4.3, -2 


e Las ecuaciones siguientes corresponden a 3 parábolas de la misma familia (el 
parámetro que varía es m): 


y = (x +4) -5 
y = (2x + 4) -5 
y = (0.5x + 4)? - 5 


Para dibujar estas parábolas las instrucciones 170 y 180 deberán escribirse así: 


170 DATA 3 
180 DATA 1, 4, 5, 2, 4, -5, 0.5, 4, -5 


4. Resolución gráfica y algebraica de la ecuación de segundo 
grado 


Gráficamente, las soluciones x, y Xx, de la ecuación de 2.* grado 
axe +bx+c=0 

vienen dadas por los puntos de corte de la parábola 
y = ax +bx+c 
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con el eje de abscisas: 


X1 xX2 X 


Es evidente que si la parábola no corta al eje de abscisas la ecuación ax? + bx + 
+ c = 0 no tiene soluciones reales. 


e Algebraicamente, las soluciones de la ecuación ax? + bx + c = 0 se obtienen 
mediante las conocidas fórmulas 


_ —b+ Vb? - ac -  —b- vb?-4ac 
ES” E: NE. ES” 


El radicando b* — 4ac que aparece en estas fórmulas se llama discriminante, 
D. Este valor determina el tipo de solución de la ecuación de 2.* grado: 


— Si D >0, se tiene que VD > 0 y, en consecuencia, las soluciones x; y X, 
serán reales y distintas. 


— Si D = 0 se tiene que VD = 0 y, en consecuencia, las soluciones x, y X> 
serán reales y coincidentes; es decir, la ecuación tendrá una sola solución 
real. 


— Si D <0, se tiene que VD no es un valor real (es imaginario) y, en conse- 
cuencia, la ecuación no tiene soluciones reales. 


El programa 


El siguiente programa resuelve gráficamente cualquier ecuación ax? + bx + c= 
= ( siempre que tenga, soluciones reales y siempre que la parábola de la función 
correspondiente y = ax? + bx + c entre en los límites de la pantalla. 


Una vez resuelta la ecuación gráficamente, el programa la discute, y si tiene solu- 
ciones reales procede a calcularlas. Si, en cambio, las soluciones no son reales, da 
el mensaje correspondiente. 
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10 REM ECUACION DE SEGUNDO GRADO 

20 DEF FNF(X) = (M*X+N)A2+P 

30 DEF FNR(X) = INT(100 * X + 0.5)/100 

40 SCREEN 2 

50 OPEN “GRP:” AS 4 1 

60 GOSUB 2000 

70 REM SOLUCION GRAFICA 

80 READA, B, € 

85 PRESET(10,0) : PRINT 4 1, A; “XA 2+"; B; “X+"; C; “=0" 
90 GOSUB 3000 
100 DATA 
110 REM RESOLUCION ALGEBRAICA 
120 D=B*B-4x*Ax*C 

130 IFD< O THEN GOTO 200 

140 X1 = (-B + SQR(D)/(2 * A) : X2 = (B — alió * A) 
150 IF X1 > X2 THEN L = X1 : X1 = X2 : X2 = 

160 IF X1 <-12 OR X2 > 12 THEN GOTO 190 

170 PRESET(129 + X1/0.1 — 14, 106) : PRINT 4 1, “X1” 

180 PRESET(129 + X2/0.1 + 5, 106) : PRINT 4 1, “X2” 

190 PRESET(10,10) : PRINT 4 1, “X1="; FNR(X1); “X2="; FNR(X2) : 

GOTO 210 

200 PRESET(10,10) : PRINT 4 1, “NO TIENE SOLUCIONES REALES” 
210 PRESET(10, 183) : PRINT 4 1, “Para salir, CTRL y STOP.” 
220 GOTO 220 
230 END 


Este bloque de instrucciones constituye el programa principal. 
— La instrucción 90 transfiere el control a la subrutina 3000 que dibuja la gráfica 
de la función y = ax? + bx + c. 


— La instrucción 120 calcula el discriminante D y la 130 averigua si la ecuación 
no tiene soluciones reales. 


— La línea 140 calcula las soluciones. 


— Las líneas 170 y 180 escriben en la pantalla de alta resolución los símbolos X1 
y X2 junto a los puntos de corte de la parábola con el eje de abscisas, siempre 
que ello sea posible (fijarse en la línea 150). 


— La línea 190 escribe en la parte superior de la pantalla las soluciones de la 
ecuación, siempre que sean reales, y la línea 200 escribe el mensaje NO TIENE 
SOLUCIONES REALES (fijarse en la línea 130). 


— La línea 210 escribe, en la parte inferior de la pantalla cómo proceder para sa- 
lir del bucle sin fin de la línea 220. 
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— La instrucción 230 da por terminado el programa. 


2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 


REM SUBRUTINA EJES 
REM EJE DE ABSCISAS 

LINE (0, 100) — (255, 100) 

FOR | = 7 TO 247 STEP 10 
PRESET(I, 97) : PRINT 4 1, “|” 
NEXT! 

PRESET(138, 106) : PRINT 4 1, 1" 
PRESET(247, 106) : PRINT 4 1, “X" 
REM EJE DE ORDENADAS 

LINE (129, 20) — (129, 180) 

FOR | = 24 TO 170 STEP 10 
PRESET(127, 1) : PRINT4 1," " 
NEXT | 

PRESET(120, 20) : PRINT 4 , Ay 
PRESET(120, 86) : PRINT 4 1, "1" 
RETURN 


Esta subrutina dibuja los ejes coordenados con divisiones, señalando con X el eje 
de abscisas y con Y el de ordenadas, y poniendo un 1 en la primera división del 
eje X y del eje Y. 


3000 
3010 
3020 
3030 
3040 
3045 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 


REM DIBUJA LA GRAFICA y 

IFA < O THEN Al = -A : B1 =-B : C1 =—C : GOTO 3030 
Al=A:B1=B:C1 =C 

M = SQR(A1) : N = B1/(2* M):P=C1-Nx*N 

E = —-N/M : IFE <-12 OR E > 12 THEN GOTO 3160 

IF7 — P<O THEN GOTO 3160 

X4 = ((21) * SQR (7 — P) — N)/M : X5 = e (7 — P) — N/M 
IFX4 > X5 THEN C=X4:X4=X5: X5 = 

IF X4 < -12 THEN X4 = -12 

IFX5 > 12 THEN X5 = 12 

FOR X = X4 TO X5 STEP 0.05 

Y = FNF(X) 

IFA < 0 THEN Y = -Y 

X3 = INT (X/0.1 + 0.5) : Y3 = INT (Y/0.1 + 0.5) 

IF 100 — Y3 <30 OR 100 — Y3 > 170 THEN GOTO 3150 
PSET (129 + X3, 100 — Y3) 

NEXT X 

RETURN 


— En esta subrutina, primero se averigua si el coeficiente A es negativo, y en 
caso afirmativo se definen las variables A1, B1 y C1 en las que se almacenan 
los valores de A, B y C, respectivamente, pero cambiados de signo (línea 
3010). Hacer esto equivale «a dar la vuelta» a la parábola de coeficientes A, B 


y C. 
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— Si el valor de A es positivo, en las variables A1, B1 y C1 se almacenan los 
mismos valores de A, B y C, respectivamente (línea 3020). 


En la línea 3030 se calculan los valores de M, N y P en función de los valores 
de A1, B1 y C1. De este modo, la función y = ax? + bx + c se pone en la 
forma y = (mx + n)? + p, y así se pueden aprovechar la mayoría de las ins- 
trucciones de la subrutina 3000 del programa anterior (representación gráfica 
de la función cuadrática f(x) = (mx + n)? + p. 

— El bloque de instrucciones 3040-3130 dibujan la parábola. 

— La instrucción 3110 cambia de signo el valor de Y en el caso de que el coefi- 
ciente A de x* sea negativo. Esto se hace para que en la pantalla aparezca la 
verdadera parábola, ya que las asignaciones hechas en la línea 3010 producen 
el efecto de “dar la vuelta” a la misma. 

La línea 3130 hace que no se representen en la pantalla los puntos de la pará- 


bola cuyas ordenadas sean mayores que 7 O menores que —7, respecto al or- 
den de coordenadas elegido: (129,100). 


Ejecuciones del programa 


e Para hallar las soluciones de la ecuación x? + 2x-8=0 
la instrucción 100 hay que escribirla así 
100 DATA 1, 2, 8 


Al ejecutar el programa se obtiene este resultado: 


E e 1 
e Para hallar las soluciones de la ecuación — 7 x? + 4x-24=0 


la instrucción 100 hay que escribirla así 
100 DATA 0.5, 4, -24 
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Al ejecutar el programa se obtiene este resultado: 


e Para hallar las soluciones de la ecuación. 
2+x+1=0 
la instrucción 100 hay que escribirla así: 
100 DATA 1,1, 1 
Al ejecutar el programa se obtiene este resultado: 


5. Matemática básica 


1. Las tablas de multiplicar 
Presentación del problema 


Se trata de practicar las tablas de multiplicar. 


El programa 


Se generan dos números al azar de una cifra y se pregunta cuál es su producto. 
Después de comparar la contestación con el resultado correcto, se anotan los 
aciertos y, en caso de error, se dan dos nuevas ocasiones. Se proponen 10 nú- 
meros utilizando un bucle (Il = 1 a 10). 


1 REM MULTIPLICAR 
10 Z = RND (-TIME) 
20 C=0:D=0 
30 FOR!=1T0 10 
40 A= INT(RND(1) * 9 + 1) 
50 B= INT(RND(1) * 9 + 1) 


60 P=A*B 

70 PRINTA;**";B;“=" 
80 INPUT X 

90 IF X =P THEN GOTO 150 
100 D=D+1 


110 PRINT “LO SIENTO: TE EQUIVOCASTE” 

120 IFD> 2 THEN GOTO 160 , 

130 PRINT “INTENTALO DE NUEVO” 

140 GOTO 70 

50 C=t+1 

160 PRINT.A; “+: B: *=": Pp 

170 NEXT | 

180 PRINT “EL NUMERO DE ACIERTOS ES”; C 
190 END 


La instrucción 10 hace que los números que generan las instrucciones 40 y 50 
sean aleatorios. Eliminar la instrucción 10 si el microordenador no es MSX; en 
todo caso, sustituirla por la que corresponda al microordenador que se utilice. 
Con la instrucción 20 se inicializan los dos contadores C y D. 
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En las instrucciones 40, 50 y 60 se generan los números de una cifra A y B al 
azar, multiplicándolos a continuación y almacenando el resultado en P. 


Se pide una contestación mediante la instrucción 80 y se compara en 90 con la 
solución correcta. Si la contestación es acertada, se suma 1 al contador C, en la 
150, y se escribe el resultado correcto. Si es errónea, se suma 1 al contador D, en 
la 100; según el valor de D, se dan nuevas oportunidades (desde D = 0 hasta 
D = 2). El proceso se repite 10 veces. 


Finalmente, en 180 se imprime el contenido final del contador C, el cual da el nú- 
mero de aciertos. 


2. Resolución algebraica de la ecuación de segundo grado 


Presentación del programa 


Una ecuación de segundo grado Ax? + Bx + C = 0 tiene siempre dos solu- 
ciones, X; y X,, dentro del cuerpo C de los números complejos. En el programa si- 
guiente se van a obtener estas dos soluciones, indicando si son reales o com- 
plejas. Primeramente se introducen como datos los coeficientes A, B y C. 


El programa 


1 REM SOLUCION DE UNA ECUACION DE SEGUNDO GRADO 
10 INPUT A, B, C 
20 D=B 1 24*A*C 
30 IF D< O THEN GOTO 90 


En la instrucción 10 se leen los coeficientes A, B y C de la ecuación. La 20 asigna 
a D el valor del discriminante B* — 4AC. Por último, la instrucción 30 pregunta si 
el discriminante es negativo y, en tal caso, envía el control a la instrucción 90, 
para calcular las soluciones complejas. Si no se verifica que D < O, las raíces son 
reales y su cálculo comienza en la instrucción 40. 


40 X1 = (B + SQR(D))/(2 * A) 

50 X2 = (8 — SQR(D))/(2 * A) 

60 PRINT “LAS SOLUCIONES SON REALES Y SUS” 
70 PRINT “VALORES SON:”; X1; “Y”; X2 

80 END 


En las instrucciones 40 y 50 se calculan las dos soluciones reales X1 y X2 (que son 


iguales cuando D = 0). En la 70 se escriben las soluciones y en la 80 se para el 
programa. 
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90 PR = -B/(2 * A) 

100 PY = SQR(-DJ/(2 * A) : PY = ABS(PY) 

110 PRINT “LAS SOLUCIONES SON COMPLEJAS Y” 

120 PRINT “SUS VALORES SON:"; PR; “+”; Pl; *1%; “Y%; PR; “-”; 
PY; Al 


La instrucción 90 calcula la parte real PR y la 100 la parte imaginaria PY. Como 
los dos números complejos son conjugados, basta con cambiar el signo de la 
parte imaginaria para obtener las soluciones en la instrucción 120. 


Ejecuciones del programa 


e Si se ejecuta el programa, para A = 1,B =-4yC = 3, se obtiene este resul- 
tado en la pantalla: 


LAS SOLUCIONES SON REALES Y SUS 
VALORES SON: 3 Y 1 


Si se ejecuta el programa, para A = 1, B =-2, € = 10, se obtienen dos solu- 
ciones complejas conjugadas: 


LAS SOLUCIONES SON COMPLEJAS Y 
SUS VALORES SON: 1 + 31 Y 1 — 31. 
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3. Producto de polinomios 
Presentación del programa 


Para multiplicar dos polinomios P(x) y Q(x) los ordenamos previamente de mayor a 
menor grado. 

PO):= ap, 1? + ax + + ax + as 

QU) = byyx! + byx07? 4... + b2x + by 


A continuación multiplicamos todos los monomios de P(x) por cada monomio de 
QU), tal como se muestra en este ejemplo: 


P(x) = 2x7 + x2+3x- 1y 
QU) = x? —.4x + 1 


Z 1 3 - 
14 1 
2 1. 3 - 
3 4 12 4 
2.1.3 4 
2-1. 1412 7 - 


Luego el polinomio producto es: 
R(x) = 2x7 - 7xé + 2 — 12x? + 7x - 1 
Para escribir un programa que realice este producto se necesitan las siguientes va- 
riables: 
P: grado del polinomio P(x) 
A(P + 1): lista para guardar los P + 1 coeficientes de P(x) 
Q: grado del polinomio Q(x) 
B(Q + 1): lista para guardar los Q + 1 coeficientes de Q(x) 


C(Q + 1,P + Q + 1): tabla para guardar las Q + 1 filas que se forman en los 
pasos intermedios de la multiplicación. Esta tabla tiene P + Q + 1 columnas 
pues este es el número de coeficientes del polinomio, producto R(x). 


R(P + Q + 1): lista para guardar los P + Q + 1 coeficientes del polinomio pro- 
ducto R(x). 


En el ejemplo anterior, tenemos: 


3  —A(21,3,-1) 
2 


P= 
Q = B (1, 4, 1) 
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2 1 3 - 
14 1 
2 1-3. 2 
3 4 12 4 Q + 1 filas (2 + 1 =3) 
2.1.3 - 


2-1. 1-12 7-4 


P+ Q + 1 columnas (3 + 2 + 1 = 6) 


El programa 


INPUT “GRADO DEL POLINOMIO P”; P 

DIM A(P + 1) 

FOR | =P + 1T0 1 STEP — 1 

PRINT “COEFICIENTE DEL TERMINO DE GRADO”; | — 1 
INPUT A(l) 

NEXT | 

INPUT “GRADO DEL POLINOMIO Q”; Q 

DIM B(Q + 1) 

FOR | = Q + 1T0 1 STEP — 1 

PRINT “COEFICIENTE DEL TERMINO DE GRADO”; | — 1 
INPUT B(1) 

NEXT | 

CLS 


En el bloque de instrucciones 10-60 se introduce el grado del polinomio P(x) y sus 
coeficientes, de mayor grado a menor. 


En el bloque de instrucciones 70-120 se realiza lo mismo que en el bloque ante- 
rior, pero en este caso para el polinomio Q(x). 


140 
150 
160 
170 
180 
190 


DIMC(Q+1,P+Q+1) 
FORI|=1T0Q +1 
FORJ=1TOP +1 
C(l,1+J-— 1) = B(I) * AQ) 
NEXT J 

NEXT | 


En este bloque se calculan los pasos intermedios de la multiplicación (la tabla C). 
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200 DIMR(P +Q + 1) 
= 1 


210 FORK=1TOP+0Q+1 
220 FORI=1TOQ+ 1 
230 R(K) = R(K) + C(l, K) 
240 NEXT | 

250 NEXT K 


Con estas instrucciones se calculan los coeficientes del polinomio producto R(x). 


260 FORK =P+0Q+1T0 1 STEP — 1 

270 PRINT “COEFICIENTE DEL TERMINO DE GRADO”; K — 1 
280 PRINT R(K) 

290 NEXT K 

300 END 


Este conjunto de instrucciones imprimen los coeficientes del polinomio producto. 


Ejecución del programa 


Si se ejecuta el programa para los polinomios: P(x) = 2x? + x? + 3x — 1 
y Q0)= x?-4x+1 
se obtiene el siguiente resultado en la pantalla: 


run 

GRADO DEL POLINOMIO P? 3 
COEFICIENTE DEL TERMINO DE GRADO 3 
22 

COEFICIENTE DEL TERMINO DE GRADO 2 
21 


COEFICIENTE DEL TERMINO DE GRADO 1 
? 


COEFICIENTE DEL TERMINO DE GRADO 0 
GRADO DEL POLINOMIO Q?2 

COEFICIENTE DEL TERMINO DE GRADO 2 
COEFICIENTE DEL TERMINO DE GRADO 1 
COEICENTE DEL TERMINO DE GRADO 0 
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Al seguir la ejecución se borra la pantalla anterior y se obtiene el resultado si- 
guiente: 


EME DEL TERMINO DE GRADO 5 
COERCIENTE DEL TERMINO DE GRADO 4 
COERCINTE DEL TERMINO DE GRADO 3 
COEFICIENTE DEL TERMINO DE GRADO 2 


2-12 

COEFICIENTE DEL TERMINO DE GRADO 1 
27 

COEFICIENTE DEL TERMINO DE GRADO 0 
2-1 

OK 

E 


4. Formación de las variaciones con repetición 
Presentación del problema 


e Las variaciones con repetición de M elementos, tomados de N en N, son 
los distintos grupos que se pueden formar de N elementos repetidos o 
no, tales que dos cualesquiera de ellos difieren en la naturaleza de algún 
elemento o en el orden de colocación de los mismos. 


Por ejemplo, con los elementos del conjunto (A, B, C) las variaciones con 
repetición de orden 3 que se pueden formar son las siguientes (en este 
ejemplo M = 3 y N = 3). 


M =3 

N =3 

AAA BAA CAA 
AAB BAB CAB 
AAC BAC CAC 
ABA BBA CBA 
ABB BBB CBB 
ABC BBC CBC 
ACA BCA CCA 
ACB BCB CCB 
ACC Bcc CEE 


e El número total de variaciones con repetición, en este caso es: 
VR33 == 27 
En general: VRyy = MP 


e La formación de las variaciones con repetición es independiente de los 
elementos. Si en el ejemplo anterior los elementos hubieran sido los del 
conjunto (1, 2, 3), bastaría sustituir A por 1, B por 2 y C por 3. 


Luego a la hora de hacer el programa lo importante es el valor de M y N, 
pues al primer elemento se le asigna el 1, al segundo elemento se le 
asigna el 2, ..., y al n-ésimo elemento se le asigna N. 


e El proceso a seguir para formar las variaciones con repetición es el que * 
describimos a continuación (tomamos como ejemplo M = 4 y N = 3), 


1.2 Se calcula el número de variaciones con repetición V = MY, 
V=43 =64 


2.” Se hace variar | de 1 a V, para indicar la variación correspondiente. 


Haremos x = | — 1. 
Si dividimos x entre N1W* = 4? y a su parte entera le sumamos 1, 
tenemos el 1.% elemento: B(3). 


Si el resto lo dividimos entre MW? = 4 y a su parte entera le 
sumamos 1, tenemos el 2.* elemento: B(2). 


Si al resto le sumamos 1, tenemos el 3. elemento: B(1). 


Como vemos, se van formando las variaciones con repetición. 
3.2 Para imprimir se realiza esta asignación: 


B(3) 1. elemento 
B(2) 2.2 elemento 
B(1) 3." elemento 
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El programa 


10 INPUT “INTRODUCE LOS CARACTERES QUE DESEES VARIAR”; A$ 
20 M= LEN(AS$) 

30 INPUT “ORDEN DE LAS VARIACIONES”; N 

40 DIM B(N) 

50 V=MIN 


La línea 10 guarda en A$ los elementos con los que se desean formar las varia- 
ciones con repetición, y la 20, el número M de dichos elementos. 


La línea 30 pide el orden N de las variaciones. 
La línea 40 dimensiona el vector B. 


La línea 50 calcula el número de variaciones. 


60 FORI=1TOV 
70 X=1-1 

80 IF N = 1 THEN GOTO 140 
90 FORJ=N- 1T0 1 STEP - 1 
100 Z=XIM ? )) 

110 B() + 1) = INT(Z) + 1 

120 X= X-INT(Z) x (M 1 )) 

130 NEXT] 

140 B(1)=X+1 


En este bloque de instrucciones se obtienen los valores de B(N), B(N-1), ..., B(1) 
que determinan los elementos de la variación |. 


150 FOR K = NTO 1 STEP — 1 
160 PRINT MID$(A$, B(K), 1); 
170 NEXT K 
180 PRINT 

- 190 NEXT | 
200 END 


En este bucle se imprimen los elementos de la variación 1. 
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Ejecuciones del programa 


Si se introducen como datos los elementos (A, B) y se hace N = 3 se obtiene 
este resultado: 


run 
INTRODUCE LOS CARACTERES QUE DESEES 
VARIAR ? AB 

ORDEN DE LAS VARIACIONES ? 3 
AAA 

AAB 

ABA 

ABB 

BAA 

BAB 

BBA 

BBB 

OK 

B 


Si se introducen como datos los elementos (R, S, T) y se hace N = 2 se ob- 
tiene este otro resultado: 


run 

INTRODUCE LOS CARACTERES QUE DESEES 
VARIAR ? RST 

ORDEN DE LAS VARIACIONES ?2 

RR 

RS 

RT 

SR 


SS 
ST 
TR 
TS 
TT 
OK 
a 
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5. Formación de las variaciones ordinarias 


Presentación del problema 


Las variaciones ordinarias de M elementos, tomados de N en N, son los distintos 
grupos que se pueden formar de N elementos sin repetir, tales que dos cuales- 
quiera de ellos difieran en la naturaleza de algún elemento o en el orden de colo- 
cación de los mismos. 


Por ejemplo, con los elementos del conjunto (A, B, C, Dj las variaciones ordina- 
rias de orden 3 que se pueden formar son las siguientes (en este ejemplo M = 4 
y N =3). 


ABC BAC CAB DAB 
ABD BAD CAD DAC 
ACA BCA CBA DBA 
ACD BCD CBD DBC 
ADB BDA CDA DCA 
ADC BDC CDB DCB 
El número total de variaciones ordinarias, en este caso, es: 
Vaz = 4.3.2 = 24 
En general: 
Van = MM — 1)... (M-— N + 1) 
El programa 


El programa que forma las variaciones ordinarias se hace a partir del anterior, 
añadiendo algunas instrucciones para eliminar la repetición. 


10. INPUT “INTRODUCE LOS CARACTERES QUE DESEES VARIAR”; A$ 
20 M= LEN(AS) 

30 INPUT “ORDEN DE LAS VARIACIONES”; N 

40 DIM B(N) 

50 V=MIN 


La línea 10 guarda en A$ los elementos con los que se desea formar las varia- 
ciones ordinarias, y la 20 el número M de éstos. 


La línea 30 pide el orden N de las variaciones. 
La línea 40 dimensiona el vector B. 
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La línea 50 calcula el número de variaciones con repetición. 


60 FORI=1TOV 

710 X=1-1 

80 IFN = 1 THEN GOTO 140 
90 FOR J=N- 1T0 1 STEP — 1 


120 x=x-— INT(Z) * (M 7 J) 

122 IFJ=N-— 1 THEN GOTO 130 

124 FORL=J+2TON 

126 IF B(J + 1) = B(L) THEN GOTO 190 
128 NEXT L 

130 NEXT) 

140 B(1) =X+ 1 


En este bloque de instrucciones se obtienen los valores de B(N), B(N — 1), ..., B(1) 
que determinan los elementos de la variación |. 


Las instrucciones 122, 124, 126 y 128 eliminan las repeticiones, saltando a otra 
variación (GOTO 190 de la línea 126). 


142 IFN = 1 THEN GOTO 150 

144 FORL=2TON 

146 IF B(1) = B(L) THEN GOTO 190 
148 NEXT L 


Este bloque de instrucciones comprueba si B(1) es igual a alguno de los otros ele- 
mentos del vector B, y en caso afirmativo elimina esta variación, saltando a la ins- 
trucción 190. 


150 FOR K = NTO 1 STEP — 1 
160 PRINT MID$ (AS, B(K), 1); 
170 NEXT K 

180 PRINT 

190 NEXT | 

200 END 


Este bucle imprime los elementos de la variación |. 
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Ejecución 


e Si se introducen como datos los elementos del conjunto (A, B, C, D) y se hace 
N = 3 se obtiene este resultado. 


run 

INTRODUCE LOS CARACTERES QUE DESEES 
VARIAR? ABCD 

ORDEN DE LAS VARIACIONES ?3 


e Si el número de elementos coincide con el orden (M = N) se obtienen las per- 
mutaciones. Por ejemplo, para (1, 2, 3) y N = 3 se tiene: 


run 

INTRODUCE LOS CARACTERES QUE DESEES 
VARIAR ? 123 

ORDEN DE LAS VARIACIONES ?3 

123 

132 

213 

231 


312 
321 
OK 
5 
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6. Comercio 


1. Tabla de amortización de un préstamo 
Presentación del problema 


Los pagos que hay que realizar anualmente a un Banco o a una entidad finan- 
ciera para amortizar un préstamo es lo que se llama anualidades A de amorti- 
zación. 

La fórmula de A es: 


ss Lali 
(1 +1 -1 
donde: 
es el capital prestado por el Banco 
es el tanto por uno anual 
es el tiempo en años 
es la anualidad que hay que pagar para amortizar el préstamo. 


DP A-=M 


El programa 


El siguiente programa calcula el valor de A y forma una tabla con cuatro co- 
lumnas que indican el número del pago (K), la parte del pago correspondiente al 
interés (X), la parte del pago correspondiente a la amortización (Y), y el capital 
amortizado hasta este pago (2). 


10 INPUT “CAPITAL”; C 

20 INPUT “INTERES (EN % ANUAL)”; | 
30 |=1/100 ] 

40 INPUT “TIEMPO (EN AÑOS)”; T 

50 B=(1+1) 9 1 

60 A=(Cx*1x*B)(B - 1) 

70 A = INT (A + 0.5) 

80 PRINT “LAS ANUALIDADES SON:”; A 


Las líneas 10, 20 y 40 introducen los valores de C, |, T. La línea 30 transforma el 
tanto por ciento anual en tanto por uno. 


Las líneas 50, 60 y 70 calculan el valor de A, redondeándolo sin decimales, y la lí- 
nea 80 imprime su valor. 


90 PRINT ”PAGO” 
100 PRINT TAB(4); “INTERESES” 
110 PRINT TAB(13); “AMORTIZACION” 
120 PRINT TAB(22); "CAPITAL AMORTIZADO” 


Estas cuatro líneas imprimen las cabeceras de las cuatro columnas. 


130 X 
140 F 
150 X 
160 Y 
170 Z=Z2+Y 

180 PRINT K; 

190 PRINT TAB(4); X; 

200 PRINT TAB(13); Y; 
210 PRINT TAB(22); Z 

220 NEXT K 

230 END 


O 
> 


El bucle de variable K calcula e imprime los elementos de la tabla. 


Ejecución del programa 


Si se ejecuta el programa con los siguientes datos: 


Capital prestado = 4 000 000 
Interés anual = 14,> % 
Tiempo en años = 15 


aparece en la pantalla el siguiente resultado: 
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2. Tabla de formación de un capital 


Presentación del programa 

Los pagos que se realizan, al principio de cada año, a un Banco o a una entidad 
financiera para formar un capital es lo que se llama anualidad A de capitalización. 
La fórmula de A es: 


Ñ Cl 
AS ADMIN 


donde: 
C  esel capital que se quiere formar 


Il es el tanto por uno anual 


T_ esel tiempo en años 


A es la anualidad que hay que pagar, al principio de cada año, para formar 
el capital C. 


El programa 


El siguiente programa calcula el valor de A y forma una tabla de tres columnas 
que indican el número de pago (K), los intereses generados (X) por este pago, y 
los pagos anteriores y sus intereses, y el capital acumulado (2) con este pago y los 
intereses. 


10 INPUT “CAPITAL”; C 
20 INPUT “INTERESES (EN % ANUAL)”; | 


30 |=1/100 A 
40 INPUT “TIEMPO (EN AÑOS)”; T 
50 B=1+!) 


60 A=(CxI/(B * (V 1 T- 1) 
70 A= INT(A + 0,5) 
80 PRINT "LAS ANUALIDADES SON:”; A 


Las líneas 10, 20 y 40 introducen los valores de C, 1, T. 

La línea 30 transforma el tanto por ciento anual en tanto por uno. 

Las líneas 50, 60 y 70 calculan el valor de A redondeándolo sin decimales. 
La línea 80 imprime el valor de A. 
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90 PRINT “PAGO” 
100 PRINT TAB(4); “INTERESES” 
110 PRINT TAB(13); “CAPITAL ACUMULADO” 
120 X=0:Z=0 
130 FORK = 1TOT 
140 X = INT ((Z + A) * | + 0.5) 
150 Z=Z+A+X 
160 PRINT K; 
170 PRINT TAB(4); X; 
180 PRINT TAB(13); Z 
190 NEXT K 
200 END 


Las líneas 90, 100 y 110 imprimen las cabeceras de las tres columnas. El bucle de 
variable K calcula e imprime los elementos de la tabla. 


Ejecución del programa 


Si se desea formar un capital de un 1000 000 de pesetas, en 10 años, y el in- 
terés que da el Banco es del 11,5 anual, las anualidades que hay que pagar al 
principio de cada año y la tabla de capitalización son las siguientes: 


run 

CAPITAL? 1 000 000 

INTERES (EN % ANUAL)? 11.5 
TIEMPO (EN AÑOS)? 10 

LAS ANUALIDADES SON: 52356 


Pago Intereses Capital acumulado 


1 6021 58377 
12734 123467 
20220 196043 
28566 276965 
37872 367193 


59818 579971 
72718 705045 
87101 844502 


2 
3 
4 
5 
6 48248 467797 
7 
8 
y 
O 103139 999997 
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3. Elección de la variable en las anualidades de amortización 
y capitalización 


Presentación del problema 


En este programa se puede elegir entre anualidades de amortización (1) y anuali- 
dades de capitalización (2). Además se puede elegir la variable que se desea cal- 
cular: 

(1) Capital C 

(2) Interés | 

(3) Tiempo T 

(4) Anualidad A 


Los cálculos se realizan por medio de las siguientes fórmulas y algoritmos: 


1. Anualidades de amortización 
1. Capital 
ds AIM +1) - 1] 
(131 
2. Interés 


1. Se estima el interés | = 10%, por ejemplo, y se inicia la última estima- 
ción del interés 11 = 0. 


2. Se calcula una estimación de la anualidad 


ar(1+ 1 
A1 = —"— 
(1+ 1-1 
3. Se hace X = ABS((l — 11)/2) 
ell =] 
4. SIA1>Aentonces | =1|-— X 


SIA1 < A entonces | =1|+X 
(en ambos casos se vuelve a 2). 


5. SiA1 = A se imprime el interés | en % anual con dos decimales. 


3. Tiempo 
A 
- LOG (lg) 
LOG (1 + I) 
4. Anualidad 
Cl (1+1 
as (1+ 1-1 
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2. Anualidades de capitalización 


1. Capital 


CAGADO +01 
o | 


2. Interés 


1. Se estima el interés | = 10 %, por ejemplo, y se inicia la última estima- 
ción del interés 11 = 0. 
2. Se calcula una estimación de la anualidad 
o] 
(1+D[0 +1" - 1] 
3. Se hace X = ABS((I — 11)/2) 


A1 = 


ell =| 
4. SIA1>Aentonces | =1| + X 
S1A1 < A entonces | = | — X 


(en ambos casos se vuelve a 2). 
5. SIA1 = A se imprime el interés | en % anual con dos decimales. 


3. Tiempo 
C-1 
E EEN 1+ A El 
y LOG (1 + l) 
4. Anualidad 
Ñ C-| 
A= 7 DIO +01] 
El programa 


10 INPUT “AMORTIZACION (1), CAPITALIZACION (2)”; N 

20 INPUT “VARIABLE: ii (1), INTERES (2), TIEMPO (3) y 
ANUALIDADES (4)"; M 

30 IFM<> 1 THEN INPUT “CAPITAL”; C 

40 IF M<> 2 THEN INPUT “INTERES (EN % ANUAL : | =1/100 

50 IF M <> 3 THEN INPUT “TIEMPO (EN AÑOS)”; T 

60 IF M <> 4 THEN INPUT “ANUALIDADES”; A 

70 IFN = 1 THEN ON M GOSUB 1000, 2000, 3000, 4000 

80 IFN = 2 THEN ON M GOSUB 5000, 6000, 7000, 8000 

90 END 
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Ejecuciones del programa 


Si se desea calcular el interés | que corresponde a un préstamo de 2 000 000 
de pesetas, con anualidades de 400 000 pesetas durante 10 años, basta seguir 
los pasos que muestra la pantalla: 


run 

AMORTIZACION(1), CAPITALIZACION(2) ?1 
VARIABLE: CAPITAL(1), INTERES(2), 
TIEMPO(3) Y ANUALIDADES(4) ?2 
CAPITAL ? 2 000 000 


TIEMPO (EN AÑOS) ? 10 
ANUALIDADES? 400 000 
INTERES (EN % ANUAL) = 15.1 
OK 

3 


Si se introducen los datos: 

(2) anualidades de capitalización, 

(3) tiempo, 

capital C = 3.000.000 

interés | = 16,5 % 

anualidad A = 900.000 

se obtiene el siguiente resultado en la pantalla: 


run 

AMORTIZACION(1), CAPITALIZACION(2) ?2 
VARIABLE : CAPITAL(1), INTERES(2), 
TIEMPO(3) Y ANUALIDADES(4) ?3 
CAPITAL ? 3 000 000 


INTERES (EN % ANUAL) ? 16.5 
ANUALIDADES ? 900 000 
TIEMPO = 2 AÑOS 6 MESES 
OK 

A 


Nota: Realmente al cabo de dos años tendría 2.270.000 ptas. más las 900.000 


pesetas de la anualidad correspondiente; es decir, 3.170.000 ptas. Luego 
el tiempo necesario es el que indica los años, despreciando los meses. 
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7. Estadistica 


1. Media, varianza y desviación típica 


Presentación del problema 
Si se tiene una serie de valores x; y sus frecuencias absolutas n 


La media aritmética, la varianza y la AESNAcIOO típica vienen determinadas por las 


fórmulas siguientes: 


k 
2 nx 
iS |=] 
M n 
k 
2 nx - m)? 
yo di 
N n 
s= VV 


donde n es la suma de todas las frecuencias absolutas n 
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El programa 


10 
20 


INPUT “NUMERO DE VALORES DISTINTOS”; K 
DIM X(K) : DIM N(K) 

FOR | = 1 TO K 

INPUT “VALOR, FRECUENCIA”; X(), N(I) 
NEXT | ' 


En este grupo de instrucciones se leen los distintos valores y sus frecuencias abso- 


lutas. 


M=0:N=0 

FOR | = 1 TO K 

M = M + N(1) * X(1) 
N=N + M(l) 

NEXT | 

M = M/N 

PRINT “LA MEDIA ES:” 
PRINT M 


En estas instrucciones se calculan la suma de frecuencias N y la media aritmética 
M, y se imprime esta última. 


140 
150 
160 
170 
180 
190 


mnm< 
> 


0 
OR|=1TOK 
V + NI) * (X() — M) Y 2 
T! 
V/N 


SQR(V) 


n<z< 
m 
x ll 


De la línea 140 a la línea 180 se calcula la varianza V y en la línea 190 se calcula 
la desviación típica S. 


200 
210 
220 
230 
240 


PRINT “LA VARIANZA ES:” 

PRINT V 

PRINT “LA DESVIACION TÍPICA ES:” 
PRINT S 

END 


Estas instrucciones imprimen los valores de la varianza y la desviación típica. 
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Ejecución del programa 


Si se ejecuta el programa para la siguiente distribución: 


se obtienen estos resultados: 


run 

NUMERO DE VALORES DISTINTOS ? 6 
VALOR, FRECUENCIA? 1,4 

VALOR, FRECUENCIA? 2,8 

VALOR, FRECUENCIA? 3,3 

VALOR, FRECUENCIA? 4,2 

VALOR, FRECUENCIA? 5,7 

VALOR, FRECUENCIA? 6,9 

LA MEDIA ES: 


3.81818181 

LA VARIANZA ES: 
3.42148760 

LA DESVIACION TIPICA ES: 
1.84972636 

OK 

al 


2. Distribución binomial 
Presentación del problema 


Si se realiza un experimento aleatorio n veces y los posibles resultados del experi- 
mento son dos, éxito y fracaso de probabilidades respectivas, cuyas probabilidades 
py q = 1- p, la distribución que se obtiene se llama binomial, de parámetros n 


y p: 
Bín, p) 
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Por ejemplo, al lanzar una moneda al aire 20 veces (n = 20), los posibles resultados 
del experimento son dos, cara y cruz (éxito y fracaso), cuyas probabilidades son 
p=05yq= 1-05 = 0,5, respectivamente. Luego este experimento res- 
ponde a una distribución binomial de parámetros 20 y 0,5: 


B(20;0,5) 


La variable binomial x de parámetros n y p es la variable que toma los valores 0, 
1, 2, ... n (número de éxitos) 


x=0,1,2,..,n 
con probabilidad: 
e px q pps 
¡de x! (n — x)! E 
La media es: M=n-:p 


y la varianza es: V=n-:p-(1-p) 


En el ejemplo anterior la probabilidad de obtener, por ejemplo, 5 caras es: 


20! 
310191 


la media es: M=>20:05=10 
la varianza: V=20:0,5-05:.=:5 
y la desviación típica: S= W5 


Ps = 0505 


El programa 


10 INPUT “PARAMETROS : N, P”; N, P 

20 INPUT “UN VALOR(1) O TODA LA DISTRIBUCION(2)”; A 
30 PRINT “VALOR X”, “PROBABILIDAD PX” 

40 IFA = 1 THEN INPUT “VALOR”; X : GOSUB 1000 

50 IFA = 2 THEN FOR X = 0TO N : GOSUB 1000 : NEXT X 
60 END 


Este es el bloque principal del programa, en el que se introducen los valores de 
los parámetros N y P, y se elige entre calcular la probabilidad de un solo valor 
(A = 1) o toda la distribución (A = 2). 


Se llama también a la subrutina 1000 para que calcule la probabilidad. 
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1000 Z =N: GOSUB 2000 : C = Y 
1010 Z=X: GOSUB 2000 : D = Y 

1020 Z=N-X: GOSUB 2000 : E = Y 

1030 PX=CU(D*E)*P ? Xx*x(1-P) 1 (N-X) 
1040 PRINT X, PX 

1050 RETURN 


En esta subrutina se calcula la probabilidad PX del valor X. Se hacen tres llamadas 
a la subrutina 2000 para que calcule los factoriales de N, X y N — X. 


2000 Y=1 
2010 FOR |= 
2020 Y=Yx* 
2030 NEXT | 

2040 RETURN 


Esta subrutina calcula el factorial de Z, devolviendo como resultado Y. 


Ejecuciones del programa 


e Si se dan como datos N = 20, P = 0.5 y X = 5 se obtiene este resultado: 


run 

PARAMETROS : N, P? 20, 0.5 

UN VALOR (1) O TODA LA DISTRIBUCION (2) 
21 

VALOR X PROBABILIDAD PX 
VALOR?5 

5 0.01478577 

OK 

ES 


Si se desea obtener todas las probabilidades de la distribución B (10, 0.5) basta 
con seguir los pasos que muestra la pantalla: 
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3. Distribuciones bidimensioñales 
Presentación del problema 


Una distribución bidimensional contempla el caso de una población de n ele- 

mentos y en cada individuo de la misma considera dos características X e Y. 

Ejemplo: Una población de 10 personas de las que se consideran la edad (X) y el 

peso (Y). , 

Sean: Xy, X2, ..., X;, ..., Xp el conjunto de valores que puede tomar X, e Y,, Y,, 
y Yjz «+, Yg el conjunto de valores que puede tomar Y. 


La distribución de frecuencias absolutas de los n elementos de la población se ex- 
presa en una tabla como ésta: 


La tabla de la distribución de frecuencias relativas f; se obtiene dividiendo cada 
frecuencia absoluta n; entre n. 


A una distribución bidimensional se le pueden asociar dos distribuciones que se lla- 
man marginales. 


— Distribución marginal respecto de X, que asocia a cada X; el número g;. 
— Distribución marginal respecto de Y, que asocia a cada Y; el número h; 
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La media y la varianza de cada una de estas distribuciones marginales son: 


p p 
MX = 2 g¡X; VX = y gi(x; = MX)? 
| =1 ii 
q p 
MY= 2 hy; W= 2  h¡(Y- MY) 


1 pe 


Para medir la variación conjunta de las dos variables X e Y se define la covarianza 
de la distribución bidimensional: 
p q 
(V= 2 2 fix - My (Y; - MY) 
¡=1l j=1 - 

La relación entre las variables X e Y viene dada por las ecuaciones de las rectas 
de regresión de Y sobre X y de X sobre Y cuyas gráficas se ajustan a las nubes 
formadas por los n puntos (X;, Y;) e (Y;, X;), respectivamente. 


Recta de regresión de Y sobre X: Coeficiente de regresión: 
- Vx E 
Y — MY = VX (X — MX) BY = VX 
Recta de regresión de X sobre Y: Coeficiente de regresión: 
_Ccv Ñ EM 
X — MX = vwY (Y —- MY) BX = Y 


Una vez ajustada una recta de regresión a la nube de puntos se define el coefi- 
ciente de correlación R entre los variables X e Y. 


Ñ CV 
iS VVX - VVY 


que sirve para determinar la dependencia entre las dos variables. 


El programa 


10 INPUT “NUMERO DE VALORES DE X”; P 
20 DIM X(P) 

30 FOR|=1TOP 

40 PRINT “VALOR DE X”; | 

50 INPUT X(1) 

60 NEXT | 

70 CLS 


108 


En este bloque de instrucciones se introducen los distintos valores de X. 


80 INPUT “NUMERO DE VALORES DE Y”; Q 
90 DIM Y(Q) 

100 FORJ=1TOQ 

110 PRINT “VALOR DE Y”; J 

120 INPUT Y(J) 

130 NEXT) 

140 CLS 


En este bloque de instrucciones se introducen los distintos valores de Y. 


150 DIMN(P,Q):N=0 

160 FOR|=1TOP 

170 FORJ=1TOQ 

180 PRINT “FRECUENCIA DE X”; |; *,Y”;J 
190 INPUT N(l, J) 

200 N=N+N(l, )) 

210 NEXT J 

220 NEXT | 

230 CLS 


En este conjunto de instrucciones se introducen los valores de las frecuencias ab- 
solutas n; y se calcula su suma total n. 


240 DIM F(P, Q) 


250 FOR!|=1TOP 
260 FORJ=1TOQ 
270 F(l, J) = NC, J/N 
280 NEXT) 
290 NEXT | 


En estas instrucciones se calculan las frecuencias relativas f;. 


300 DIM G(P) : DIM H(Q) 
310 FORI|=1TOP 

320 FORJ=1T0OQ 

330 G(1) = G(1) + FI, J) 
340 H() = H() + Fl, J) 
350 NEXT J 

360 NEXT | 
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En estas instrucciones se calculan las frecuencias g; y h; de las distribuciones mar- 
ginales X e Y, respectivamente. 


3710 MX=0:MY=0:VX=0:VY=0:CV=0 
380 FOR!|=1TOP 

390 MX = MX + G(l) * X(l) 

400 NEXT 


Las líneas 380-400 calculan la media MX. 


410 FORJ=1T0Q 
420 MY = MY + H()) * Y()) 
430 NEXT J 


Las líneas 410-430 calculan la media MY. 


440 FOR|=1TOP 
450 VX = VX + Gl) * (X() — MX) 1 2 
460 NEXT | 


Las líneas 440-460 calculan el valor de la varianza VX. 


470 FORJ=1TOQ 
480 VY = VY + H()) * (YO) - MY) 1 2 
490 NEXT J 


Las líneas 470-490 calculan el valor de la varianza VY. 


500 FORI=1TOP 

510 FORJ=1TOQ 

520 CV = CV + FI, J) * (X() — MX) * (Y()) — MY) 
530 NEXT J 

540 NEXT | 


Las líneas 500-540 calculan el valor de la covarianza CV. 


550 BY = CVWVX 
560 BX = CV/VY 
570 R= CV(SQR(VX) * SQR(VY)) 
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La línea 550 calcula el coeficiente de regresión BY de la recta de Y sobre X. La lí- 
nea 560 calcula el coeficiente de regresión BX de la recta de X sobre Y. La línea 
570 calcula el coeficiente de correlación R. 


580 
590 
600 
610 
620 
630 
640 
650 
660 


DEF FNR(Z) = INT(Z * 100 + 0.5)/100 
MX = FNR(MX) 

MY = FNR(MY) 

VX = FNR(VX) 

VY = FNR(VY) 

CV = FNR(CV) 

BY = FNR(BY) 

BX = FNR(BX) 

R = FNR(R) 


La línea 580 define la función FNR(Z) que redondea a dos decimales el valor de Z. 


Las líneas 590-660 utilizan esta función para redondear a dos decimales todos los 
valores obtenidos. 


670 
680 
690 
700 
710 
720 
730 
740 


PRINT “VARIABLE MARGINAL X” 
PRINT “MEDIA ="; MX 

PRINT “VARIANZA ="; VX 
PRINT “VARIABLE MARGINAL Y” 
PRINT “MEDIA ="; MY 

PRINT “VARIANZA ="; VY 
PRINT “COVARIANZA” 

PRINT CV 


Estas instrucciones imprimen las medias y varianzas de las distribuciones margi- 
nales, y la covarianza. 


750 
760 
770 
780 
790 
800 
810 


PRINT “RECTA DE REGRESION DE Y SOBRE X” 
PRINT rd de —MY; A BY; "e y MX; e 
PRINT “RECTA DE REGRESION DE X SOBRE Y” 


PRINT X "+ =MG "=*0 BXx: “e (Y, MY)" 
PRINT “COEFICIENTE DE CORRELACION” 
PRINTR 

END 


Este conjunto de instrucciones imprimen las ecuaciones de las dos rectas de re- 
gresión y el valor del coeficiente de correlación. 
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Ejecución del programa 


Sea la tabla de distribución de frecuencias: 


Al ejecutar el programa y al ir Introduciendo estos datos se obtienen sucesiva- 
mente las siguientes pantallas: 


run 

NUMERO DE VALORES DE X ? 5 
VALOR DE X1 
2 

VALOR DE X2 
22 

VALOR DE X3 
23 

VALOR DE X4 
24 

VALOR DE X5 
5 


NUMERO DE VALORES DE Y?4 
VALOR DE Y1 

22 

VALOR DE Y2 

24 

VALOR DE Y3 

26 

VALOR DE Y4 

28 


FRECUENCIA DE X1, Y1 
FRECUENCIA DE X1, Y2 
FRECUENCIA DE X1, Y3 
FRECUENCIA DE X1, YA 
E ECUENOA DE X2, Y1 
FRECUENCIA DE X2, Y2 
FRECUENCIA DE X2, Y3 
FRECUENCIA DE X2, Y4 
-RÉCUENCIA DE X3, Y1 
FRECUENCIA DE X3, Y2 
FRECUENCIA DE X3, Y3 
FRECUENCIA DE X3, Y4 
FRECUENCIA DE X4, Y1 
FRECUENCIA DE X4, Y2 
FRECUENCIA DE X4, Y3 


FRECUENCIA DE XA, Y4 
FRECUENCIA DE X5, Y1 
FRECUENCIA DE X5, Y2 
FRECUENCIA DE X5, Y3 
FRECUENCIA DE X5, Y4 
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VARIABLE MARGINAL X 

MEDIA = 3.25 

VARIANZA = 2.03 

VARIABLE MARGINAL Y 

MEDIA = 5.42 

VARIANZA = 5.01 . 

COVARIANZA 

0.08 

RECTA DE REGRESION DE Y SOBRE X 


Y — 5.42 = 0,04 x (X — 3.25) 

RECTA DE REGRESION DE X SOBRE Y 
X — 3.25 = 0.02 + (Y — 5.42) 
COEFICIENTE DE CORRELACION + 
0.03 

OK 

ñ 


Al ser el coeficiente de correlación casi cero (R = 0,03) las dos variables apenas si 
están relacionadas. 


Como los coeficientes de regresión son prácticamente cero, el ajuste es malo. 


Observación 


Dado que en el programa la tabla N(l, J) sólo se utiliza para calcular la tabla F(l, J), 
se puede ahorrar bastante memoria sustituyendo ambas por la segunda. Para ello 
basta con modificar las líneas siguientes del programa: 


150 DIMF(P,Q):N=0 
190 INPUT F(I, J) 

200 N=N +FK(l, )) 

240 (se suprime esta línea) 
270 F(l, J) = F(l, J/N 


De esta forma, en las líneas 150-230, F(I, J) representa las frecuencias absolutas. 


En la línea 270 se realiza el cambio a frecuencias relativas, y a partir de esta línea 
F(I, J) representa las frecuencias relativas. 
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8. Matemática avanzada 


1. Área de un polígono 
Presentación del problema 


El área de un polígono, en función de sus vértices (x;, y;), viene dada por la fór- 
mula: 


n 
2 (x - A 
=2 


$= > 
donde y,+: = Y; 
(Xn  Yn) 
Y 
dl 
> (x1,y1) 
Y 
(Xs, ys) 5 

2 (x2,y2) 

(X4, Ya) (x3,y3) 


Por ejemplo, el área de este cuadrilátero 
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se obtiene de la siguiente forma: 


5= (1 — %2) * (Ya — y) + (X1 — %3) > (ya — y2) + (1 — %4) > (Ys — ya) 
2 


Q-3):-4(1-3+(Q+1):(4+ )+(Q+4+2):(G4+2) 
2 


[are _ 0 


El programa 


10 INPUT “NUMERO DE VERTICES”; N 

20 DIM X(N) : DIM Y(N + 1) 

30 FORI|I=1TON 

40 PRINTT “COORDENADAS DEL VERTICE”; | 
50 INPUT X(1), Y(l) 

60 NEXT | 

70 Y(N + 1) = Y(1) 


En este grupo de instrucciones se introduce el número de vértices y sus valores. 


80 S=0 
90 FORI=2TON 
100 S=S+ (X(1) — X(1)) * (Y(I + 1) —- Y(Y — 1) 
110 NEXT | 
120 S = ABS(S)/2 
130 PRINT “EL AREA DEL POLIGONO ES:” 
140 PRINTS 
150 END 


El bloque de instrucciones 80-120 calcula el área, y las instrucciones 130 y 140 
imprimen el resultado. 
Ejecución del programa 


Si se ejecuta el programa para los datos N = 4, (2,3), (3,1), (-1,-2) y (22,4) se 
obtiene el siguiente resultado en la pantalla: 
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NUMERO DE VERTICES? 4 

COORDENADAS DEL VERTICE 1 
COORDENADAS DEL VERTICE 2 
COORDENADAS DEL VÉRTICE 3 
COORDENADAS DEL VERTICE 4 


?-2,4 
EL AREA DEL POLIGONO ES: 


2. Operaciones con vectores 
Presentación del problema 


Vamos a estudiar un programa que permita, conociendo las componentes de dos 
vectores, resolver las siguientes operaciones entre dichos vectores: 


a) suma 
b) resta 
c) producto escalar 
a) producto vectorial 


teniendo en cuenta que si las componentes del vector A son (A1, A2, A3) y las 
del vector B (B1, B2, B3), los resultados de estas operaciones se obtienen a partir 
de las expresiones: 


ES 


= (A1 + B1, A2 + B2, A3 + B3) 
= (A1 — B1, A2 — B2, A3 - B3) 
A1 * B1 + A2 * B2 + A3 * B3 (producto escalar) 


= (A2 + B3 — A3 * B2, A3 * B1- Al x B3, Al * B2— A2 + B1) (pro- 
ducto vectorial) 


DD |ID>iDI 
| 

oo] 00/00] ul 
| 


Xx 
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El programa 


1 REM OPERACIONES CON VECTORES 
10 PRINT “COORDENADAS VECTOR A” 
20 INPUT A1, A2, A3 
30 PRINT “COORDENADAS DEL VECTOR B”. 
40 INPUT B1, B2, B3 


50 PRINT 
60 PRINT “A + B = (”; Al + B1; pd A2 + B2; es A3 + B3;”)" 
70 PRINT “E _ B E Ds Al E B1; Pi A2 Ss B2: ¡Pu A3 _ B3;”)" 


80 PRINT“A-B=";A1* B1 + A2* B2 + A3 * B3 
90 PRINT“A*B=("¡A2*B3-A3%B2;",";A3xB1*—A1*B3; 
"/"¡A1*B2-A2xB1;")" 


En estas instrucciones, después de introducir las coordenadas de los vectores (ins- 
trucciones 10 a 40), se obtienen e imprimen los resultados de las distintas opera- 
ciones (instrucciones 60 a 90). 


100 PRINT “MAS DATOS (1 = SI, O = NO)” 
110 INPUT X 

120 IFX= 1 THEN GOTO 10 

130 END 


En esta última parte se pregunta si se introducirán nuevos datos y, en caso afir- 
mativo, se inicia el proceso de nuevo. En caso contrario, se termina el programa. 


Ejecución del programa 


Si se introducen como datos los vectores A(2, 3, 1) y B(1, 0, 2) se obtiene como 
resultado lo que muestra la pantalla: 
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3. Módulos, cosenos directores y ángulo de dos vectores 


Presentación del problema 


Dados dos vectores A y B, si (A1, A2, A3) son las coordenadas del primer vector 


A y (B1, B2, B3), las coordenadas del segundo vector B, son módulos respectivos 
M1 y M2, se obtienen aplicando estas fórmulas: 


M1 =V A1? + A2? + A3? 
M2 = V B1? + B2? + B3? 


Los cosenos del ángulo que forman con los ejes se obtienen calculando estos co- 
cientes: 


A AZ A 
M1 M1 

Br 82 83 
M2 M2 M2 


Y el coseno del ángulo que forman se determina así: 


HZ, _ A1xB1+4+A2 + B2 + A3 x B3 
cos (A, B) - Mi * M2 
Los vectores son perpendiculares si este coseno es cero, y son paralelos si su valor 
es16-1. 


MN 


>> 
Observación. Para calcular el ángulo Z = (A, B) se utiliza la fórmula: 


ME Ep 
Z = arctg > Siendo C = cosí(A, B) 
Esto es necesario tenerlo en cuenta porque algunos microordenadores no dispo- 
nen de la función arco coseno. 


El programa 


a) El módulo de cada uno de ellos. 

b) Los cosenos de los ángulos que forman con los ejes X, Y, Z (cosenos direc- 
tores). 

Cc) El ángulo que forman los dos vectores. 
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El programa averigua, además, si los vectores son paralelos o perpendiculares. 


PRINT “COORDENADAS DEL VECTOR A” 
INPUT A1, A2, A3 
PRINT “COORDENADAS DEL VECTOR B” 
INPUT B1, B2, B3 


A z —> > 
En estas líneas se han introducido las coordenadas de los vectores A y B. 


50 
60 
70 


M1 =SQR(A1 1 2+A2 1 2+A3 1 2) 
M2 = SQR(B1 1 2 +8B2 1 2+B3 1 2) 
PRINT “MODULOS” : PRINT MÍ : PRINT M2 


En las líneas 50 y 60 se calculan los módulos de los dos vectores, y en la línea 70 
se imprimen. 


IF M1 = O THEN GOTO 110 

PRINT “COSENOS DIRECTORES DEL VECTOR A” 
PRINT A1/M1 : PRINT A2/M1 : PRINT A3/M1 

IF M2 = O THEN GOTO 140 

PRINT “COSENOS DIRECTORES DEL VECTOR B» 
PRINT B1/M2 : PRINT B2/M2 : PRINT B3/M2 


En el bloque de instrucciones 80-130 se calculan los cosenos directores de los dos 
ss» 
vectores, A y B. 


140 


150 
160 
170 
180 
190 
200 
210 
220 
230 


IF M1 = 0 OR M2 = O THEN PRINT “NO HAY ANGULO” : GOTO 
240 

C =(A1 x B1 + A2 * B2 + A3 * B3) / (M1 * M2) 

C = INT (C * 1000 + 0.5)/1000 

IFC = O THEN Z = 90 : GOTO 200 

Z = ATN (SQR (1 — C f 2) /C) 

Z=Z+* 180/3.1416 

PRINT “ANGULO ENTRE LOS DOS VECTORES” 

PRINT Z; “GRADOS” 

IFC = O THEN PRINT “LOS VECTORES SON PERPENDICULARES” 
IFC = 1 OR C = -1 THEN PRINT “LOS VECTORES SON PARA- 
LELOS” 
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La línea 150 calcula el coseno C del ángulo que forman y la línea 160 lo redon- 
dea a 3 decimales. 


La línea 170 comprueba si el valor de C es cero, para evitar en la línea 180 se 
produzca un error al calcular el ángulo Z por medio de la función arcotangente. 


La línea 190 pasa el ángulo Z de radianes a grados sexagesimales. 


Las líneas 200-230 imprimen el valor del ángulo e informes en el caso de que los 
vectores sean perpendiculares o paralelos. 


240 PRINT “MAS VECTORES (1 = SI, 0 = NO)” 


250 INPUT | 
260 IFI= 1 THEN GOTO 10 
270 END 


En esta última parte se pregunta si hay más vectores, para reiniciar o no el pro- 
ceso de cálculo. 


Ejecuciones del programa 


s o > 
e Si se introducen como datos los vectores A(1, 1, 1) y B(8, 8, 8), paralelos, se 
obtiene el siguiente resultado: 


Al pulsar O termina la ejecución. 
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e Si se introducen como datos los vectores A(S, 0, 0) y B(0, 7, 8) perpendicular, 
se obtiene el siguiente resultado: 


Al introducir O termina la ejecución. 
> > 


e Si se introducen como datos los vectores A(2, 3, 4) y B(6, 4, 2) se obtiene el 
siguiente resultado: 


Al introducir O termina la ejecución. 


4. Suma de matrices 
Presentación del problema 


Una matriz es una tabla de números que están distribuidos en filas y columnas. 
Por ejemplo, la siguiente distribución de números: 


1.2 columna 2.* columna  3.* columna 


za a) 
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es una matriz de 2 filas y 3 columnas. Si a la matriz anterior se le llama A, el ele- 
mento que ocupa la primera fila y la segunda columna, que es 3, suele indicarse 
as2, y en trabajos de programación, A(1, 2). El primer número indica la fila y el 
segundo la columna. Luego los elementos de la matriz A son: 


A(1, 1) =2 A (1,2) =3 A (1, 3) =4 
A(2, 1) =1 A (2, 2) =0 A (2, 3) = 2 


La matriz A se dice que es de dimensión 2 X 3, pues tiene 2 filas y 3 columnas. 
Para sumar dos matrices A y B, ambas tienen que ser de la misma dimensión 
(igual número de filas e igual número de columnas). La matriz suma se obtiene 
sumando los términos que ocupan la misma posición. Ejemplo: 


(o 2)+ (42 5)= (323) 


El programa 


En el programa, primeramente se debe indicar la dimensión de las mismas, M 
y N. A continuación, hay que definir las matrices A(M, N); B(M, N); S(M, N); 
siendo A-B las matrices sumandos (datos) y S la matriz suma. 


1 REM SUMA DE MATRICES 
10 INPUT “FILAS”, M 
20 INPUT “COLUMNAS”, N 
30 DIM A(M, N) 
40 DIM B(M, N) 
50 DIM S(M, N) 


En las instrucciones INPUT se introducen. las dimensiones (M = número de filas y 
N = números de columnas). En las instrucciones 30, 40 y 50 se dimensionan las 
matrices A, B y S. 


60 PRINT “MATRIZ A” 
70 FORI=1TOM 

80 PRINT “FILA”; | 

30 FORJ=1TON 

100 PRINT “COLUMNA”; J 
110 INPUT A (l, J) 

120 NEXT J 

130 NEXT | 

140 CLS 
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Mediante el bucle 70-130 se hace variar la fila y mediante el 90-120, la columna. 
La instrucción 110 va pidiendo los distintos valores de la matriz A. La 140 borra la 
pantalla. 


150 PRINT “MATRIZ B” 
160 FORI=1TOM 

170 PRINT “FILA”; | 

180 FORJ=1TON 

190 PRINT “COLUMNA”; J 
200 INPUT B (l, J) 

210 NEXT J 

220 NEXT | 

230 CLS 


Este bloque repite el mismo proceso que el anterior, sólo que ahora se introducen 
los valores de la matriz B. 


240 FORI=1TOM 

250 FORJ=1TON 

260 S (Il, J) = A (1, J) + B (Il, J) 
270 NEXT J 

280 NEXT | 


Por medio de los bucles 240-280 (que hacen variar las filas) y 250-270 (que hacen 
variar las columnas) se suman los elementos que ocupan la misma posición en A y 
B (instrucción 260), almacenándose el resultado en SÍ, J). 


290 PRINT “LA MATRIZ SUMA ES” 
300 FORI=1TOM 

310 PRINT “FILA”; | 

320 FORJ=1TON 

330 PRINT S(I, J) 

340 NEXT) 

350 NEXT | 

360 END 


En este grupo de instrucciones se imprime la matriz suma S por filas, dándose por 
terminado el programa. 
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Ejecución del programa 


Si se ejecuta el programa para M = 2 y N = 3, siendo las matrices 


A=(h 0 2) (425) 


se obtiene este resultado: 


En consecuencia, la matriz suma es: 


5. Producto de matrices 
Presentación del problema 


Para poder multiplicar dos matrices A y B es necesario que el número de co- 
lumnas de A coincida con el número de filas de B. 


Si la dimensión de A es M X N, la dimensión de B tiene que ser N x P; es decir, 
el número de columnas de la primera matriz tiene que coincidir con el número de 
filas de la segunda. 


La matriz producto R = A x B tiene por elementos: 


R(1, 1) = suma de los productos de los términos de la 1.* fila de A por los de la 
1.2 columna de B. 


R(1, 2) = suma de los productos de los términos de la 1.* fila de A por los de la 
2.* columna de B. 


R(1, P) = suma de los productos de los términos de la 1.? fila de A por los de la 
P-ésima columna de B. 
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R(2, 1) = suma de los productos de los términos de la 2.* fila de A por los de la 
1.* columna de B. 


R(2, 2) = suma de los productos de los términos de la 2.* fila de A por los de la 
2.? columna de B. 


Ejemplo: 
Sean 


6 3 
13 
Aja ( )y e-(1 0) 
405 de 
E 
ear (2 39( o) - 
28 


1*r6+2*1+3*2 1*3+2*0+3*5 14 18 
TX 4xr6+0:*14+5*x2 4x3+0*0+5x5 ) — 34 37 


Al ser la dimensión de A, 2 X 3, y la de B, 3 x 2, la dimensión de P es 2 x 2 (el 
primer 2 corresponde al número de filas de A y el segundo 2 al número de co- 
lumnas de B). 


El programa 


El programa que efectúe el producto de dos matrices tiene que leer sus dimen- 
siones, dimensionar las matrices A, B y R e introducir los elementos de las ma- 
trices A y B. 


1 REM PRODUCTO DE MATRICES 
10 INPUT “FILAS DE LA PRIMERA MATRIZ”; M 
20 INPUT “COLUMNAS DE LA PRIMERA MATRIZ O FILAS DE LA SE- 
GUNDA MATRIZ”; N 
30 INPUT "COLUMNAS DE LA SEGUNDA MATRIZ”; P 
40 DIM A(M, N) 
50 DIM BÍN, P) 
60 DIM R(M, P) 
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En las instrucciones 10, 20 y 30 se introducen las dimensiones de las dos ma- 
trices, teniendo en cuenta que el número de columnas de la primera es igual al 
número de filas de la segunda. Las instrucciones 40, 50 y 60 dimensionan las ma- 


trices A, B y R. 
70 PRINT “MATRIZ A” 
80 FORI=1TOM 
90 PRINT “FILA”; | 
100 FORJ=1TON 
110 PRINT “COLUMNA”; J 
120 INPUT A(l, J) 
130 NEXT J 
140 NEXT | 
150 CLS 


Este bloque sirve para introducir los elementos de la matriz A, por filas. Se utilizan 
dos bucles: el de variable | para variar las filas, y el interno, de variable J, para va- 
riar las columnas (la fila | fija). Por último, la instrucción 150 borra la pantalla. 


160 
170 
180 
190 
200 
210 
220 
230 
240 


PRINT “MATRIZ B” 
FOR | = 1TO N 
PRINT “FILA”; | 

FOR J = 1 TOP 
PRINT “COLUMNA”; J 
INPUT B(l, J) 

NEXT J 

NEXT | 

CLS 


Este bloque es análogo al anterior; solamente se diferencia en que introduce los 
elementos de la matriz B. 


250 
260 
270 
280 
290 
300 
310 


R(l, J) = R(l, J) 
NEXT K 
NEXT J 

NEXT | 


A(l, K) * B(K, J) 
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El bucle 250-310, de variable |, va fijando la fila de R y la de A. 
El bucle 260-300, de variable J, fija la columna de R y la de B. 


El bucle 270-290, de variable K, calcula el valor del elemento R(l, J) por medio de 
la instrucción 280. 


Por tanto, en este bloque se calculan los distintos elementos de la matriz pro- 
ducto R. 


320 PRINT “LA MATRIZ PRODUCTO ES” 
330 FOR|=1TOM 

340 PRINT “FILA”; | 

350 FORJ=1TOP 

360 PRINT R(I, J) 

370 NEXT J 

380 NEXT | 

390 END 


Este grupo de instrucciones imprime los valores de las distintas filas de la matriz 
producto R. Con la instrucción 390 termina el programa. 


Ejecución del programa 


Si se ejecuta el programa dando como datos: 


M=2 N=3 P=2 


6 3 
(133) o0( o) 
2 5 


el resultado que se obtiene en la pantalla es éste: 
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es decir: 
_ (14 18 
ES En a) 


6. Resolución de sistemas de ecuaciones lineales 


Presentación del programa 


Vamos a analizar el caso de los sistemas lineales en los que el número de ecua- 
ciones coincide con el número de incógnitas. Trataremos primero de determinar si 
tiene solución y, en caso afirmativo, hallarla. 


Aunque el programa es válido para un sistema cualquiera de N ecuaciones con N 
incógnitas, estudiaremos el procedimiento refiriéndonos a un sistema de 3 ecua- 
ciones con 3 incógnitas: 


A(1, 1) * X(1) + A(1, 2) * X(2) + A(1, 3) * XB) = A(1, 4) 
A(2, 1) * X(1) + A(2, 2) * X(2) + A(2, 3) * X(3) = A(2, 4) (1) 
A(3, 1) * X(1) + AG, 2) * X(2) + AB, 3) * XB) = AG, 4) 


donde A(1, 1), ..., A(3, 3) son los coeficientes de las incógnitas X(1), X(2) y X(3), y 
A(1, 4), A(2, 4) y A(3, 4), los términos independientes. 


Vamos a transformar este sistema en otro equivalente, reemplazando las ecua- 
ciones 2.* y 3.*, y dejando fija la 1.* Para ello realizamos las siguientes opera- 
ciones: 


— Multiplicamos la 2.* ecuación por A(1, 41) y la 1.2 por A(2, 1), y restamos 
ambas miembro a miembro. Obtenemos así una nueva ecuación, que reem- 
plaza a la 2.*, con el primer término nulo. 


— Multiplicamos la 3.? ecuación por A(1, 1) y la 1.2 por A(3, 1), y restamos 
ambas miembro a miembro. Obtenemos de este modo una nueva ecuación, 
que reemplaza a la 3.?, con el primer término nulo. 

Obtenemos así el siguiente sistema, equivalente al anterior: 

A(1, 1) * X(1) + A(1, 2) * XQ) + A(l, 3) * = A(1, 4) 
A(2, 2) * X(2) + A(2, 3) * XG) = A(2, a (2) 

) * X(2) + A ) * = AB, 4) 
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Transformamos este sistema en otro equivalente, dejando fijas la 1.? ecuación y la 
2.?, y reemplazando la 3.* Para ello multiplicamos la 3.* ecuación por A'(2, 2) y la 
2.* por A'(3, 2), y restamos ambas miembro a miembro. Obtenemos así una 
nueva ecuación, que reemplaza a la 3.*, con el primer y segundo términos nulos. 
Queda, entonces, el siguiente sistema equivalente: 


A(1,1)* X()+A(1, 2)* X(2) +A (1,3) * X(3)=A (1, 4) 
A(2, 2) * X(2) + A' (2, 3) * X(3) = A' (2, 4) (3) 
A”"(3, 3) * X(3) = A“, 4) 


En este sistema, a partir de la 3.? ecuación obtenemos el valor de X(3), dividiendo 
A"(3, 4) entre A”(3, 3). Reemplazando el valor de X(3) en la 2.* ecuación obte- 
nemos X(2) y, finalmente, reemplazando los valores de X(3) y X(2) en la 1.* ecua- 
ción, hallamos el valor de X(1). 


Este procedimiento es posible siempre que algún coeficiente de X(1) sea distinto 
de cero; en el ejemplo, estos coeficientes son A(1, 1), A(2, 1) y AZ, 1). 


Si alguno de estos coeficientes es distinto de cero; el procedimiento conduce al 
sistema equivalente (2). En éste se repite el proceso con las ecuaciones del sis- 
tema, con exclusión de la primera. Así se llega al sistema equivalente (3). 


Ahora bien, para que el sistema tenga solución es necesario que el coeficiente 
A(N, N) (en el ejemplo, A”(3, 3)) sea distinto de cero. 


El programa 


10 PRINT “NUMERO DE ECUACIONES O INCOGNITAS” 
20 INPUT N 

30 DIMA(M, N + 1) 

40 FORI=1TON 

50 PRINT “ECUACION”; | 

60 FORJ=1TON 

70 PRINT “COEFICIENTE”; J 

80 INPUT A(l, J) 

90 NEXT J 
100 PRINT “TERMINO INDEPENDIENTE” 
110 INPUT A(I, N + 1) 
120 NEXT | 
130 CLS 


Esta parte del programa introduce en la memoria los coeficientes de las incógnitas 
y los términos independientes de las N ecuaciones. 
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La instrucción 20 introduce el número de ecuaciones del sistema. 


140 FORK =1TON-— 1 

150 IF A(K, K) <> O THEN GOTO 250 
160 FOR|=K+1TON 

170 IF A(I, K) <> O THEN GOTO 200 
180 NEXT | 

190 GOTO 430 


Mediante estas instrucciones se analiza si algún coeficiente A(K, K) es distinto de 
cero y se busca uno que lo sea. 


200 FORM=KTON + 1 
210 B=A(l, M) 

220 A(l, M) = A(K, M) 
230 A(K, M) = B 

240 NEXT M 


Con las instrucciones 200 a 240 se permutan las ecuaciones (tener en cuenta que 
todavía no se ha cerrado el bucle de variable K). 


Para el cambio de las ecuaciones se utiliza la variable B en la que se almacenan 
los resultados. 


250 FOR!|=K+1 
260 FORJ=K+1 
270 All, J) = All, J) 
280 NEXT J 
290 NEXT | 
300 NEXT K 


TON. 
TON +1 
+ A(K, K) — A(K, J) * A(l, K) 


En esta parte del programa se obtienen sucesivos sistemas equivalentes hasta lle- 
gar a uno en el que figura una ecuación gon una sola incógnita. Una vez conse- 
guido este sistema se cierra el bucle de variable K. 


310 IF A(N, N) = O THEN GOTO 430 

320 PRINT “SISTEMA COMPATIBLE. LA SOLUCION ES:” 
330 FOR K = N TO 2 STEP -1 

340 B=A(K,N + 1)/A(K, K) 

350 PRINT *X”; K; “="; B 

360 FOR|=1TOK- 1 

370 A(I, N + 1) = A(I, N + 1) — A(l, K) * B 
380 NEXT | 

390 NEXT K 

400 X1 =A(1,N + 1/A(1, 1) 

410 PRINT "X 1 ="; X1 
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Mediante estas instrucciones se han ido obteniendo y sustituyendo el valor de las 
incógnitas hasta llegar al valor de X1. 


420 GOTO 440 
430 PRINT “EL SISTEMA NO ES COMPATIBLE DETERMINADO” 
440 END 


La instrucción 420 transfiere el control a la 440, después de ser impreso el valor 
de la incógnita X?. 


En el caso de que los coeficientes A(K, K) sean todos nulos, la instrucción 190 
transfiere el control a la 430, y si el coeficiente A(N, N) es nulo, la instrucción 310 
también transfiere el control a la instrucción 430. 


Con la instrucción 440 se termina el programa. 


Ejecución del programa 


Si al ejecutar el programa se introduce la dimensión N = 4 y los coeficientes y 
términos independientes del sistema 


2x1 — 12 + 1 — 2Xa =-5 
2X1 + 2X2 — 3X3 + 1Xy = -1 
=1X1 + 1X, — 1X3 + 0x4 = -1 
AX, — 3X2 + 2X3 — 3X4 = 8 


se obtiene en la pantalla el siguiente resultado: 


SISTEMA COMPATIBLE. LA SOLUCION ES: 


x4=3 
x3=2 
x2=1 
x1=0 
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7. Resolución de una ecuación por el método del punto medio 
Presentación del problema 
Supongamos que f(x) es una función continua en el intervalo [a, b], y los signos 


de f(a) y de f(b) son distintos. Con estas condiciones, se puede demostrar que la 
ecuación 


TO) =0 


tiene solución en el intervalo (a, b] 


Solución 


El procedimiento de cálculo numérico que hay que seguir para obtener la solución 
con una precisión E determinada es el siguiente: 


— se divide el intervalo por 2, obteniéndose así el punto medio c del intervalo 
la, b] (de ahí el nombre del método); 

— se verifica si la longitud del intervalo la, b] es menor que E; es decir, se com- 
prueba si la — b| < E. Si es menor, la solución es c; 


— se verifica si f(c) > 0; en caso afirmativo, se hace b = c y en caso negativo se 
hace a = c. De esta forma se obtiene un nuevo intervalo [a, b] más pequeño 
que sigue cumpliendo las condiciones iniciales; 


— se repite el proceso hasta lograr que c (solución) < E. 


El programa 


Para elaborar el programa hay que tener en cuenta que: 


FNF(x) es la función que hay que definir en cada caso. 

AyB son los extremos del intervalo. Estos valores deben cumplir la condición de 
que FNF(A) y FNF(B) tengan distinto signo. 

E es la precisión con que se desea obtener la solución. 

E es el punto medio de los sucesivos intervalos que se van obteniendo. 
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Aunque el programa que transcribimos empieza en la instrucción 20, en cada 
caso concreto hay que añadir una instrucción 10 que defina la función correspon- 
diente a la ecuación cuya solución se desee hallar. 


20 INPUT “EXTREMOS DEL INTERVALO”; A, B 

30 INPUT “PRECISION”; E 

40 C=(A + BJ/2 : 

50 IF ABS(B — A) < E THEN PRINT “SOLUCION ="; C : END 
60 IF FNF(C) > O THEN B = CELSEA = C 

70 GOTO 40 


En las instrucciones 20 y 30 se introducen los valores de A, B y E. 
La línea 40 calcula el punto medio C del segmento AB. 


La línea 50 pregunta si el valor absoluto de la diferencia de A y B es menor que la 
precisión E; en caso afirmativo imprime la solución C y termina el proceso. 


La línea 60 pregunta si la función en C es positiva; en caso afirmativo hace B = C 
y en caso contrario A = C. 


La línea 70 envía el control a la 40 para que se repita el proceso. 


Ejecución del programa 


Si se desea hallar la solución de la ecuación 
+x=x-1=0 
previamente hay que definir la función 
(0) =x + x=x-=1 
Luego la instrucción 10, deberá ser ésta: 

10 DEF FNF(x) =x3+x2 2-x- 1 
Elegimos como extremo del intervalo a = 0 y b = 3. 
Dado que 

fla) = f(0) = -1 

f(b) = f(3) = 32 
son de distinto signo, la solución se encuentra en el intervalo [0, 3]. 
Al ejecutar el programa y dar los datos: 


m0 > 
IT 


0 
3 
0.001 
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se obtiene el resultado que aparece en la pantalla. 


EXTREMOS DEL INTERVALO? 0.3 
PRECISION? 0.001 
SOLUCION = 1.0001220703125 


Observación 


A este programa se le puede añadir una línea 25 para comprobar si la función 
tiene el mismo signo en los extremos A y B, en cuyo caso, imprime, un informe 
que expresa que estos extremos no valen y después envía el control a la instruc- 
ción 20, para que se introduzcan nuevos valores de A y B. 


25 IF FNF(A) * FNF(B) > = O THEN PRINT “NO VALEN” : GOTO 20 


8. Resolución de la ecuación cúbica 
Presentación del problema 


La ecuación general de tercer grado es de la forma 
Ax? + Bx?+Cx+D=0 


Para resolverla es necesario realizar una serie de cambios de variable. Los pasos 
que hay que dar son los siguientes: 


1. Hacer este cambio de variable: 


Desarrollando y simplificando queda: 


B? ) (E BC ) 
5 — —— o — — = 
AY +( a FO + A o 


Y dividiendo por A se obtiene: 


ec 22.2 BCO, 0D 
3 = .=— — = — — = 
, +( A) e 3Az * o) y 


¡ A Y AA 
Haciendo o Cial 
3 Fr 
Q= 28 BCO, D 


2783 3? A 


La ecuación anterior queda así: 
Y? +3PY+Q=0 


Hacer este otro cambio de variable: 
Y=U+V 
Desarrollando 
U? + 3U%V + 3UV? + V? + 3PU+V+0Q=0 
y simplificando queda: 
U* + Y? + (3UV + 3P)(U+VW+0Q=0 
Entre U, V y P imponemos esta relación: 
3UV + 3P =0 
de la que se deduce 


UV =-— Po su equivalente V = o 


UÚ 
Sustituyendo V por = se obtiene la expresión 
p3 
y? E 1 + Q = 0 


o su equivalente 
Uf + QUÉ - P?=0 
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que es una ecuación cuadrática en U? ((U”) + Q(U?) — P? = 0) cuya solución 
es: 
-Q+ VQ + 4p> 
Z 


Si se hubiera despejado V en función de U se habrían obtenido las mismas 


soluciones. En consecuencia se puede elegir: 
-Q- V_Q? + ap? 


ya - -Q+ VQ? + Ap? 
2 2 


Luego se obtienen tres soluciones para U y otras tres para V, que combinadas 
dan nueve soluciones, pero de ellas sólo se deben tomar las que cumplan la 
condición: 


U3 = 


Ves 


U-V=-P 


Con esta restricción, las tres raíces de la ecuación cúbica vienen dadas por la 
fórmula de Cardano: 


SES Y QUEAP” y al NO RAP” 
2 2 


en la que el producto de las raíces cúbicas es igual a —P. 


Para calcular las soluciones se tiene que tener en cuenta: el valor del discrimi- 
nante de la ecuación U* + QU? — P?=0 


Tar | 


e SiT<O0 existen tres raíces reales: 


¡ERA TT IA E VQ) -” - 
2 A 2 


= V Lí(cos Z + ¡sen Z) 
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NT E A A: 
Va = VI (cos (3 + n) -isen (3 + n)) 


_ YT CS. 14 E Z 4 
Ya = VI (cos (3 + 1") ¡ sen (2 + 1) 


Y, = WT cos (2 + 21) -7 -P cos (2 + an) 


yn 
E 
=2 VW L c0s (2 + +1) =2 VW cos (2 + ES) 


e SiT = 0 existen dos raíces reales, una de ellas doble. 
En este caso U = V, y las raíces son: 
Y, = + 2U 
Y, = —U (raíz doble) 


e SiT>0 existen una raíz real y dos complejas conjugadas. 
Y, = U + V (raíz real) 


_  U4+V ESTE" 
Y, = 2 A “a (U V)i 
Ed A 
e = > > (U — Vi 


Por último, hay que deshacer el cambio X = Y — => para obtener las solu- 


ciones de la ecuación Ax? + Bx? + CX+D=0. 
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El programa 


10 INPUTA, B, C, D 
20 P=-B 1 2(9*A 1 2) + 0(3*A) 

30 Q=(2*B 1 3)/(27*A 1 3)-(8*0)/(3*A f 2) + D/A 
40 T=Q12+4+*P13 

50 IF T< O THEN GOSUB 1000 : END 

60 U=(- Q + SQR(T)/2 

70 IF U<O THEN U = —((-U) 1 (1/3) 

80 IF U>O0 THEN U =U ) (1/3) 

90 IF T = 0 THEN GOSUB 2000 : END 

100 V =(-Q - SQR(T))/2 

110 IF V<OTHEN V = - ((-V) 1 (1/3) 

120 IF V>OTHEN V = V % (1/3) 

130 IF T> O THEN GOSUB 3000 : END 

140 END 


Éste es el bloque principal del programa, que llama a tres subrutinas, 1000, 2000 
y 3000, para resolver la ecuación según los valores del discriminante T. 

— La línea 10 pide los coeficientes de la ecuación. 

— Las líneas 20 y 30 calculan los valores de P y Q. 

— La línea 40 calcula el valor del discriminante T. 


— Si T<O0 la línea 50 envía el control a la subrutina 1000 para que encuentre 
las soluciones. 


— Las líneas 60, 70 y 80 calculan el valor U. Distinguiendo entre los casos de 
U<0U>0yU=0. 


— Si T =0 la línea 90 envía el control a la subrutina 2000 para que encuentre 
las soluciones. 


— Las líneas 100, 110 y 120 calculan el valor de V. 


— SIT > 0 entonces la línea 130 envía el control a la subrutina 3000 para que 
encuentre las soluciones. 


1000 PRINT “TRES RAICES REALES” 

1010 IF -0.01 < Q AND Q < 0.01 THEN Z = 1.5708 : GOTO 1030 
1020 ee ATN (SQR (-Q + Q/4 — P 1 3) / ABS 0/2) 

1030 IF (-Q/2) < O THEN Z = 3.1416 — Z 

1040 PRINT 2 + SQR(-P) + cos (2/3) — B/(3 + A) 

1050 PRINT 2 * SQR(-P) + cos (2/3 + 6.2832/3) — B/(3 * A) 

1060 PRINT 2 * SQR(-P) + cos (2/3 + 2 * 6.2832/3) — B/(3 * A) 
1070 RETURN 


Esta subrutina se ejecuta cuando T < 0 y da como resultado tres raíces reales. 
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Las líneas 1010, 1020 y 1030 calculan el valor del argumento Z, averiguando si 
está en el primer cuadrante o en el segundo, o si es un ángulo recto (línea 1010), 
en cuyo caso no se calcula en la línea 1020 para que no dé error (ya que no 
existe la tangente de un ángulo recto). 


Las líneas 1040, 1050 y 1060 imprimen los resultados. 


2000 PRINT “DOS RAICES REALES (UNA DOBLE)” 
2010 PRINT 2 * U — B/(3 * A) 

2020 PRINT — U — B/(3 * A); “RAIZ DOBLE” 
2030 RETURN 


Esta subrutina se ejecuta cuando T = 0, y da como resultado dos raíces reales, 
una de ellas doble. 


3000 PRINT “UNA RAIZ REAL Y DOS COMPLEJAS” 
3010 PRINT U + V — B/(3 * A) 

3020 R=-— (U + V)/2 — B/(3 + A) 

3030 | = ABS (SQR(3)/2 * (U — V)) 

3040 PRINTR; “+"; 1; “1” 

3050 PRINTR; “—”; 1; %i” 

3060 RETURN 


Esta subrutina se ejecuta cuando T > 0 y da como resultado una raíz real y dos 
complejas conjugadas. 
La línea 3010 imprime la raíz real. 


La línea 3020 calcula la parte real de las raíces complejas, y la 3030 la parte ima- 
ginaria en valor absoluto. 


Las líneas 3040 y 3050 imprimen las raíces complejas. 


Ejecuciones del programa 


e La ecuación X? — 6x? + 11x — 6 = 0 tiene como soluciones x = 1; x = 2; 
x = 3 como muestra la ejecución del programa 


run 

21,6, 11,6 

TRES RAICES REALES 
2.99999929 


1.00000031 
1.99999733 
OK 

5 
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e La ecuación x? + 3x? + 3x + 1 = 0 tiene como soluciones x = 1: x =-1; 
x = —1. En la pantalla aparece una raíz doble: 


run 

21,3, 3,1 

DOS RAICES REALES (UNA DOBLE) 
-1 

—1 RAIZ DOBLE 


OK 
Y 


e La ecuación x? — 1 = 0 tiene una solución real x = 1 y dos complejas 


E] 


X= 7; a ¡. Al ejecutar el programa se obtiene: 


run 

21,0, 0, -1 
UNA RAIZ REAL Y DOS COMPLEJAS 
1 : 

0.5 + 0.86602540 ¡ 


0.5 — 0.86602540 ¡ 
OK 
5 
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9. El calendario 


1. Fechas de la Semana Santa de cualquier año 
Presentación del problema 


e La Semana Santa viene determinada por la fecha de la Pascua de Resurrección. 
La fecha de esta conmemoración religiosa es celebrada el primer domingo des- 
pués del plenilunio (luna llena) que ceincide o sigue inmediatamente al equi- 
noccio de primavera (21 de marzo). 


El algoritmo que se expone a continuación permite determinar la fecha de la 
Pascua (día y mes), y es válido desde el año 1583, siguiente a la implantación 
del calendario gregoriano (15-10-1582). 


El calendario gregoriano acumula el error de un día al cabo de 3314 años, 
luego el equinoccio de primavera (21 de marzo) del año 1582 + 3314 = 4896 
tendrá lugar el día 20 de marzo (un día antes). 


Por eso el algoritmo no es válido para calcular la Pascua de los años poste- 
riores al 4896. 


Las fechas extremas de la Pascua son el 22 de marzo (cae en esta fecha si el 
22 es domingo) y el 25 de abril (cae en esta fecha si el 20 de marzo es plenilu- 
nio). 


Luego la Pascua puede caer en cualquier día que se encuentre dentro de 
ambas fechas mencionadas (variación: 35 días). 


e Algoritmo: 


El año x 
X 
B 
B+8 
B-F+1 
MA+B-D-G+15 
C 
32 + 2E+21-H-) 
A+ 11H + 22K 
H+K-7L+ 114 


EZ T|-—-]|aonuoow]| 
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Como resultado se tiene: M, que da el mes y N + 1, que da el día de la Pas- 
cua el año x. 


La Semana Santa será del lunes anterior, al domingo N + 1. 


La Semana Santa del año x se extiende, entonces, desde el lunes anterior al 
domingo N + 1, hasta dicho domingo. 


El programa 
10 INPUT “AÑO”; X 
20 IFX< 1583 OR X > = 1582 + 3314 THEN GOTO 10 
30 Z=X:Y=19: GOSUB 1000: A=R 
40 Z=X:Y= 1000 : GOSUB 1000:B=Q:C=R 
50 Z=B:Y=4:GOSUB 1000:D=Q:E=R 
60 Z=B+8:Y=25: GOSUB 100:F=0Q 
710 2Z=B-F+1:Y=3: GOSUB 1000:G=Q 
80 Z=19*A+B-D-G+15:Y=30: GOSUB 1000: H =R 
9 Z=C:Y=4:GOSUB 1000 :1=Q:J]=R 
100 Z=32+2*E+2*!|-H-J:Y=7: GOSUB 1000 : K =R 
110 Z=A+11+H+22*K: Y = 451 : GOSUB 1000 :L =Q 
120 Z=H+K-7*L+ 114: Y = 31: GOSUB 1000: M=0Q:N=R 


En este conjunto de instrucciones se introduce el valor de año X, y realizan las di- 
visiones que indica el algoritmo, llamando a la subrutina 1000 que calcula el co- 
ciente Q y el resto R de la división entera Z/Y. 


130 
140 
150 
160 
170 
180 


190 


DIM M$(12) 

FOR | = 1 T0 12 

READ M5(1) 

DATA “ENERO”, “FEBRERO”, “MARZO”, “ABRIL” 

DATA "MAYO”, “JUNIO”, “JULIO”, “AGOSTO” 

DATA SEPTIEMBRE”, “OCTUBRE”, “NOVIEMBRE”, “DICIEM- 
BRE” 

NEXT | 


El bucle 140-190 asigna a la lista M$(1) los meses del año. 


200 PRINT “LA PASCUA DEL AÑO”; X; “ES EL DIA”; N + 1; “DEL 
MES DE”; M$(M) 
210 END 
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La línea 200 imprime el día y el mes en que cae el domingo de Pascua del año X. 


1000 Q = INT(Z/Y) 
1010 R=Z-YxQ 
1020 RETURN 


La subrutina 1000 calcula el cociente Q y el resto R de la división entera de Z/Y. 


Ejecuciones del programa 


e Para averiguar las fechas de la Semana Santa del año 1986 basta con seguir 
los pasos que indica la pantalla 


run 

AÑO ? 1986 

LA PASCUA DEL AÑO 1986 ES EL DIA 30 
DEL MES DE MARZO 

OK 


Luego la Semana Santa del año 1986 tendrá lugar del lunes 24 de marzo al 
domingo 30 de marzo. 


e La Semana Santa del año 2000 caerá en las fechas que se obtienen a conti- 
nuación 


run 

AÑO 22000 

LA PASCUA DEL AÑO 2000 ES EL DIA 23 
DEL MES DE ABRIL 

OK 


Luego, la Semana Santa del año 2000 tendrá lugar del 17 al 23 de abril. 
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2. 


Día de la semana y día juliano, DJ 


Presentación del problema 


José Scaliger propuso, en 1582, una escala de tiempo continua, de 7980 años ju- 
lianos de 365,25 días, que comienza el 1 de enero del año 4713 a. de C. 


A esta escala de tiempo la llamó período juliano en honor de su padre Julio. 


Para obtener el día juliano DJ de un año A de nuestro calendario basta con apli- 
car las siguientes reglas: 


1: 
Z: 


4, 


Se calcula el producto (4712 + A) - 365,25. 


Si el resultado es un número entero se le resta 1, y si es decimal, se toma su 
parte entera. 


Si el año es siguiente al de la reforma gregoriana (4-10-1582) se tienen que 
realizar las siguientes operaciones a continuación del anterior: 


1583 <A < 1700 : se resta 10 (*) 
1701 <A < 1800 : se resta 11 
1801 <A < 1900 : se resta 12 
1901 <A < 2100 : se resta 13 
2101 <A < 2200 : se resta 14, etc. 


Se suma el resultado, el número de días transcurridos desde el comienzo del 
año A. 


Para determinar el día de la semana basta con dividir por 7 el día juliano DJ, y se- 
gún sea el resto, se harán estas asignaciones: 


0 => LUNES 

1 —> MARTES 

2 => MIERCOLES 
3 > JUEVES 

4 > VIERNES 

5 > SABADO 

6 => DOMINGO 


(*) Se debe restar 10 cuando A > 1583, y también si A = 1582, siempre que se verifique 
M>106M=10yD>= 15. 
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El programa 


10 INPUT “DIA, MES, AÑO”; D, M, A 

20 Z= (4712 + A) * 365.25 

30 IF Z = INT(Z) THEN Z = Z — 1 ELSE Z = INT(Z) 

40 IF(A >= 1583 OR (A = 1582 AND (M > 10 OR (M = 10 AND 
D >= 15)))) ANDA <= 1700 THENZ=Z- 10 

50 IF 1701 < = A AND A < = 1800 THEN Z = Z — 11 

60 IF 1801 < = A AND A < = 1900 THEN Z = Z — 12 

70 IF 1901 < = AAND A < = 2100 THEN Z = Z — 13 

80 IF 2101 < = AAND A < = 2200 THEN Z = Z — 14 


En este conjunto de instrucciones se introduce el día D, el mes M y el año A, y se 
calcula Z, que es el número de días transcurridos según el período juliano, hasta 
que comienza el año A. 


90 FOR!I=1TOM 

100 READ B 

110 DATA O, 31, 59, 90, 120, 151, 181, 212, 243, 273, 304, 334 
120 NEXT | 


Este bucle asigna a B es número de días de los meses anteriores a M. 


130 IF (A/4 <> INT (A/4)) THEN GOTO 180 

140 IFA < 1582 THEN GOTO 160 

150 IF (A/100 = INT(A/100)) AND (A/400 <> INT(A/400)) THEN 
GOTO 180 

160 IFM< = 2 THEN GOTO 180 

170 B=B+1 


Las líneas 130-170 determinan si el año es bisiesto: 


— La línea 130 comprueba si A no es divisible por 4, en cuyo caso el año no es 
bisiesto. 


— La línea 140 comprueba si el año pertenece al período juliano. 


— La línea 150 elimina los bisiestos seculares que no sean divisibles por 400, se- 
gún la reforma gregoriana. 


— La línea 160 comprueba si el mes es enero o febrero, en cuyo caso no influye 
que el año sea bisiesto. 
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— La línea 170 añade un día al computo, por ser el año bisiesto. 


180 Z=Z+B+D 
190 Q = INT(Z/7) 
200 R=Z-7*x0Q 


La línea 180 calcula el día juliano Z y, la 200, el resto de su división entre 7. 


210 RESTORE 240 

220 FORI=0TOR 

230 READ R$ 

240 DATA “LUNES”, “MARTES”, “MIERCOLES”, “JUEVES”, 
“VIERNES”, “SABADO”, “DOMINGO” 

250 NEXT | 


El bucle 220-250 lee el día de la semana R$ correspondiente al resto R. 


260 PRINT “EL DIA”; D; “DEL MES”; M; “DEL AÑO”; A; “ES”; R$ 
270 PRINT “EL DIA JULIANO DJ CORRESPONDIENTE ES”; Z 
280 END 


La línea 260 imprime el día de la semana y la 270 el día juliano DJ. 


Ejecuciones del programa 


e ¿En qué día de la semana nació Carolina si la fecha de su nacimiento fue el 13 
de abril de 1983? 


run 

DIA, MES, AÑO? 13, 4, 1983. 

EL DIA 13 DEL MES 4 DEL ANO 1983 

ES MIERCOLES 

EL DIA JULIANO DJ CORRESPONDIENTE ES: 


2445438 
OK 
ES 
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e ¿En qué día de la semana comenzará el año 2000? 


run 

DIA, MES, AÑO? 1, 1, 2000 

EL DIA 1 DEL MES 1 DEL AÑO 2000 ES 

SABADO 

EL DIA JULIANO DJ CORRESPONDIENTE ES: 
2451545 


OK 
$ 


3. Determinación del número de días que hay entre dos fechas 
Presentación del problema 


El procedimiento que permite determinar el mismo número de días que hay entre 
dos fechas se reduce a calcular los días julianos correspondientes a las dos fechas 
y a restarlos. Sugerimos, en consecuencia, repasar los conceptos expuestos en el 
programa anterior. 


El programa 


10 INPUT “DIA, MES, AÑO”; D1, M1, A1 

20 INPUT “DIA, MES, AÑO”; D2, M2, A2 

30 D=D1:M=M1:A=A1: GOSUB 100 : Z1 =Z 

40 D=D2:M=M2:A=A2: GOSUB 100 : Z2 =Z 

50 PRINT “EL NUMERO DE DIAS ENTRE EL”; D1; ””; M1; “”; A1; 
“VENDA MZA, AZ AES> 

60 PRINT ABS(Z2 — Z1) 

70 PRINT “EL DIA JULIANO CORRESPONDIENTE AL”; D1; ”,”; M1; 
CAT": PES:"; 21 

80 PRINT “EL DIA JULIANO CORRESPONDIENTE AL”; D2; ”.”; M2; 
cid A2; "ES:": 22 

90 END 
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Las líneas 10-90 constituyen el programa principal que permite introducir las dos 
fechas, llama a la subrutina 100 para que calcule los dos números julianos corres- 
pondientes, y por último, imprime el número de días entre las dos fechas y sus 
días julianos. 


100 Z = (4712 + A) * 365.25 

110 IF Z = INT(Z) THEN Z = Z — 1 ELSE Z = INT(Z) 

120 IF (A >= 1583 OR (A = 1582 AND (M < 10 OR (M = 10 AND D 
< = 15)))) AND A < = 1700 THEN Z =Z — 10 

130 IF 1701 <= A AND A < = 1800 THEN Z = Z — 11 

140 IF 1801 < = A AND A < = 1900 THEN Z = Z — 12 

150 IF 1901 < = A AND A < = 2100 THEN Z = Z — 13 

160 IF 2101 < = A AND A < = 2200 THEN Z = Z — 14 


En la línea 100 comienza la subrutina que calcula el número juliano Z. En el blo- 
que de instrucciones 100-160 se calcula el número de días julianos que hay hasta 
el principio del año A. 


170 RESTORE 200 

180 FOR!|=1TO M 

190 READ B 

200 DATA O, 31, 59, 90, 120, 151, 181, 212, 243, 273, 304, 334 
210 NEXT | 


El bucle 180-210 calcula el número B de días transcurridos del año A antes del 
mes M. 


220 IF (A/4 <> INT (A/4)) THEN GOTO 270 

230 IFA < 1582 THEN GOTO 250 

240 IF (A/100 = INT (A4/100)) AND (A/400<>INT (A/400)) THEN 
GOTO 270 

250 IFM< = 2 THEN GOTO 270 

260 B=B+1 

2710 2Z=Z2+B+D 

280 RETURN 


La línea 220 comprueba si el año no es bisiesto. La línea 230 comprueba si el año 
es anterior al de la reforma gregoriana del calendario (15-10-1582). La línea 240 
elimina de los bisiestos los años seculares no divisibles por 400. La línea 250 com- 
prueba si el mes M es enero o febrero. 
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La línea 260 suma 1 al número de días transcurridos del año A si es bisiesto y si 
el mes M es distinto de enero o febrero. 


La línea 270 calcula el día juliano Z. 


Ejecución del programa 


Si se desea saber el número de días que hay entre el 15 de febrero de 1963 y el 
29 de julio de 1985 basta con ejecutar el programa introduciendo los datos, 
como muestra la pantalla: 


run 

DIA, MES, AÑO? 15, 2, 1963 

DIA, MES, ANO? 29, 7, 1985 

EL NUMERO DE DIAS ENTRE EL 15, 2, 1963 y 
EL 29, 7, 1985 ES: 

8200 


EL DIA JULIANO CORRESPONDIENTE AL 
15, 2, 1863 ES: 2438076 

EL DIA JULIANO CORRESPONDIENTE AL 
29, 7, 1985 ES: 2446276 

OK 

a 


Observación 


Una posible aplicación del programa es calcular el número de días vividos hasta la 
fecha. Para ello basta con introducir la fecha de nacimiento y la fecha de hoy. 
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10. Programas para estudiar 
conceptos de fisica y química 


1. Cálculo de la aceleración de un sistema de dos masas 


Presentación del problema 


Sea el sistema de dos masas representado en esta figura. 


ma 


rm8 


mag 


Se trata de hallar la aceleración a que adquiere el sistema. Para ello aplicamos la 
segunda ley de Newton a las masas m1 y m2: 


m,9g-T=m,a 
T-rm29g=m,a 


siendo m, y m, los valores de las masas, T la tensión de la cuerda, r el coeficiente 
de rozamiento de la masa m, con la superficie horizontal, y g = 9,8 m/s”. 


Resolviendo el sistema, se obtiene el valor de la aceleración a. 
_ (m,-—rm)):9,8 
m; + m> 


Si el rozamiento no es elevado, la aceleración a puede ser positiva, con lo que el 
sistema caería de forma acelerada. 


Si, en cambio, el rozamiento es lo suficientemente elevado pueden presentarse 
estas situaciones: 


a) Quea<0a=0 y v (velocidad inicial) sea nula. En este caso, el sistema no 


se movería. 
b) Que a =0 y v = 0. En este caso, el sistema se movería con velocidad cons- 
tante. 
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El programa 


El programa resuelve el problema planteado, teniendo en cuenta los distintos 
casos que se pueden presentar. 


1 REM ACELERACION 
10 PRINT “MASA 1” 
20 INPUT M1 
30 PRINT “MASA 2” 
40 INPUT M2 
50 PRINT “COEFICIENTE DE ROZAMIENTO” 
60 INPUT R 
70 PRINT “VELOCIDAD INICIAL” 
80 INPUT V 


En estas instrucciones se introducen los datos para resolver el problema. 


90 A= (M1 — R * M2) * 9.8 / (M1 + M2) 
En la instrucción 90 se calcula el valor de la aceleración. 


100 IFA = O THEN GOTO 140 

110 IFA <O THEN GOTO 160 

120 PRINT “LA ACELERACION ES”; A; “M/(S 1 2)” 

130 GOTO 170 

140 IFV = O THEN GOTO 160 

150 PRINT “LA VELOCIDAD ES SIEMPRE IGUAL A”; V; “M/S” 
160 PRINT “EL SISTEMA NO SE MUEVE” 

170 END 


En las instrucciones 100-110 se averigua si la aceleración es positiva; en caso afir- 
mativo se imprime su valor (instrucción 120) y, en caso contrario, se tienen en 
cuenta las distintas alternativas. 


Ejecución del programa 


Si después de teclear RUN y | RETURN | se introducen estos valores 
M1 = 5 kg 
M2 = 1 kg 
R = 0.85 
V = 0.2 m/s 
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se obtiene el resultado que muestra la pantalla: 


2. Cálculo del período de revolución de un satélite 


Presentación del problema 


El período T de revolución de un satélite en órbita alrededor de la Tierra viene ex- 


presado por la fórmula 
a / (R +A)? 
T=2n R2 


g 


siendo R el radio de la Tierra, que se supone igual a 6370000 metros, y la acele- 
ración de la gravedad (9,8 m/s”, aproximadamente) y A la altura sobre la superfi- 
cie terrestre. (Para que la fórmula dé resultados correctos A debe ser mayor que 
200 km y menor que 50 000 km.) 


El programa 


El programa calcula el período T e imprime en la pantalla el período de revolución 
en horas y minutos. 


1 REM PERIODO DE REVOLUCION 


10 G=9.8 
20 R= 6370000 
30 Pl = 3.14159 


40 INPUT “ALTURA DEL SATELITE EN KM” 

60 IFA > 200 AND A < 50000 THEN GOTO 90 
70 PRINT “DATO INCORRECTO” 

80 GOTO 40 
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En estas instrucciones se asignan valores a las variables G. R y Pl (no es necesario 
definir el valor de Pl si el microordenador tiene esta función) y se introduce des- 
pués el valor de la altura, comprobándose si es una altura permitida (instruc- 
ción 60). 


9 A=A* 1000 

100 T=2*Pl* SQR ((R + A) $ 3/(G * R + R)) 
110 H = INT (T/3600) 

120 MN = INT ((T — 3600 * H)/60) 

130 PRINT “PERIODO DE REVOLUCION” 

140 PRINT H; “HORAS”; MN; “MINUTOS” 


En esta parte del programa se calcula el-período de revolución del planeta (ins- 
trucción 100). Previamente, en la instrucción 90 se expresa la altura en metros. En 
las instrucciones 110 y 120 se expresa el período de revolución en horas y mi- 
nutos, pues en la instrucción 100 está expresado en segundos. 


150 PRINT “HAY OTRA ALTURA” 


160 INPUT S$ 
170 1FS$ = “SI” THEN 40 
180 END 


Finalmente se pregunta si se va a introducir una nueva altura y, en caso afirma- 
tivo, se empieza de nuevo el cálculo. 


Ejecución del programa 


Si después de teclear RUN y | RETURN | se introduce 2300 para A, se obtiene el 
valor 2 horas 14 minutos para T. 


3. Análisis de las ondas de una cuerda 
Presentación del problema 


Cuando en un medio lineal, como por ejemplo una cuerda, se establece en un 
punto (el foco) una vibración armónica, ésta se propaga por la cuerda y cada 


154 


punto de ella vibra con un movimiento vibratorio armónico pero desfasado con 
respecto al foco. 


La ecuación que describe el movimiento de todos los puntos de la cuerda es: 


Ñ DM eo DL 
Y =Asen (2 t L 3 


siendo Y la elongación de cada punto, A la amplitud, P el período (tiempo trans- 
currido para realizar una oscilación), L la longitud de onda, t el tiempo transcu- 
rrido y R la distancia al foco. 


En la ecuación se aprecia una doble periodicidad: una con respecto al tiempo, ya 
que cuando t = P (período) el movimiento en un punto se repite, y otra con res- 
pecto a la distancia, ya que los puntos situados a 1 L, 2 L ... vibran en fase. 


A esta distancia L se llama longitud de onda. 


Se llama frecuencia F al número de oscilaciones que experimenta un punto del 
medio por unidad de tiempo, verificándose que F = 1/P. 


Por otro lado, existe una relación entre la longitud de onda y el período, y es la 
siguiente: 


L=vT=VHF 


en la que v es la velocidad de propagación de la perturbación (que depende de 
las características del medio y no del movimiento del foco). 


El programa 


El programa siguiente analiza la doble periodicidad respecto al tiempo y a la dis- 
tancia. 


1 REM ONDAS 
10 INPUT “AMPLITUD”; A 
20 INPUT “FRECUENCIA”; F 
30 INPUT “VELOCIDAD”; V 
40 L=VF:P=1/F 
50 DEF FN F(T, X) = A + SIN (2 * 3.14 * (T/P — XIL)) 


En este conjunto de instrucciones se introducen los datos y se calculan la longitud 
de onda y el período. 
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En la instrucción 50 se define la función que describe a la onda. 


60 T1=0:T2=3x*P:Y1=-A:Y2=A 
70 DT = (T2 — T1) / 245 : DY = (Y2 — Y1) / 180 


La instrucción 60 determina los valores máximo y mínimo del tiempo t y de la 
elongación y. 


La instrucción 70 determina el valor de las escalas que se van a utilizar en la re- 
presentación gráfica. 


80 SCREEN 2 ; 
90 FOR X =0TO 2 + L STEP L/2 
100 LINE (0, 90) — (255, 90) 
110 LINE (15, 0) — (15, 180) 
120 FORT =0TO 3 * P STEP DT 
130 Y =FN F(T, X) 
140 X3 = INT ((T — 0) / DT + 0.5) 
150 Y3 = INT ((Y — Y1) / DY + 0.5) 
160 PSET (15 + X3, 180 — Y3) 
170 NEXT T 
180 CLS 
190 NEXT X 
200 GOTO 200 
210 END 


Este grupo de instrucciones dibuja las gráficas de los movimientos vibratorios de 
puntos situados a distancias X = 0, 1/2, L, L + U2 y 2L del foco. Para ello se uti- 
liza el bucle 90-190. Cada uno de estos movimientos vibratorios se representa 
para el tiempo transcurrido entre O y 3 * P (bucle 120-170). 


Ejecución del programa 

Si después de teclear RUN y | RETURN | se introducen estos valores 
A=0.5m 
F = 800 Hz 
V= 1200 ms” 
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se obtiene este resultado: 


En las gráficas de esta pantalla se comprueba que para cada longitud de onda L 
se repite el mismo movimiento. 


4. Estudio de las interferencias 
Presentación del problema 


Un punto de un medio, por ejemplo de una cuerda, cuando es alcanzado simul- 
táneamente por dos movimientos, vibra con un movimiento que es suma del que 
producirían cada uno de ellos independientemente. 


Consideremos el caso de dos movimientos ondulatorios de la misma amplitud y 
frecuencia. 


La interferencia se produce debido a la diferencia de caminos recorridos. 


La ecuación que representa el movimiento de cada punto de la cuerda puede ob- 
tenerse sumando las funciones de onda de cada uno de los movimientos: 


Y = Y1 + Y2 = A sen (2n(T/P — X1/L)) + A sen (2n(T/P — X2/L)) 


El programa 


Mediante el siguiente programa se demuestra que si X2 — X1 = KL (K = 1, 2, 
3, ...), siendo L la longitud de onda, el efecto de la interferencia corresponde a un 
máximo y si X2 — X1 = (2K — 1) U2 (K = 1, 2, 3, ...) se produce una interferen- 
cia destructiva. 
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REM INTERFERENCIAS 


CLS 


INPUT “FRECUENCIA”; F 

INPUT “VELOCIDAD”; V 

INPUT “DISTANCIA AL FOCO 1”; X1 

INPUT “DISTANCIA AL FOCO 2”; X2 

L=VF:P=1/F 

DEF FN F(T) =:A « SIN (2 * 3.14 * (T/P — X1/L)) + A + SIN 


(2 * 3.14 x (T/P — X2/L)) 


Con estas instrucciones se introducen los datos, se calculan la longitud de onda y 
el período y se define la función que se representa gráficamente. 


80 T1=0:T2=3*P:Y1=-2*A:Y2=2*A 


90 


DT = (T2 — T1) / 240 : DY = (Y2 — Y1)/ 180 


Mediante estas instrucciones se determinan los valores máximo y mínimo del 
tiempo t y de la elongación y. 


100 
110 
120 
130 
140 
150 
160 
170 
180 
200 
210 


SCREEN 2 

LINE (0, 90) — (250, 90) 

LINE (15, 0) — (15, 180) 

FOR T = OTO 3 * P STEP DT 
Y = FN F(T) 

X3 = INT ((T — 0) / DT + 0.5) 
Y3 = INT ((Y — Y1) / DY + 0.5) 
PSET (15 + X3, 180 — Y3) 
NEXT T 

GOTO 200 

END 


Las instrucciones 110 y 120 dibujan los ejes coordenados. Las restantes instruc- 
ciones dibujan la gráfica de la elongación e respecto al tiempo. 


Ejecución del programa 


Después de teclear RUN y | RETURN | introducir, por ejemplo, estos valores: 


A = 


0.5 m 


F = 50 HZ 
V = 340 m/s 
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Se obtiene, entonces, una longitud de onda de 6.8 m, y según el valor de X2 — X1 
aparecerán máximos y mínimos. 


5. Valoración ácido fuerte-base fuerte 
Presentación del problema 


La ecuación (no simplificada) que expresa el equilibrio ácido fuerte-base fuerte 
viene dada por la expresión 


x? + (Cg — Ca)x-Kg=0 


siendo x la concentración molar de iones H20* presentes en la disolución, Cg y Ca 
las concentraciones molares de la base y el ácido en cada etapa de disolución, y 
Kg el producto iónico del agua, que depende de la temperatura. 


Se trata de estudiar la variación del pH en el curso de una valoración de un ácido 
fuerte, por ejemplo, ácido clorhídrico con una base fuerte, como puede ser hi- 
dróxido de sodio, o a la inversa. 


De las dos raíces de la ecuación sólo interesa la positiva. Sin embargo, debido al 
pequeño valor de Kg (107** a 25 %C) y al número de cifras significativas que ma- 
neja un microordenador, podría encontrarse que la raíz cuadrada llegara a tomar 
el valor cero, dando una solución de x negativa; ésta, al calcular el logaritmo 
(pH = -log [H30*)), produciría error, dado que no existen logaritmos decimales 
de números negativos. 


Para solucionar este problema, calculamos la mayor de las dos raíces en valor ab- 
soluto. Si es positiva, es la que buscamos, y si es negativa, calculamos el producto 
de las raíces (x1 + x2 = c/a = K£) para obtener la otra. 
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El programa 


10 PRINT “CARACTERISTICAS DE LA DISOLUCION A VALORAR” 
20 PRINT “SI ES ACIDA ESCRIBIR 1, SI ES BASICA ESCRIBIR 2” 
30 INPUT P 

40 IFP= 1 0ORP = 2 THEN GOTO 60 

50 GOTO 30 


En estas instrucciones se determinan las características de la disolución a valorar. 


po] 


60 PRINT “MOLARIDAD DEL PRODUCTO INICIAL” 


70 INPUT M1 
80 IF M1 <o0.001 OR M1 > 1 THEN 100 
90 GOTO 120 ] 
100 PRINT “DATOS INCORRECTOS” 
110 GOTO 70 , 
120 PRINT “VOLUMEN INICIAL EN MILILITROS” 
130 INPUT V1 


140 IFV1 <0 OR V1 > 100 THEN 160 

150 GOTO 180 

160 PRINT “DATO INCORRECTO” 

170 GOTO 130 

180 PRINT “DISOLUCION UTILIZADA EN LA VALORACION” 
190 PRINT “MOLARIDAD” 

200 INPUT M2 

210 IF M2 < 0.001 OR M2 > 1 THEN 230 

220 GOTO 250 

230 PRINT “DATO INCORRECTO” 

240 GOTO 200 

250 PRINT “VOLUMEN TOTAL VERTIDO EN MILILITROS” 
260 INPUT VT 

270 IFVT <= 1 ORVT > 10 + V1 THEN 290 

280 GOTO 310 

290 PRINT “DATO INCORRECTO” 

300 GOTO 260 

310 PRINT “CUANTOS MILILITROS SE AÑADEN CADA VEZ” 
320 INPUT IV 

330 IF IV < VT/20 OR IV > VT/5 THEN 350 

340 GOTO 370 

350 PRINT “DATO INCORRECTO” 

360 GOTO 320 


Hasta aquí se han introducido los datos necesarios verificando su corrección. 
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En las instrucciones 80, 140, 270, 330 se han establecido las condiciones habi- 
tuales en un laboratorio. 


370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 


PRINT 

PRINT “V”, “PH” 

FOR V = 0 TO VIT STEP IV 
C1 = M1 * V1NV1 + V) 
C2 = M2 * V/(V1 + V) 
IFP = 1 THEN 440 
GOTO 470 

CA = C1 

CB = C2 

GOTO 490 

CA = C2 

CB = C1 

GOSUB 600 

PRINT V, PH 

NEXT V 

STOP 


Con este grupo de instrucciones se calculan las concentraciones de las sustancias, 
teniendo en cuenta que sean ácido o base, y se envía el control del programa 
(cada vez que es necesario) a la subrutina 600 para calcular el pH. 


600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 


KE =1E - 14 

D = SQR ((CB — CA) 1 2 + 4 + KE) 
H1 = (- (CB - CA) + D)/2 
H2 = (- (CB - CA) - D)/2 
IF (CB — CA) < O THEN 670 
H = -KE/H2 

GOTO 680 

H=H1 

PH = - LOG (H) / 2.3 
RETURN 

END 


Mediante esta subrutina se calculan las soluciones de la ecuación, teniendo en 
cuenta las indicaciones que se hicieron previamente. 
Posteriormente se calcula el pH (pH = — log [H307] = — LOG [H] / 2.3). 


En la instrucción 680 LET PH = — LOG (H) / 2.3, la función LOG (X) calcula el 
logaritmo neperiano, de ahí que haya que dividir por la constante 2.3 para obte- 
ner el logaritmo decimal (*). 


(*) Ver libro. 


Ejecución del programa 


Después de teclear RUN y | RETURN | introducir, por ejemplo, estos datos: 
P=1 


M1 = 0.05 
VÍ = 85 
M2 = 0.1 
VT=5 
IV = 0.4 


El resultado que se obtiene es el que muestra la pantalla. 


6. Gráficas de puntos de las temperaturas de fusión, de ebullición 
y de las energías de ¡onización de los 40 primeros elementos 
del sistema periódico 


Presentación del problema 


La temperatura de fusión y de ebullición, así como la energía de ionización de los 
elementos son propiedades periódicas cuya variación depende de su número ató- 
MICO. 


La variación de estas propiedades periódicas se puede apreciar representándolas 
gráficamente. 
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El programa 


El siguiente programa representa gráficamente: 


— el punto de fusión en función del número atómico; 
— el punto de ebullición en función del número atómico, y 
— la energía de ionización en función del número atómico. 


90 
100 
110 
120 
130 


DIM Y (40), M (40), B (40), E (40), H (40), V (40), X (40) 

FOR | = 1 T0 40 

READ M (1), B (1), E (1) 

NEXT | 

PRINT “ESCRIBE:” 

PRINT : PRINT “1 — GRAFICA PUNTO DE FUSION — NUMERO 
ATOMICO” 

PRINT “2 — GRAFICA PUNTO DE EBULLICION — NUMERO ATO- 
MICO” 

PRINT “3 — GRAFICA ENERGIA DE IONIZACION — NUMERO ATO- 
MICO” 

PRINT “SI QUIERES SALIR DE UNA GRAFICA” 

PRINT “PULSA UNA TECLA” 

INPUT A 

IFA < 1 OR A > 3 THEN GOTO 50 

ON A GOTO 1000, 2000, 3000 


Este conjunto de instrucciones lee los datos correspondientes a los distintos ele- 
mentos y presenta las opciones que dispone el programa. 


200 
210 
220 
230 
240 
250 
260 
265 
270 
280 
290 
300 


REM DIBUJA 

CLS 

X1=0:X2=40 

Y2 = Y (1) : Y1 = Y(1) 

FOR | = 2T0 40 

IF Y2 < Y (1) THEN Y2 = Y (1) 
IF. Y1 > Y (1) THEN Y1 = Y (1) 
X(1) =1 

NEXT | 

IFY1 > O THEN Y1 =0 

DX = 40 : DY = Y2 — Y1 

XS = 245/DX : YX = 175/DY 


Estas instrucciones permiten escoger los valores máximo y mínimo de la variable 
Y (1), así como establecer las escalas de acuerdo a las cuales se realizará el dibujo. 
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310 
320 
330 
340 
350 
360 
370 
375 
380 
390 
400 


SCREEN 2 

FOR | = 1T0 40 

H.(I) = X (1) + XS : V (1) = (Y (1) + ABS (Y1)) * YS 
NEXT | 

PSET (H (1), 175 — V (1)) 

FOR | = 2T0 40 


LINE (H (1 — 1), 175 — V (1 — 1)) — (H (1), 175 — V (1) 
NEXT | 

A$ = INKEY$ 

IF AS = ”” THEN GOTO 380 

CLS : SCREEN O : GOTO 50 


Este conjunto de instrucciones representa la gráfica Y (I) — | en la pantalla del mi- 
croordenador, uniendo con tramos rectos los distintos puntos. 


1000 
1010 
1020 
1030 
2000 
2010 
2020 
2030 
3000 
3010 
3020 
3030 


En estos 


FOR | = 1T0 40 
Y (1) = M (1) 
NEXT | 

GOTO 200 

FOR | = 1 TO 40 
Y (1) = B (1) 
NEXT | 

GOTO 200 

FOR | = 1 TO 40 
Y (1) = E (1) 
NEXT | 

GOTO 200 


bloques de instrucciones la variable Y (I) adquiere el valor de M (1), B (1) o 


E (1), según la opción. 


4000 
4001 


REM DATOS 

DATA -259, -253, 13.6 
DATA -270, -269, 24.6 
DATA 180, 1360, 5.4 
DATA 1285, 2470, 9.3 
DATA 2030, 3700, 8.3 


Estas instrucciones DATA contienen sucesivamente los puntos de fusión, ebulli- 
ción y energía de ionización de los primeros 40 elementos. 
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Ejecución del programa 


Si después de teclear RUN y | RETURN | se introduce el valor 3 para A, se obtiene 
la gráfica de la energía de ionización en función del número atómico. 


7. Valores energéticos de los alimentos 
Presentación del problema 
La determinación del valor energético de un alimento se consigue quemando en 


un crisol una determinada cantidad del mismo dentro de una cámara rodeada de 
agua y por la que pasa una corriente de oxígeno. 


TERMÓMETRO 
TEMPERATURA 
=20C VOLUMEN DE AGUA 
= 1000 cc 
= 1000 g =1 kg 
ALIMENTO 


OXÍGENO 


Los gases calientes procedentes de la combustión pasan a través de un serpentín 
rodeado de agua, lo que hace que la temperatura del agua se eleve. 
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El valor energético se calcula mediante la fórmula 


AQ=mAt:0,42 


en la que A Q es el contenido energético del alimento, m su masa en gramos, 
A t la variación de temperatura (siendo la temperatura inicial 20C) y 4.2 el 
equivalente mecánico del calor (1 caloría = 4.18 julios). 


El programa 


CLS 

PRINT “EL PESO PUEDE VARIAR ENTRE 0.5 y 3.0 GRAMOS" 
PRINT “PULSA UNA TECLA PARA CONTINUAR” 

AS = INKEYS 

IF AS = "” THEN GOTO 40 . 

INPUT “TECLEA EL CODIGO NUMERICO PARA EL ALIMENTO”; F 
RESTORE 

FOR | = 1TO F 

READ F$ 

READ J 

NEXT | 

PRINT FS 


La instrucción 100 pide el código numérico del alimento. 


200 


210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 


INPUT "ESCRIBE LA MASA DEL ALIMENTO EN GRAMOS 
(0.5 — 3.0)"; M 

T1 =J * M/(4.2 * 1000) + 20 

T = INT (T1 * 10 + 0.5)/10 

PRINT *T ="; T 

CLS 

PRINT "PULSA UNA TECLA PARA CONTINUAR” 

A$ = INKEYS 

IF A$ = "" THEN GOTO 260 

PRINT “LISTA DE DATOS” 

PRINT : PRINT “ALIMENTO: ”; F$ 

PRINT “MASA DE LA COMIDA ="; M; “G” 

PRINT “VOLUMEN DE AGUA = 1000 cc” 

PRINT “MASA DE AGUA = 1000 G” 

PRINT “TEMPERATURA INICIAL = 20 C” 

PRINT “TEMPERATURA FINAL = ”; T 

PRINT “CAMBIO EN LA TEMPERATURA ="; T — 20 


Con estas instrucciones se calcula la temperatura final y se imprimen los resul- 


tados. 


166 


PRINT “PULSA UNA TECLA PARA CONTINUAR” 
A$ = INKEY$ 
IF AS = ”” THEN GOTO 370 


PRINT "1 GRAMO DE AGUA REQUIERE 4.2 JULIOS (J) DE ENER- 
GIA” 

PRINT “PARA ELEVAR SU TEMPERATURA UN GRADO CENTI- 
GRADO” 

PRINT : PRINT “VALOR ENERGÉTICO () =" 

PRINT : PRINT “4.2 * ELEVACION DE TEMPERATURA (C) +* 
MASA DE AGUA (G)”; ”/”; “MASA DE ALIMENTO (G)” 

PRINT : PRINT “EL VALOR ENERGÉTICO (EN JULIOS)”; F$; “="; 
q 

PRINT “EL VALOR ENERGÉTICO (EN KJ) DE ”; F$;” ="; J/1000; 
ao 

PRINT “PULSA UNA TECLA PARA CONTINUAR” 

A$ = INKEYS 

IF A$ = ”” THEN GOTO 470 

CLS 

PRINT “PULSA E PARA EMPEZAR” 

PRINT “PULSA T Y RETURN PARA ACABAR” 

INPUT 15 

IF 15 = “E” THEN GOTO 10 

IF 15 = “T” THEN END 

DATA MANZANAS, 2000 

DATA BACON, 20000 

DATA MANTEQUILLA, 38000 

DATA ZANAHORIAS, 800 

DATA QUESO, 20000 

DATA POLLO, 6000 


La instrucción 1000 y las siguientes almacenan los nombres de los alimentos y el 
valor energético de una masa fija de ellos. 


A cada uno de estos alimentos se le asocia un código numérico que corresponde 
al lugar que ocupan en las instrucciones DATA. 


Así, en el orden arbitrario supuesto a las manzanas les correspondería el código 1, 


al bacon el 2, a la mantequilla el 3, y así sucesivamente. Este código es impor- 


tante para localizar el elemento cuando se lee. 


Ejecución del programa 


Si después de teclear RUN y | RETURN | se introduce el valor 3 para el código F, 


en la pantalla se obtiene el valor energético de la mantequilla. 
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11. Dos programas para utilizar 


en el laboratorio 


1. Pérdida de masa de un sólido 


El programa siguiente indica los pasos que deben seguirse en el laboratorio para 
realizar un experimento, registra los resultados y realiza las operaciones matemá- 
ticas necesarias para obtener conclusiones. 


90 
100 
110 
120 
130 
140 


150 
160 
170 
180 
190 
200 


REM LABORATORIO 

PRINT “PERDIDA DE MASA DE UN SOLIDO CUANDO SE CA- 
LIENTA” 

PRINT : PRINT “PULSA UNA TECLA PARA CONTINUAR” 
A$ = INKEYS$ 

IF A$ = ”” THEN GOTO 30 

CLS 

PRINT “INSTRUCCIONES” 

FOR K = 1 TO 1000 : NEXT K 

PRINT : PRINT “1 — COGE UN CRISOL VACIO” 

FOR K = 1 TO 1000 : NEXT K 

PRINT : PRINT "2 — PESA EL CRISOL VACIO” 

FOR K = 1 TO 1000 : NEXT K 

PRINT : PRINT “3 — PESA EL CRISOL + EL SOLIDO” 

FOR K = 1 TO 1000 : NEXT K 

PRINT : PRINT “4 — CALIENTA HASTA QUE NO HAYA CAM- 
BIO...” 

FOR K = 1 TO 1000: NEXT K 

PRINT : PRINT 5 — DEJA ENFRIAR...” 

FOR K = 1 TO 1000 : NEXT K 

PRINT : PRINT “6 — PESA EL RECIPIENTE + RESIDUO” 
FOR K = 1 TO 1000 : NEXT K 

PRINT : PRINT “7 — REGISTRA TODOS LOS RESULTADOS” 


En este conjunto de instrucciones se dan las indicaciones necesarias para realizar 
el experimento. Después de cada indicación se introduce un tiempo de espera con 
la ayuda de un bucle. 
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205 
210 
220 
230 
240 
300 
310 
320 


330 
340 
350 
360 


370 
400 
410 
420 
430 
440 
450 
460 


470 


PRINT : PRINT : PRINT 

PRINT “PULSA UNA TECLA PARA CONTINUAR” 

A$ = INKEY$ 

IF A$ = ”” THEN GOTO 220 

CLS 

PRINT “MEDIDAS” : PRINT 

INPUT “CUAL ES LA MASA DEL CRISOL VACIO”; A 

IFA < = 0 THEN PRINT “EL CRISOL DEBE PESAR ALGO” : FOR 
K = 1 TO 1000 : NEXT K : GOTO 310 

CLS 

INPUT “CUAL ES LA MASA DEL CRISOL + SOLIDO”; B 

A1=A 

IF B < = A1 THEN PRINT “EL CRISOL + SOLIDO DEBEN SER 
MAS PESADOS : FOR K = 1 TO 1000 : NEXT K : GOTO 310 

CLS 

PRINT “CALIENTA HASTA QUE NO HAYA CAMBIOS...” 

FOR K = 1 TO 1000 : NEXT K 

PRINT “Y DEJA ENFRIAR” 

FOR K = 1 TO 1000 : NEXT K 

CLS 

PRINT “CUAL ES LA MASA DEL” : INPUT “CRISOL + RESIDUO”; C 

IF C > = B THEN PRINT “EL SOLIDO DEBE SER MAS LIGERO 
ANTES QUE DESPUES DE CALENTAR” : FOR K = 1 TO 1000 : 
NEXT K : GOTO 450 - 

IFC < = A1 THEN PRINT “EL CRISOL Y EL RESIDUO DEBEN PE- 
SAR MAS QUE EL CRISOL VACIO” : FOR K = 1 TO 1000 : NEXT 
K : GOTO 450 


Mediante estas instrucciones se introducen los resultados de las medidas y se re- 
chazan los incorrectos. 


500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 


CLS 

PRINT “RESULTADOS” , 

PRINT “1 — MASA DEL CRISOL ="; A; “G” 

PRINT “2 — MASA DEL CRISOL + SOLIDO ="; B; “G” 

$1 =B-— A: S=INT ($1 * 100)/100 

PRINT “3 — MASA DEL SOLIDO ="; B; "—"; A; "="; S; “G” 
FOR K = 1 TO 1000 : NEXT K 

CLS 

FI =C-A:F=INT (F1 * 100)/100 

11 =S-—F:T= INT (T1 * 100)/100 

PRINT “MASA DEL RESIDUO ="; F; “G” 

PRINT “PERDIDA DE MASA ="; T; “G” 

W1 = (T/S) * 100 : W = INT (W1 * 100)/100 

PRINT “% PERDIDA DE MASA”; T; */”; S; * * 100”; “="; W 
FOR K = 1 TO 1000 : NEXT K 
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Este bloque de instrucciones realiza las operaciones necesarias, redondeando a 
dos cifras e imprime en la pantalla los resultados. 


700 CLS 

710 PRINT “RESULTADOS FINALES” : PRINT : PRINT 

720 PRINT “MASA DEL CRISOL ="; A; “G" 

730 PRINT : PRINT “MASA DEL SOLIDO ="; S; “G” 

740 PRINT : PRINT “MASA DEL RESIDUO ="; F; “G” 

750 PRINT : PRINT “PERDIDA DE MASA POR CALENTAMIENTO ="; 
T, "6" 

760 PRINT “PORCENTAJE DE PERDIDA DE MASA ="; W; “%" 

770 END 


Estas últimas instrucciones imprimen los resultados de las medidas tomadas en el 
experimento y calculan el porcentaje de pérdida de masa. 


2. Análisis gráfico de datos experimentales 
Presentación del problema 


Antes de realizar un determinado experimento en el laboratorio, se establece la 
hipótesis de la existencia de una relación entre dos variables X, Y. Estas variables 
pueden ser, por ejemplo, la presión y el volumen de una masa de un gas, la ve- 
locidad y el tiempo en la caída de un cuerpo, la diferencia de potencial y la 
intensidad en un circuito, etc. 


Para encontrar la relación entre las variables X, Y se registra una serie de datos en 
el curso del experimento, correspondientes a ambas variables, y después se dibuja 
la gráfica, Y — X. Esta permite observar si entre ambas variables existe o no rela- 
ción y, en caso afirmativo, determinar la expresión matemática de la misma. Por 
ejemplo, si los puntos de la gráfica pertenecen a una línea recta, la relación entre 
las variables X, Y es una proporcionalidad, cuya expresión matemática es Y = K X + 
+ B, siendo K la constante de proporcionalidad y B el término independiente. 


Cuando la gráfica Y — X no es una línea recta, se tantean otras posibilidades, di- 
bujando por ejemplo las gráficas Y — X?, Y — 1/X, Y — 1/X?, ... Si la gráfica obte- 
nida es una recta, la relación entre X, Y será Y = K X? Y = K (1/K), Y = K 
(1/X?), ..., respectivamente, siendo K la constante de proporcionalidad. 
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El programa 


El programa que se transcribe permite encontrar la relación, si existe, entre las va- 
riables X, Y, analizando distintas gráficas. 


PRINT “ANALISIS DE DATOS” 

PRINT : PRINT 

INPUT “CUANTOS DATOS”; N 

DIM X (N) : DIM Y (N) : DIM Z (N) : DIM H (N) : DIM V (N) 
FOR | = 1TON 

INPUT *X ="; X (1) : INPUT “Y ="; Y (1) 
NEXT | 

CLS 

PRINT EA nd 

PRINT 

FOR | = 1TON 

PRINT X (1), Y (1) 

NEXT | 

PRINT “PULSA UNA TECLA” 

A$ = INKEY$ 

IF AS = ”” THEN GOTO 150 
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Con estas instrucciones se introducen los datos (medidas del experimento), y se 
presentan a continuación en la pantalla en forma de tabla. 


170 
180 
190 
200 
210 
220 
230 
240 
250 


CLS 

PRINT “PUEDES ELEGIR ENTRE LAS GRAFICAS: ” 
PRINT : PRINT “1) Y — X” 

PRINT : PRINT “2) Y — XA 2” 

PRINT : PRINT "3) Y — 1/X” 

PRINT : PRINT “4) Y - 1/X A 2” 

INPUT A 

IFA < 1 ORA > 4 THEN GOTO 170 

ON A GOTO 1300, 1400, 1500, 1600 


Con las instrucciones 190, 200, 210 y 220 se elige la opción que se desee y con 
las instrucciones 1300, 1400, 1500 y 1600 se obtienen los valores de la variable 
independiente que se van a representar (según la opción, la variable indepen- 
diente, Z (X), será X, X*, 1/X Ó 1/X?). 


300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 


REM DIBUJA 

CLS 

22=Z(1):21 =2 (1): Y2 =Y (1): Y1 =Y (1) 
FOR|=2TON 

IF 22 < Z (1) THEN 22 = Z (1) 
IF 21 > Z (1) THEN Z1 = e 
IF Y2 < Y (1) THEN Y = Y (1) 
IFY1 > Y (1) THEN Y1 = Y (1) 
NEXT | 

IF Z1 > O THEN Z1 =0 

IF Y1 > 0 THEN Y1 =0 

DZ = Z2 - 21 : DY = Y2 — Y1 
ZS = 240/DZ : YS = 175/DY 


Estas instrucciones escogen el mayor y menor valor de las variables y con ellos se 


establecen las escalas que se van a utilizar. 
La instrucción 
430 OPEN “GRP:” FOR OUTPUT A$ 4 1 
abre el canal 1, permite imprimir textos en la pantalla de alta resolución. 
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440 
450 
460 
470 
480 
490 
500 
510 


SCREEN 2 

FOR | = 1TON 

H (1) = (Z (1) + ABS (Y1)) * Z5 : V (1) = (Y (1) + ABS (Y1)) * YS 
NEXT | 

PSET (15 + H (1), 175 — V (1)) 

FOR | =2TON 

LINE (15 + H (1 — 1), 175 — V (1 — 1)) — (15 + H (1), 175 — V (1) 
NEXT | 


Con estas instrucciones se dibujan los puntos obtenidos y se unen entre sí con 
tramos rectos. 


520 
530 
535 
537 
540 
550 
560 
570 


PSET (15, 175) 

PRINT 4 1, "PULSA E PARA ELEGIR” 

PSET (15, 185) 

PRINT 4 1, "PULSA FPARA ACABAR” 

A$ = INKEY$ 

IFA$ = "E" ORA$ = “e” THEN CLOSE 4 1: SCREEN O : GOTO 170 
IFAS ="F” OR A$ = “f” THEN SCREEN O : END 

GOTO 540 


La instrucción 520 coloca el cursor en el punto (15, 175) de la pantalla de alta re- 
solución, a partir del cual la instrucción 530 imprime el texto. 


Las instrucciones 550 y 560 analizan la opción elegida. 


1300 
1310 
1320 
1330 
1400 
1410 
1420 
1430 
1500 
1510 
1520 
1530 
1540 
1600 
1610 
1620 
1630 
1640 


FOR | = 1 TO N 

Z (1) = X (1) 

NEXT | 

GOTO 300 

FOR | = 1 TO N 

Z (1) =X(1A2 

NEXT | 

GOTO 300 

FOR | = 1 TO N 

IF X (1) = O THEN PRINT “NO ES POSIBLE” : GOTO 170 
Z (1) = 1/X (1) 

NEXT | 

GOTO 300 

FOR | = 1TO N 

IF X (1) = O THEN PRINT “NO ES POSIBLE” : GOTO 170 
Z (1) = 1X(YA2 

NEXT | 

GOTO 300 
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Estos distintos bloques de instrucciones calculan los valores de la variable indepen- 
diente Z (1) en función de los datos X (I) introducidos al principio del programa. 


Dos ejecuciones del programa 


e Supongamos que las medidas registradas en un experimento son las indicadas 
en la tabla. 


Jos [5[=[5[0]5]2 
y] 06 [29/6210] 0052 5| 05056570 


Después de teclear RUN y | RETURN | se procede a introducir estos datos. 
Si se elige la opción 1 se obtiene el resultado que muestra la fotografía. 


Como esta gráfica no se aproxima a una línea recta, en lugar de pulsar F para 
acabar, se elige E para elegir la opción 2. Se obtiene este otro resultado: 
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Esta gráfica se aproxima mucho a una recta. Se deduce, entonces, que la rela- 
ción es del tipo Y = K * X?, siendo K la constante de proporcionalidad (*). 


e Supongamos, ahora que las medidas registradas en otro experimento son 
éstas: 


OOOO cane 
> 100]550[305[200]225/ 190155 140[125[ 150 


Después de teclear RUN y | ENTER | se procede a introducir los datos regis- 
trados en la tabla. 


Si se ejecutan las cuatro opciones (pulsar sucesivamente las teclas 1, 2, 3 y 4) 
se obtienen estas gráficas. 


o 


UNA o PARA ELEGIR Y F PARA Puga E PARA ELEGIR Y F PARA 


- 


CUBA E PORO ELEGIR Y F PARA PULGAR E PARA ELEGIR Y F PARA 


Observando estas gráficas, se llega a la conclusión que la relación entre X, Y es 
del tipo Y = K (1/X), siendo K la constante de proporcionalidad (**), 


(*) Esta constante se puede determinar dividiendo un valor de Y por el cuadrado de su corres- 
pondiente valor de X: Y = K + X? => K = Y/X? 


(**) Esta constante se puede determinar multiplicando un valor de Y por su correspondiente 
valor de X: Y = K * (1/X) > Y * X = K 
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12. Simulación de fenómenos 
aleatorios en biología y física 


1. Estudio de leyes genéticas 
Presentación del problema 


Vamos a analizar una de las leyes genéticas: la segunda ley de Mendel. 
Recordaremos esta ley mediante un ejemplo. 


Si se cruzan flores rojas (RR) con flores blancas (BB), en la primera generación sur- 
gen flores híbridas, de color rosa (RB) (*). Al cruzarse dos flores de color rosa de 
la primera generación, en la segunda generación pueden surgir estas posibili- 
dades: 


Ja 
generación 
posibilidades 

ZN 

generación 
roja rosa rosa blanca 
—__ "Y" 7 


Es decir, si se parte de dos individuos híbridos de la primera generación, en la se- 
gunda generación existe 25 % de posibilidades de que la flor sea roja, 50 % de 
que sea de color rosa y 25 % de que sea blanca. O, pasando a términos de fre- 
cuencia, de 100 000 individuos de esta segunda generación, unos 25 000 serán 
rojos, unos 50 000 serán de color rosa y unos 25 000 serán blancos. 


La comprobación experimental de esta ley lleva mucho tiempo y exige mucho tra- 
bajo. Pero con ayuda del microordenador se pueden simular los experimentos y 
comprobar rápidamente la ley. Para ello simularemos el lanzamiento de dos datos 
X, Y. Cada uno de ellos representa una flor híbrida, en cuanto que contiene, por 
decirlo así, dos genes diferentes: impar (R) y par (B). 


(*) RR y BB representan genotipos de dos razas puras, cada uno de ellos con dos genes R, R y 
B, B. El genotipo RB es híbrido, por tener dos genes diferentes R y B. 
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Al lanzar los dos dados pueden presentarse cuatro posibilidades: 


dado 1 dado 2 . 
lanzamiento 
E 
MS dados 
A 
j posibilidades 


25% 25% + 25% 25% 


A 
50% 
El programa 


El número N de «lanzamientos» a que hacen referencia las instrucciones 10 y 20 
no hay que entenderlo como el número total de cruzamientos que se van a simu- 
lar en el programa entre los individuos híbridos X e Y. El número total de cruza- 
mientos viene dado por el valor de la variable T (instrucción 205). 


10 PRINT “NUMERO DE LANZAMIENTOS” 


20 INPUT N 

25 T=0 

30 C=0 

40 D=0 

50 E=0 

60 FORI=1TON 

70 X = INT (RND (-TIME) + 6 + 1) 
80 Y = INT (RND (TIME) * 6 + 1) 
90 FORJ=1T0 3 
100 FOR K = 1 TO 3 


110 IFX=(2*)J— 1) AND Y = (2 * K— 1) THEN 180 
120 IF X =(2*J- 1) AND Y = 2 + K THEN 200 
130 IF X=2* J AND Y = (2 * K-— 1) THEN 200 
140 1F X= 2*J AND Y = 2 + K THEN 160 

150 GOTO 210 - 

160 C=C+1 

170 GOTO 210 

180 D=D+1 

190 GOTO 210 

200 E=E+1 

205 T=T+1 

210 NEXT K 

220 NEXT) 

230 NEXT | 


En los bucles anidados se simula el lanzamiento de un dado N veces y se obtienen 
los resultados X e Y de cada lanzamiento. 


En los bucles correspondientes a las instrucciones 90 a 220 se cuenta el número 
de veces que sale cada combinación, teniendo en cuenta que la combinación 
pares - impares se repite. 


240 PRINT C, C * 100/N 
250 PRINT D, D * 100/N 
260 PRINT E, E * 100/N 
270 END 


Finalmente, con estas instrucciones se imprimen los resultados C, D y E y los por- 
centajes que representan, respectivamente. 


Ejecución del programa 


Si después de teclear RUN y | RETURN | se introduce 1000 para N se obtienen 
estos resultados. 


2. Vida 


Planteamiento del juego 


Este es un juego muy conocido que simula el nacimiento, crecimiento y muerte 
de una colonia de células. 
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Las células iniciales se colocan al azar en una cuadrícula de 10 X 10 y nacen, vi- 
ven y mueren según las siguientes reglas: 


— Cada célula de la retícula tiene 7 vecinas. 

— Cada célula con 2 ó 3 vecinas sobrevive a la siguiente generación. 

— Si hay 3 y sólo 3 células próximas nace una nueva célula. 

— Cualquier célula con 4 o más vecinas muere por haber exceso de población. 


El programa 


A continuación se simula lo que les ocurre a diferentes generaciones de células 
hasta llegar a una situación estable (una situación puede consistir en la muerte de 
todas las células). 


10 CLS 

20 G=0 

30 DIM A (10, 10), B (10, 10) 
40 FOR X= 2T0 9 

50 FORY =2T09 

60 A = RND (-TIME) 

70 IFA > 0.35 THEN A (X, Y) = 1 
80 B(X,Y)=A(X, Y) 
90 NEXT Y 
100 NEXT X 
105 PRINT 


Con estas instrucciones se plantea una situación inicial. 


110 GOSUB 1000 
120 G=G+1 

130 FORX=2T09 
140 FORY=2T09 
145 C=0 


150 IFA(X-1,Y- 1)= 1THENC=C+1 
160 IFA(X- 1, Y) = 1THENC=C+ 1 
170 IFA(X-1,Y+1)= 1THENC=C+1 
180 IFA (X, Y - 1) =1THENC=C+ 1 
190 IFA (X, Y + 1) = 1THENC =C+ 1 
200 IFA(X+1,Y- 1) = 1THENC=C +1 
210 IFA(X+ 1, Y) = 1THENC=C+1 
220 IFA(X+1,Y+1)= 1THENC=C+1 


La instrucción 110 transfiere el control a una subrutina que dibuja la situación 
inicial, la instrucción 120 establece el paso de las generaciones y las siguientes, 
analizan la situación de cada célula. 


230 IFA (X, Y) = 1 AND C <> 3 AND C <> 2 THEN B (X, Y) = 0 
240 IFA (X, Y) = 0 AND C = 3 THEN B (X, Y) = 1 

250 NEXT Y 

260 NEXT X 

280 GOTO 105 


La instrucción 230 produce la muerte de una célula, mientras que la 240 produce 
su nacimiento. 


1000 LOCATE 9, 3 : PRINT “GENERACION”; G 

1010 FOR X= 1TO 10 

1020 FOR Y = 1 TO 10 

1030 A (X, Y) = B (X, Y) 

1040 IF A (X, Y) = 1 THEN LOCATE X + 10, Y + 4 : PRINT “0” 
1050 IF A (X, Y) = O THEN LOCATE X + 10, Y + 4: PRINT*” ” 
1060 NEXT Y 

1070 NEXT X 

1080 RETURN 


Esta subrutina construye la figura adecuada a las distintas situaciones. 


Ejecución del problema 


Tecleando RUN y | RETURN | se obtiene este resultado. 
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3. Simulación de la desintegración radiactiva 
Planteamiento del problema 


Supongamos una determinada cantidad de monedas que lanzamos al aire. Una 
vez que éstas han caído, quitamos las que presentan «cruz» anotando las 
«caras». Lanzamos al aire las monedas restantes (las «caras») retirando las 
«cruces» y anotando de nuevo las que han salido «cara». Repetimos el procedi- 
miento varias veces. 


La variación en cada lanzamiento del número de caras puede utilizarse como mo- 
delo para analizar la variación en el tiempo del número de átomos que no se de- 
sintegran en un proceso radiactivo. 


El programa 


El programa que se incluye a continuación hace esta simulación con la ayuda de 
un microordenador. Permite obtener el número de átomos que no se han desinte- 
grado todavía en cada instante, y pulsando CONT se obtiene la representación 
gráfica de este número en función del tiempo. 


2 CLS 
5 DIM N (500), H (500), V (500), D (500) 
10 INPUT “NUMERO INICIAL DE MONEDAS”; N 
20 D(0)=0:J=0:N(0)=N 
30 FOR!|=1TON ()) 
40 S= INT (RND (-TIME) + 0.5) 
50 IF S = 0 THEN GOTO 70 
60 D(J) =D (Y) +1 
70 NEXT | 
80 N(+1)=N ()) - D ()) 
90 IFN (J) > = 2 THEN J) =J+ 1: GOTO 30 


Estas instrucciones simulan el lanzamiento de N monedas, asignando el valor 1 a 
las caras y el O a las cruces. En cada lanzamiento se contabilizan las caras. 


A continuación se plantea un nuevo lanzamiento retirando las caras que se ha- 
bían obtenido anteriormente (instrucción 80). 


100 FOR U=0T0O J 
110 PRINT U, N (U) 
120 NEXT U 

130 STOP 
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La instrucción 110 imprime en la pantalla el número de átomos sin desintegrar (o 
monedas que presentan cruz), en función del tiempo. 


140 DX=J-0:DY=N(0) — N ()) 

150 XS = 245/DX : YS = 180/DY 

160 SCREEN 2 

170 FORU=0T0J 

180 H(U) = U*XS : V (U) = (N (U) + N (J)) + YS 
190 NEXT U 

200 PSET (H (0), 175 — V (0)) 

210 FORU=1T0J 

220 LINE (J (U — 1), 175 — V (U — 1)) — (H (U), 175 — V (U)) 
230 NEXT U 

240 GOTO 240 


Una vez que se pulsa CONT se obtiene la representación gráfica del número de 
átomos sin desintegrar en función del tiempo. 
Ejecución del programa 


Si después de teclear RUN y | RETURN | se introduce 5000 para N se obtiene esta 
pantalla. 


4. Equilibrio térmico de un sólido 


Presentación del problema 


El equilibrio térmico de un sólido puede simularse con el microordenador me- 
diante el siguiente juego: 
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Se construye un tablero de 800 casillas, por ejemplo, y en cada una de ellas se 
coloca una ficha. 


Se lanzan dos dados hipotéticos de 800 caras. El resultado del lanzamiento del 
primer dado señalará de qué casilla se retira una ficha, mientras que el del se- 
gundo, señalará el lugar en el que debe colocarse dicha ficha. 


Si la casilla obtenida en el primer lanzamiento no contiene ficha, se repite el lan- 
zamiento, pero contabilizándolo entre los realizados. 


Si se hace un número suficiente de lanzamientos la distribución del número de 
casillas que contienen 0, 1, 2, 3, 8 fichas es exponencial. 


Los resultados concuerdan en cierta manera con los que se obtienen en la distri- 
bución cuántica de un sólido en equilibrio térmico. 


El programa 
A continuación se simula este juego. 


5 KEY OFF 
10 CLS 
20 DIM C (8), H (8), V (8) 
30 INPUT “CUANTOS MOVIMIENTOS”; A 
40 X= 50000 
50 FOR!|=0T0 800 
60 POKE X + 1, 49 
70 NEXT | 


Con la instrucción POKE se coloca un cuanto en las 800 posiciones de memoria 
siguientes a la 50 000 (el código ASCII del número 1 es precisamente 49). 


80 GOSUB 2000 


Mediante esta subrutina se simula el lanzamiento de un dado. 


90 Z = PEEK (X + D) 

100 IF Z< 49 THEN GOTO 80 
110 POKEX+D,Z-1 

120 GOSUB 2000 

130 B = PEEK (X + D) 

140 POKEX+0D,B +1 

150 J=J+1 

160 IF J= A THEN GOTO 1000 
170 GOTO 80 
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En la instrucción 90 se coge una ficha de una casilla elegida al azar (subrutina 
2000). Si la casilla no contiene ficha se repite el lanzamiento (instrucción 80). 


Después de simular de nuevo el lanzamiento del dado (instrucción 120), de 
acuerdo con el resultado que salga se añade una ficha a la posición obtenida (ins- 
trucción 140). 


El proceso se repite hasta alcanzar el valor de A (número de lanzamientos). 


1000 FORT = 0 TO 800 

1010 FOR!|=0T0 8 

1020 IF PEEK (X + T) = 48 + | THEN C (1) = C (1) + 1 
1030 NEXT | 

1040 NEXT T 


Con este bloque de instrucciones se contabilizan las casillas que contienen 0 fi- 
chas (48 es el código ASCII de 0), 1 ficha (49 es el código ASCII de 1), y así suce- 
sivamente hasta 9 fichas (48 + 8 es el código ASCII de 8). 


1100 FOR!I=0T08 

1110 PRINT “EL NUMERO DE ATOMOS CON”; |; “CUANTOS ES”; C (1) 
1120 NEXT | 

1130 STOP 


Estas instrucciones imprimen en la pantalla el número de átomos con 0, 1, 2, ..., 
8 cuantos. 


1200 Y2 = C (0) : Y1 = C (0) 

1210 FOR!|=1T0 8 

1220 IF Y2 < C (1) THEN Y2 = C (1) 
1230 IF Y1 > C (1) THEN Y1 = C (1) 

1240 NEXT | 

1250 IF Y1 > O THEN Y1 =0 

1260 DX=8:DY=Y2-Y1 

1270 XS = 245/DX : YS = 175/DY 

1280 SCREEN 2 

1290 FORI=0T0 8 

1300 H(I) =1*XS: V (1) = (C (1) + ABS (Y1)) * YS 

1310 NEXT | 

1320 PSET (H (0), 175 — V (0)) 

1330 FOR!|=1T08 

1340 LINE (H (1 — 1), 175 — V (1 — 1)) — (H (1), 175 — V (1) 
1350 NEXT | 

1360 GOTO 1360 
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Utilizando un procedimiento que se ha usado varias veces a lo largo de este libro, se 
representa el número de casillas en función del número de cuantos que hay en di- 
chas casillas. La gráfica aparece después de pulsar la tecla CONT. 


2000 D = INT (800 * RND (—TIME)) 
2010 RETURN 


Ésta es la subrutina que elige al azar la casilla. 


Ejecución del programa 


A pesar de utilizar instrucciones POKE y PEEK, la ejecución del programa es bas- 
tante lenta, debido al gran número de casillas que elige y al número de movi- 
mientos que realiza. Así, para 5000 movimientos, la ejecución puede durar 
aproximadamente unos diez minutos. Si se introduce, entonces, 5 000 para A, 
después del tiempo indicado se obtiene este resultado. 
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13. Cálculo numérico en Física 
y Biología 


1. Análisis del concepto de velocidad instantánea (derivada de 
una función en un punto) 


Presentación del programa 


Consideremos un móvil que se mueve sfguiendo una trayectoria cualquiera. Su 
velocidad media entre dos instantes t1 y t2, en los que ocupa dos posiciones y1 
e y2 (con respecto a un punto de referencia 0), se obtiene dividiendo el espacio 
recorrido por el intervalo de tiempo que tarda en recorrerlo: 


vm = (y2 — y1) / (t2 — t1) 
yy 
Ya 
t, : 
0) 5 
Otra velocidad que se puede considerar en este movimiento es la velocidad ins- 
tantánea (la que marcaría el velocímetro de un coche), que puede obtenerse ha- 


llando el límite al que tienden la sucesión de velocidades medias cuando los ins- 
tantes t1 y t2 tienden a confundirse en uno solo. 


La sucesión se expresa así: 


El programa 


En el programa que se incluye a continuación se realiza este proceso para com- 
probar finalmente que los valores de las velocidades medias se aproximan a uno 
determinado, que es la velocidad instantánea (o derivada de la función que da el 
espacio con respecto al tiempo en ese instante). 
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El procedimiento permite calcular, con aproximación suficiente, la derivada de 
cualquier función en un punto. 


1 REM VELOCIDAD INSTANTANEA 
10 CLS 
20 DEFFNS(T=T7 2+2+*T+3 
30 INPUT “INTRODUCE EL INSTANTE”; T1 


En estas instrucciones se define la función que relaciona el espacio con el tiempo 
y se introduce el instante en el que se va a obtener la velocidad. 


35 CLS 
40 FORN=1T020 

50 T=T1+05 7 N 

60 V=(FENS(T) - ENS(T1)/T- T1 
70 PRINT TAB (1); T; TAB (17); V 

80 NEXT N 

90 END 


En esta última parte se van obteniendo los términos de la sucesión de velocidades 
medias. Las distintas aproximaciones de T a T1 se obtienen sumando a T1 los va- 
lores cada vez más pequeños obtenidos al elevar 0.5 a los diversos valores de N. 


El procedimiento puede aplicarse a cualquier función matemática que se desee, 
por ejemplo, a y = el, y = e* y = 35 ent y a otras más complicadas. 


Ejecución del programa 


Si después de teclear RUN y | RETURN | se introduce el valor 5.7 para T1, se ob- 
tiene este resultado. 


2. Cálculo del espacio recorrido en un movimiento 
(Cálculo numérico de integrales definidas) 


Presentación del problema 


En un movimiento rectilineo y uniforme, el espacio recorrido entre dos instantes 
a y b viene expresado por la fórmula 


S =V (b - a) 


donde S es el espacio recorrido, V es la velocidad (constante) y b — a, el intervalo 
de tiempo transcurrido. 


En un diagrama velocidad — tiempo el espacio recorrido corresponderá al área del 
rectángulo de altura v y base b — a. 


Si el movimiento no es uniforme, la gráfica velocidad — tiempo no es una línea 
recta paralela al eje t, con lo que el cálculo del espacio recorrido se hace más 
complejo. Para hacer este cálculo se divide el intervalo a — b en intervalos par- 
ciales c — d, y en cada uno de estos intervalos se aproxima el movimiento real a 
un movimiento uniforme con velocidad constante e igual a la del móvil real en la 
mitad del intervalo. 
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El espacio recorrido en cada intervalo c — d, en el movimiento ficticio, es igual al 
área del rectángulo destacado en color. 


Me d 


La suma de los espacios recorridos en todos los intervalos c — d se aproximará 
tanto más al espacio total recorrido en el movimiento real cuanto menor sea la 
amplitud de los intervalos c — d, lo cual equivale a que estén más próximos entre 
sí los tramos escalonados. 


Formando una sucesión de sumas, así calculadas, para intervalos de tiempo cada 
vez más pequeños, dicha sucesión tiende a un valor límite: el espacio recorrido 
por el móvil en el movimiento, o sea la integral definida de la velocidad en el 
intervalo de tiempo a — b. 


El programa 


Aplicando los conceptos que se acaban de exponer, se logra calcular la integral 
definida con suficiente aproximación. 


1 REM ESPACIO RECORRIDO 
10 DEFFNV(T)=2*T 7 2+3 
20 INPUT “LIMITE INFERIOR”; A 
30 INPUT “LIMITE SUPERIOR”; B 
40 LETD=B-A 
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En esta parte del programa se define la función velocidad y se introducen los 
tiempos entre los cuales se va a calcular el espacio recorrido, A y B. 


Se define también una variable D, que se inicia con el valor B — A y que en el 
curso de la ejecución del programa se irá dividiendo por 2 para obtener su mitad. 


50 CLS 

60 FOR!I= 1 TO 10 

70 LETM=0  ' 

80 FORT =A + D/2 TO B STEP D 
90 LETM=M + EFN V (T) 
100 NEXT T 

110 LETM=M*D 

120 PRINTM 

130 LET D = D/2 

140 NEXT | 

150 END 


Este bloque de instrucciones calcula diez términos de la sucesión de sumas, las 
cuales se van acumulando en la variable M. Cada suma se obtiene en el bucle ¡n- 
terno (instrucciones 80 a 100) y en la instrucción 110. 


Para cada valor de | que determina el bucle externo (60 a 140), se divide D por la 
mitad. | 


El último valor de M que se imprime en la pantalla es el más aproximado al espa- 
cio recorrido por el móvil real. 


Ejecución del programa 


Tecleando RUN y | RETURN | e introduciendo el valor 2 para A y el valor 15 para 
B, se obtiene este resultado. 
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3. Cálculo del número de células nacidas en un determinado 
tiempo 


Presentación del problema 


El número de células de una población puede duplicarse si dichas células se repro- 
ducen dividiéndose en dos, y siempre que el alimento sea el adecuado. 


., ON e 
La relación a N, que expresa que la variación del número de células por uni- 


dad de tiempo es igual al número de células existentes en un instante (número 
inicial de células), representa bien este fenómeno. 


Escribiendo = = N en forma de incrementos finitos se obtiene: 


AN 


AN _ on N(t+ At) — N(t) _ 
2 =N - obien At = N 


Despejando N (t + At) se consigue esta ecuación: 
N (t+ At) =N (t) x At+N (t) 
Si, por ejemplo, N (0) = 1 y At = 1, se tiene: 


N(1)=N(0)* 1+N(0)=1+1=2 
N(2)=N(1)*1+N(1)=2+2=4 
N(B)=N(2)*1+N()=4+4=8 


Se ve cómo, dando valores a t, el número de células se duplica en cada unidad de 
tiempo. Así: 

— al término de la primera unidad de tiempo el número de células es 2; 

— al término de la segunda unidad de tiempo el número de células es 4; 

— al término de la tercera unidad de tiempo el número de células es 8; 

— etcétera. 


El programa 
El programa que se transcribe a continuación imprime el número de células que 


hay en cada instante y después de teclear CONT dibuja la gráfica número de cé- 
lulas — tiempo. 
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10 DIMT (100), N (100), H (100), V (100) 

20 CLS 

30 INPUT “NUMERO INICIAL DE CELULAS”; NO 
40 INPUT “TIEMPO”; TM 


En estas instrucciones se introducen los datos necesarios para que el programa 
pueda ejecutarse. 


50 E=0:3N.(0)=NO:M=1:T10)=.0 

60 N(I+1)=N()x*IT+N (1) 

70 T(+1)=T()+IT 

80 IFT(1+ 1) <= TMTHEN | = 1 + 1: GOTO 60 


Este grupo de instrucciones determina el número de células existentes en un ins- 
tante a partir de las existentes en el instante anterior. En la instrucción 70 se in- 
crementa el tiempo en IT y se comprueba si se ha alcanzado el tiempo previsto. 


90 FORJ=0T01+ 1 
100 PRINTT (3), N ()) 
110 NEXT J 
120 STOP 


Mediante estas instrucciones se imprimen en la pantalla los tiempos y el número 
de células. 


En la línea 120 se detiene la ejecución. Para que la ejecución continúe hay que te- 


clear CONT : 


130 X2=T(0): X1 =T(0): Y2 =N (0): Y1 = N (0) 
140 FORJ=1TO1+ 1 

150 IF X2 <T (J) THEN X2 = T ()) 
160 IF X1 >T (J) THEN X1 =7T (0) 
170 IF Y2<N (J) THEN Y2 = N ()) 
180 IF Y1 >N (J) THEN Y1 =N ()) 
190 NEXT J 

200 IFX1> O THEN X1 = 0 

210 IF Y1>O THEN Y1 =0 

220 DX=X2- X1: DY = Y2 — Y1 
230 XS = 245/DX : YS = 175/DY 
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Cuando se escribe CONT se buscan los valores mínimos y máximo de T ( ) yN ( ), y se 
establecen las escalas adecuadas para dibujar la gráfica N — T. 


240 SCREEN 2 

250 FORJ=0TO1+1 

260 H (1) = (T (1) + ABS (T (0))) + XS : V (1) = (N (1) + ABS (N (0))) + 
YS 

270 NEXT) 

280 PSET (H (0), 175 — V (0)) 

290 FORJ=1TO1+1 

300 LINE (H (J — 1), 175 — V (Y — 1)) — (H (9), 175 — V (1) 

305 NEXT) 

310 GOTO 310 


Este bloque de instrucciones calcula y dibuja los puntos, y a continuación los une 
mediante tramos rectos. 


Ejecución del programa 


Si después de teclear RUN y | RETURN | se introduce 100 para N y el valor 18 
para T, se obtiene en la pantalla este resultado. 


En la línea 120 la ejecución del programa está detenida. Tecleando CONT y 
RETURN | sigue la ejecución. Al final de la misma se obtiene en la pantalla la 
gráfica N — T. 


Observación 


El programa que se acaba de ejecutar se puede utilizar, con ligeras modifica- 
ciones, para estudiar otros procesos, como la desintegración radiactiva, la des- 
carga de un condensador, etc.; es decir, aquellos en que dos magnitudes F y X 
están relacionadas por una ecuación del tipo 


dF 


- 


siendo K una constante. 


4. Un problema de ecología 


Presentación del problema 


Estudiaremos ahora la evolución de dos especies de animales que coexisten en un 
mismo hábitat, alimentándose una especie de la otra. 


Un ejemplo de estas características lo constituye una población de conejos y otra 
de zorros, donde los segundos se alimentan de los primeros. 


Si x es el número de conejos e y el número de zorros, se admite que el tamaño 
de las dos poblaciones evoluciona según estas leyes: 


de = Ax — Bxy 
Y Zo - 


Si no hay zorros, los conejos aumentan proporcionalmente a su número Ax y si 
no hay conejos, los zorros se morirían de acuerdo con la expresión — Dy. Los tér- 
minos en xy reflejan la interrelación entre ambas especies, y representan un creci- 
miento en la población de zorros y una disminución en la de conejos. 


En forma de incrementos finitos, las relaciones anteriores pueden escribirse así: 


"25 A 
MEA Ax — Bxy 
Ay Ñ Ml 
Mt Cxy — Dy 


Despejando x (t + A t) e y (t+ A 1), se obtiene: 


Xx (t + A t) = (Ax (t) — Bx (t) : y (t)) At + x (t) 
y (t + At) = (Cx (t) y (t) — Dy (t)) At + y (t) 
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El programa 


Mediante estas expresiones se puede averiguar la población X de conejos y la po- 
blación Y de zorros en cada instante. 


Esto es lo que calcula el siguiente programa, que además dibuja en una misma 
pantalla (después de teclear CONT y | RETURN |) las gráficas X — T e Y — T, lo 
que permite comparar la evolución de ambas especies. 


En este programa se eligen como valores adecuados los siguientes valores para las 
constantes: 


A=1,B=2C=2yD=3. 


1 REM ECOLOGIA 
10 DIMT (500), X (500), Y (500), H (500), V (500), Z (500) 
DO ASAS ZAS 
30 CLS 
40 INPUT “NUMERO INICIAL DE CONEJOS”; XO 
50 INPUT “NUMERO INICIAL DE ZORROS”; YO 
60 INPUT “TIEMPO”; MT 
70 IN=0.0X10)=X0.5Y:10).= Y0:T-10)=.00 0 =0 


Aquí se plantean las condiciones iniciales y se introducen los datos necesarios 
para que el programa pueda ejecutarse. 


80 X(1+ 1) =(A*X (1) — B * X (1) x Y (1)) x IT + X (1) 
90 Y (1+ 1) = (C * X (1) * Y (1) — D * Y (1) + IT + Y (1) 
100 T(1+1)=T (1) + IT 

110 IFT(I +1) <= MTTHEN | =1+ 1: GOTO 80 


Estas instrucciones se calculan los valores X e Y a partir de sus valores anteriores. 
El tiempo aumenta en intervalos IT hasta que se alcance el valor MT. 


120 FORJ=1TO!I+ 1 

130 PRINT TAB (1); X (1); TAB (16); Y (3) 
140 NEXT J 

145 STOP 


El bucle 120-130 imprime en la pantalla los valores sucesivos de X e Y, y la ins- 
trucción 145 detiene la ejecución del programa. 
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150 


160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 


X= M0) Y2=V 101: 712.2 TF (0) XT =X10)<1 ="W10): 
T1=T(0) 

FOR J = 17101 +1 

IFT2 <T (3) THEN T2 =7 ()) 


IF T1 > T (J) THEN T1 =7 ()) 

IF X2 < X (J) THEN X2 = X ()) 

IF X1 > X (1) THEN X1 = X (9) 

IF Y2 < Y (1) THEN Y2 = Y ()) 

IF Y1 > Y (3) THEN Y1 = Y ()) 

NEXT J 

IF X1 > 0 THEN X1 =0 

IF Y1 > 0 THEN Y1 =0 

IFT1 > O THEN T1 =0 

DT =T2 — T1: DX = X2 — X1 : DY = Y2 — Y1 
TS = 245/DT : XS = 180/DX :'YS = 245/DY 


Después de escribir CONT se buscan los valores máximo y mínimo de las distintas 
variables y se establecen las escalas para dibujar. 


290 
300 
310 


320 
330 
340 
350 
360 
370 
380 
390 
400 
410 


SCREEN 2 

FORJ=0T0O1+ 1 

H (3) = (T (1) + ABS (T (0))) * TS : V (1) = (Y (1) + ABS (Y (0))) * 
YS : Z (1) = (X (1) + ABS (X (0))) * XS 

NEXT J 

PSET (H (0), 175 — V (0)) 

FOR J = 1T01+ 1 

LINE (H (J — 1), 175 — V QU) — 1)) — (H (9), 175 — V (9) 
NEXT J 

PSET (H (0), 175 — Z (0)) 

FOR J=1T01+ 1 

LINE (H (J — 1), 175 — Z (Y — 1)) — (H (9), 175 — 2 (1) 
NEXT J 

GOTO 410 


Este conjunto de instrucciones dibuja en una misma pantalla las gráficas número 
de conejos — tiempo y número de zorros — tiempo. 


Los valores iniciales de conejos XO y los valores iniciales de zorros YO, así como el 
tiempo MT, deben ser números pequeños para que no se produzca error (OVER 


FLOW). 
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Ejecución del programa 


Después de teclear RUN y | RETURN | introducir los valores que se indican a con- 
tinuación para las correspondientes variables: 

XO = 3, YO = 0.5, MT = 4.2 (*) 
En la pantalla se obtiene este resultado: 


En la línea 145 se detiene la ejecución del programa. Para que continúe hay que 


teclear CONT | RETURN |, obteniéndose las gráficas que se observan en esta pan- 
talla. 


(*) Los valores 3 y 0.5 no indican que se parta inicialmente de 3 conejos y de 0.5 zorros, sino que 
la relación inicial del número de conejos a zorros es 3 a 0.5. El tiempo MT = 4.2 tampoco es el 
real, sino el que corresponde a una cierta escala. 
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5. Análisis del movimiento de un cohete de masa variable 


Consideremos un cohete cuya masa en un instante es m y que, quemando una 
masa Am durante un tiempo At, aumenta su velocidad en Av. Supongamos 
también que los gases procedentes de esta combustión tienen una velocidad u 
con respecto al cohete. 


La fuerza media que actúa sobre el cohete debe ser igual a la variación de la can- 
tidad de movimiento del combustible expulsado. 


Es decir, 
mA vAt=-v AmAt 


Si la masa inicial del cohete es My y el combustible se quema a velocidad cons- 
tante y = Am/At, después de un tiempo t, la masa será: 


m= Mo — pt 
y sustituyendo en la ecuación anterior, se obtiene: 
(Mo — put) E = Up 
O bien 
A o 
At Mo — pt 


Si la masa inicial del cohete es My = 10000 kg, la velocidad de combustión de 
los mismos es y = 50 kgs”, y la velocidad de éstos con respecto al cohete es 
u = 2000 ms”', la expresión se convierte en 


Av 2000-50 
At 10000 — 50t 


Si el cohete está despegando de la Tierra y si suponemos que no se aleja excesi- 
vamente de ésta como para considerar variable la aceleración de la gravedad, po- 
demos escribir: 


Av _ 100000 
At “10000 — 50t 


Teniendo en cuenta que — se puede escribir como 


v(t+ At) — v (t) 
At 


9,8 
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despejando, se obtiene: 


100000 


v(t+ Ad = v() + (50000 - 50t 


- 98) At 


La altura viene dada por 
y (t+ At) =y()+v() At 


Si la velocidad inicial es V (0) = 0, considerando intervalos de tiempo de un se- 
gundo, la velocidad en la mitad del primer intervalo es 


100000 
10000 — 50t 


y la altura alcanzada en el primer segundo viene dada por 
Y (1) = V (0.5) x 1 + y (0) 
De la misma manera se obtienen las expresiones v (1.5) e y (2): 


100000 
10000 — 50t 


y (2) = v (1.5) x 1 + y (1) 
y análogamente se procede para los sucesivos valores del tiempo. 


v (0.5) = V (0) + ( — 9.8) * 0.5 


v (1.5) = V (0.5) + ( — 9.8) 


El programa 


El programa siguiente permite obtener una representación gráfica de la función 
que da la altura alcanzada en función del tiempo. 


5 CLS 

10 DIM T(100), Y(100), H(100), V(100) 

20 INPUT “ALTURA”; YO á 

30 CLS 

20- "TO =0211=+1:MTF= 20. V=0:1=:0 

50 V= (100000 / (10000 — 50 x IT) — 9.8) * IT/2 + V 
60 Y (+1) =Y() + Vx*IT 

70 V=V+ (100000/(10000 — 50 * T) — 9.8) * IT + V 
80 T(1+1)=T() +17 

90 IFT(I+ 1) <= MT THEN | =1+ 1: GOTO 40 


Mediante estas instrucciones se introducen los datos y se obtienen los valores de 
la altura a partir de las alturas y velocidades anteriores. El tiempo se incrementa 
en intervalos | hasta llegar al límite MT (que aquí se ha fijado en 20). 
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100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
205 
210 
220 


230 
240 
250 
260 
270 
280 


X2 = T (0) : X1 =7T (0) : Y2 = Y (0) : Y1 = Y (0) 
FOR J=1TO1+ 1 

IF X2 < T (3) THEN X2 = T ()) 
IF X1 > T (3) THEN X1 =T ()) 
IF Y2 < Y (J) THEN Y2 = Y (0) 
IF Y1 > Y (1) THEN Y1 = Y ()) 
NEXT J 

IFX1 > O THEN X1 =0 

IF Y1 > 0 THEN Y1 =0 

DX = X2 — X1 : DY = Y2 — Y1 
XS = 245/DX : YS = 175/DY 


SCREEN 2 

FORJ=0TO1+ 1 

H (3) = (T (1) + ABS (T (0))) + XS : V (3) = (Y (1) + ABS (Y (0))) * 
XS ] 

NEXT J 

PSET (H (0), 175 — V (0)) 

FOR J=1TO1+ 1 

LINE (H (J — 1), 175 — V (Y — 1)) — (H (9), 175 — V (1) 


NEXT J 
GOTO 280 


Con estas instrucciones se seleccionan los valores máximo y mínimo de las varia- 
bles, se establece la escala y se realiza el dibujo uniendo los puntos con tramos 


rectos. 


Ejecución del programa 


Si después de teclear RUN y | RETURN | se introduce el valor 100 para YO, se ob- 
tiene esta gráfica. 
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6. Estudio de las oscilaciones en un circuito eléctrico de corriente 
alterna 


Presentación del problema 


Consideremos un circuito con un condensador C, una resistencia R y una autoin- 
ducción L, dispuestos en serie. 


R 


“0 
1 2 


El voltaje V (t) aplicado en los bornes 1 y 2 es igual al trabajo necesario para tras- 
ladar la unidad de carga desde el punto 1 al 2, a través de todos los elementos 
del circuito. Este trabajo se obtiene sumando los siguientes valores: 


— La caída de potencial en los extremos de la bobina: 


— La caída de potencia en los extremos de la resistencia: 
Va=RI=RdQ/dt 
— La caída de potencial en los extremos del condensador: 
Ve = Q/C 
Luego V (t) viene dado por 
V (t) = L d? Q/d t* + Rd Q/d t + Q/C 
Si el voltaje varía sinusoidalmente en función del tiempo, se tiene 
V (t) = Vo sen (2 n Ft) 


donde V, es el voltaje máximo, F la frecuencia y t el tiempo transcurrido. 
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Reemplazando V (t) por Vy sen (2 mn F t) en la ecuación anterior se obtiene esta 
expresión: 


dl Oe. 
L na + Rl + E = Vo sen (2 nm F T) 
. dl Busta 
Despejando 41 5 tiene: 
e = Vol sen (2 m1 Ft) — RI/L — QULO 


Y considerando incrementos finitos se tendrá: 
A WA t = VyL sen (2 n Ft) — RI/L — Q/LC 
O bien, teniendo en cuenta que 
SAN MEE = HH 
At At 
se obtiene: 
I(t+ A t) =1 (t) + (VyL sen (2 n Ft) — RIL — Q/LC) A t 


Pero como | = A Q/ A t, se llega finalmente a la expresión: 
Q(+A0=Q()+1:At 

Teniendo en cuenta las expresiones anteriores se calcula la carga para intervalos 

de tiempo sucesivos y la intensidad en la mitad de dichos intervalos. Así, por 

ejemplo, para los intervalos 0 — 0.1 y 0.1 — 0.2, se obtiene: 


| (0.05) = 1(0) + VyL * sen (2 rn FT x 0.1) — (RI(OL) - (0.1) — Q(O/LC) x (0.1) 


y 
Q (0.1) = Q (0) + 1 (0.05) x (0.1) 


| (0.15)=1 (0.05) +VyL sen (2 ni F tx0.1)—(RI (0.05)//L)x0.1—Q (0.05/LC) x 0.1 


Q (0.2) = 1 (0.15) x (0.1) + Q (0.1) 


El programa 


El programa tiene en cuenta los conceptos que se acaban de exponer y dibuja la 
gráfica intensidad — tiempo. 
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CLS 
INPUT “RESISTENCIA”; R 

INPUT “INDUCTANCIA”; L 

INPUT “CAPACIDAD”; C 

INPUT “VOLTAJE EFICAZ"; V 
INPUT “FRECUENCIA”; F 
W=2+*3,14 + F 

P=1/f 

XC = 1/(C * W) 

XL =L*W 

Q=0 

Z =SQRÍR 1 2 + (XL — XC) 1 2) 
VO = V * SQR (2) 

lO = VO/R 


Con estas instrucciones se determinan, a partir de los datos introducidos, los va- 
lores de la impedancia Z, del voltaje máximo VO en función del voltaje eficaz 
(Vmáx. = Ver - W2) y de la intensidad máxima. 


120 
125 
130 
140 
150 


160 
170 
173 
176 
180 
190 
200 


210 
220 
250 


T11=0:712=3*P. 

Y1 =-10 : Y2 =10 

DT = (T2 — T1)/245 : DY = (Y2 - Y1)/180 

SCREEN 2 

| = VO x* SIN (2 x 3.14 * Fx T) + DT/(2*L)-—Rx*IxDT/(2*L) - 
-Q*DT/(2xL*C)+1 

X3 = INT ((T - T1)/DT + 0.5) 

Y3 = INT ((l - Y1)/DY + 0.5) 

LINE (0,90) — (250,90) 

LINE (15,0) — (15,180) 

PSET (15 + X3, 180 — Y3) 

Q=0Q+Q+*DT 

| = VO + SIN (2 x 3.14 * Fx 1) x DT/L—Rx*I|*DT/L-Q* 
DT/(L*C) +! 

T=T+DT 

IF T < = T2 THEN GOTO 160 

GOTO 250 


En esta parte del programa se van trazando los distintos puntos, obtenidos por el 
procedimiento recursivo que ha sido descrito en los comentarios que anteceden al 
programa. Las instrucciones 175 y 176 trazan los ejes de coordenadas. 
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Ejecución del programa 


Proponemos se ejecute el programa para dos situaciones: la primera con valores 
cualesquiera de R, L, C, V y F, y la segunda, introduciendo el valor de F que co- 


a ¡ 
rresponda a la resonancia del circuito, es decir, para F = 67 VLC. 


e Primera ejecución 


Después de teclear RUN y | RETURN | si se introducen los valores: 


300 para R 

98 para L 

0.5 para C 

100 para V 

5 para F 

se obtiene esta gráfica: 


Segunda ejecución 


En esta ejecución introducimos estos valores: 


300 para R 

98 para L 

0.5 para C 

100 para V 

y para F, la frecuencia de resonancia: 


A A A 
3 ias 231487 00023 
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La gráfica que se obtiene es la que muestra la pantalla. 
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14. Geografía e Historia 


1. Repaso de capitales de países 


Presentación del problema 


Se puede utilizar el microordenador para repasar ciertos conocimientos, como, por 
ejemplo, capitales de países, ríos, montes y otros accidentes geográficos. También 
puede utilizarse en Literatura para repasar autores y obras, en Química para repa- 
sar nombres y símbolos de elementos químicos, etc. 


El programa 


En primer lugar se introducen en la memoria, como. variables alfanuméricas con 
índices, los países y capitales: P$ (1), C$ (1). 


Para dar cada contestación se introduce otra variable alfanumérica, B$. El número 
N de países que se introducen varía a voluntad. 


De todos los países y capitales se extraerá una al azar escogiendo su índice (5). 


La respuesta se comparará con el contenido de la variable C$ (l). En caso de error 
se dará el resultado correcto, y en caso de acierto se sumará una unidad a un 
contador. 


1 REM AVERIGUA CAPITALES 
10 Z = RND (TIME) 
15 INPUT N 
20 DIM P$ (N) 
25 DIM C$ (N) 


La instrucción 10 hace que la instrucción 90 genere números aleatorios. Si el mi- 
croordenador no es MSX la instrucción 10 se deberá suprimir o, en todo caso, ha- 
brá que sustituirla por otra, según el tipo de microordenador que se esté utili- 
zando. 


La instrucción 15 pide el número de países y capitales, y las instrucciones 20 y 25 
dimensionan las variables de cadena P$ y C$. 
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30 FORI=1TON 
40 INPUT P$ (1) 

50 INPUT C$ (1) 

60 NEXT | 


Hasta aquí se han introducido los N países y N capitales en la memoria del mi- 
croordenador. 


70 LETT=0 
La instrucción 70 inicializa un contador. 


80 FORJ=1T0 10 

90 S= INT (RND (1) * N + 1) 

100 PRINT “CAPITAL DE”; P$ (S) 
110 INPUT B$ 

120 PRINT B$ 

130 IF B$ = C$ (S) THEN GOTO 160 
140 PRINT “TE EQUIVOCASTE; EL RESULTADO ES:”; C$ (S) 
150 GOTO 180 

160 T=T+1 

170 PRINT “ACERTASTE” 

180 NEXT J 


En esta parte se obtiene un ndice S al azar y se compara la capital correspon- 
diente con la respuesta dada. Con el bucle 80-180 se repite el proceso 10 veces. 


190 PRINT “EL NUMERO DE ACIERTOS ES”; T 
200 END 


Se imprime finalmente el número de aciertos (contenido del contador T). 


2. Pirámides de población 
Presentación del problema 


En Geografía se plantea frecuentemente el problema'de la distribución de la po- 
blación de un país por edades y sexo. La representación gráfica de esta distribu- 
ción es lo que llamamos pirámide de población. 


Así, por ejemplo, la población de Guatemala en el año 1960 era de 3.965.000 
habitantes aproximadamente, y su distribución por edades, de cinco en cinco, era 
la siguiente: 
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Grupo de edades Número de hombres Número de mujeres 
(años) (miles) (miles) 


0a4 
5a9 
10a14 
15a19 
20 a 24 
25a29 
30 a 34 
35 a 39 
40 a 44 
45 a 49 
50 a 54 
55 a 59 
60 a 64 
65 a 69 
710274 
75 y más 


El total de hombres era igual a 1.903.000 y el total de mujeres 2.062.000. A par- 
tir de estos datos se calculan las frecuencias relativas correspondientes, dividiendo 
el número de hombres (o mujeres) de un grupo determinado de edades entre el 
total de hombres (o mujeres). La tabla de porcentajes que se obtiene es ésta: 


Grupos de edades Frecuencia relativa Frecuencia relativa 
(años) (hombres) (mujeres) 


0a4 
520 
10a 15 
15 a 19 
20a 24 
25a29 
30 a 34 
35 a 39 
40 a 44 
45 a 49 
50 a 54 
55 a 59 
60 a 64 
65 a 69 
710a74 
75 y más 
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Con estos datos se puede construir la pirámide de población correspondiente, co- 
locando en la parte izquierda los nombres y en la derecha las mujeres. Al eje ver- 
tical se ha llevado la edad de cinco en cinco años y al eje horizontal la frecuencia 
relativa porcentual (frecuencia relativa x 100). 


75 y más 
HOMBRES MUJERES 70274 
65 a 69 
60 a 64 
55a59 
S0 a 54 
45 a 49 
40 a 44 
31939 
30 a 34 
24 a 29 
20a 24 
15a19 
10a14 
5a9 
0a4 


161412108 6 4 2 0 2 4 6 8 10121416 
El programa 


Los elementos necesarios para diseñar un programa que dé la pirámide de pobla- 
ción son los siguientes: 


1.2 Dos matrices H y M, de dimensión 16 cada una (una dimensión por cada 
grupo de edad), para almacenar los datos de los grupos de edad de hom- 
bres y mujeres, respectivamente. Además, dos totalizadores de hombres y 
mujeres TH y TM. 


2.2 Dos bucles para ir asignando los valores de los distintos grupos de edad: pri- 
mero, a la matriz H y, segundo, a la matriz M. En estos bucles se contaliza- 
rán todos los hombres, TH, y todas las mujeres, TM. 


3.2 Un bucle en el que se calculen los porcentajes (H/TH y M/TM); como estos 
nuevos valores van a estar entre O y 0.20 aproximadamente, los multipli- 
camos por 50 para que se encuentren entre O y 10, pues de esta forma es 
más fácil representarlos gráficamente. 

H(K) = (H(K)/TH) * 50 
M(K) = (M(K)/TM) * 50 
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Por último, en este mismo bucle, estos valores se redondean al entero más 
próximo. Es decir, si, por ejemplo, H(3) = 5.6 tomamos H(3) = 6 y si H(4) = 
= 3.5 tomamos H(4) = 3. 


4.2 El paso siguiente consiste en dibujar los ejes, y en la parte inferior izquierda 
poner HOMBRES y en la parte inferior derecha MUJERES. 


5.2 Por medio de dos bucles anidados se dibujan en la parte izquierda de la grá- 
fica los puntos correspondientes a los hombres y, en la parte derecha, los 
correspondientes a las mujeres. 


1 
10 
20 
30 
40 
50 
60 


REM PIRAMIDE DE POBLACION 
DIM H(16) d 
DIM M(16) 

TH=0 

TM=0 

I=0 
j=4 


En las instrucciones 10 y 20 se dimensionan con el valor 16 las matrices de 
grupos de hombres y de grupos de mujeres H y M. En las 30 y en la 40 se ponen 
a cero los contadores total hombres TH y total mujeres TM. Por último, en la 50 
se pone a cero el contador izquierdo de los grupos de edad y la 60 inicia en cua- 
tro el contador derecho (0 a 4 años). 


70 
80 
90 
100 
110 
120 
130 
140 


FOR K = 1 TO 16 

PRINT “HOMBRES DE”; I; “A”; J; “AÑOS” 
INPUT H(K) 

TH = TH + H(K) 

I=1+5 

J=J+5 

NEXT K 

CLS 


La instrucción 80 imprime el grupo de edad que hay que introducir, introducción 
que se hace en la instrucción 90. En la 100 el total de hombres TH se incrementa 
con el valor introducido en 90. Las instrucciones 110 y 120 incrementan los lí- 
mites del grupo de edad en cinco. La 130 repite el bucle hasta que se introduz- 
can los 16 valores. La 140 borra la pantalla. 
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150 
160 
170 
180 
190 
200 
210 
220 
230 
240 


| = 
J = 
FOR K = 1 TO 16 A 
PRINT “MUJERES DE”; I; “A”; J; “AÑOS” 
INPUT M(K) 

TM = TM + MÍ(K) 

I=1+5 

ales 

NEXT K 

CLS 


Las instrucciones 150 a 160 vuelven a iniciar los límites de grupos de edad | y J, 
para los valores de mujeres. Desde la instrucción 170 a la 230 se repite el bucle 
anterior, sólo que para mujeres. Se termina borrando la pantalla con la instrucción 
240. 


250 
260 
270 
280 
290 
300 
310 
320 


FOR K = 1 TO 16 

H(K) = (H(K)/TH) * 50 

IF H(K) — INT (H(K)) > 0.5 THEN H(K) = H(K) + 1 
H(K) = INT (H(K)) 

M(K) = (M(K)/TM) * 50 

IF M(K) — INT (M(K)) > 0.5 THEN M(K) = M(K) + 1 
M(K) = INT (M(K)) 

NEXT K 


En la instrucción 250 comienza el bucle que termina en la 320. En la 260 se cal- 
cula el porcentaje de cada grupo de hombres respecto del total de hombres; 
como este valor que se obtiene está entre 0 y 1, se multiplica por 50 para luego 
poder hacer la representación gráfica. En la 270 se observa si la parte decimal de 
H(K) es mayor que 0.5; en caso afirmativo, se suma una unidad. La instrucción 
280 deja a H(K) con su parte entera, eliminando la parte decimal. En las instruc- 
ciones 290, 300 y 310 se realiza el mismo proceso para los grupos de mujeres. 


330 
340 
350 
360 
370 
380 
390 
400 


LOCATE 7,21 : PRINT “HOMBRES” 
LOCATE 17,21 : PRINT “MUJERES” 
FOR K = 0 TO 31 

LOCATE K, 19 : PRINT “WN” 

NEXT K 

FOR K = 0 TO 18 

LOCATE 15,K : PRINT “Pl” 

NEXT K 
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Las instrucciones 330 y 340 escriben las palabras HOMBRES y MUJERES en la 
parte inferior del gráfico. El bucle 350-370 dibuja el eje horizontal y el 380-400 el 
eje vertical. 


410 
420 
430 
440 
450 


FOR K = 1T0 16 

FOR L = 0 TO HI(K) 

LOCATE 15-L, 19-K : PRINT “N” 
NEXT L 

NEXT K 


Estos bucles anidados dibujan la parte de la pirámide de población correspon- 
diente a los 16 grupos de edades de hombres. 


460 
470 
480 
490 
500 


FOR K = 1T0 16 

FOR L = 0 TO MI(K) 

LOCATE 15 + L, 19-K : PRINT “P" 
NEXT L 

NEXT K 


Análogamente al bloque anterior, éste dibuja la pirámide correspondiente a los 16 
grupos de edades de mujeres. 


Ejecución del programa 


Si se ejecuta el programa introduciendo como datos los del ejemplo anterior 
(Guatemala - 1960) se obtiene el siguiente resultado en la pantalla: 
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Observación 


Si se quieren mantener los datos junto con el programa, eliminar las instrucciones 
80, 90, 180 y 190, e introducir las siguientes: 


80 READ H(K) 

180 READ M(K) 

510 DATA 321.6, 262.6, 232.2, 201.7, 182.7, 148.4, 114.2, 106.2, 76.1, 
73.3, 56.1, 41.9, 38.1, 19.0, 11.4, 17.1 

520 DATA 344,3, 278.3, 233, 233, 212.4, 165.0, 117.5, 115.5, 90.7, 76.3, 
63.9, 41.2, 39.2, 18.6, 14.4, 18.6 


3. Ordenar cronológicamente 


Presentación del problema 


El siguiente programa sirve para comprobar el acierto en una ordenación cronoló- 
gica de acontecimientos. Serviría también en Literatura para situar en el tiempo a 
ciertos autores, así como para comprobar la ordenación de elementos según ener- 
gías de ionización, radios atómicos, etc. 


El programa 


Inicialmente se introducen estos acontecimientos ordenadamente en la memoria. 
El número de acontecimientos introducidos N puede determinarse a voluntad. El 
número de acontecimientos que salen en pantalla al azar es tres, distintos en- 
tre sí. 


Para evitar escribir demasiado, se colocan las letras A, B y C delante de ellos y 
éstas son las respuestas que deben darse. 


Las respuestas se comparan con el contenido de la memoria y, en caso de error, 
el programa da la respuesta correcta. 


10 INPUT N 
20 DIM AS (N) 
30 DIM C (3) 
40 DIM D$ (3) 


Se declaran variables que se van a utilizar y se introduce el número de aconteci- 
mientos. 


213 


50 FORI=1TON 
60 INPUT AS (1) 
70 NEXT | 


Con este bucle se introducen ordenadamente los distintos acontecimientos. 


75 D$ (1) 
74 D$ (2) 
76 D$ (3) 


"A" 
”g" 
se 0d 


Se almacenan los valores A, B y C en D$ (1), D$ (2), D$ (3), respectivamente. Más 
adelante se compararán con las contestaciones. 


90 S= INT (RND (1) * N + 1) 

100 T= INT (RND (1) * N + 1) 

105 IFT= S THEN GOTO 100 

110 X= INT (RND (1) * N + 1) 

115 IFX=SOR X = T THEN GOTO 110 


Se eligen tres respuestas al azar que sean distintas. 


120 PRINT “A”; A$ (5) 
130 PRINT “B”; A$ (T) 
140 PRINT “C”; A$ (X) 


Se imprimen los tres acontecimientos elegidos con las letras A, B y C delante. 


150 FOR J= 1T0 3 

160 INPUT B$ ()) 

170 PRINT B$ ()) 

180 IF B$ (1) = D$ (1) THEN GOTO 220 
190 IF BS (1) = D$ (2) THEN GOTO 240 
200 IF BS (J) = DS (3) THEN GOTO 260 
210 GOTO 550 

220 B$ (1) = AS (5) 

230 GOTO 270 

240 B$ (1) = AS (T) 

250 GOTO 270 

260 B$ (J) = A$ (X) 

270 NEXT J 
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Pide una contestación A, B o C y después se asigna la contestación A$ (S), A$ (T) 
o A$ (X) a B$ (1), B$ (2) y B$ (3), respectivamente. 


280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 


C(1) = 5 

C(2) = T 

C(3) = X 

FOR M = 1T0 2 
FOR N = 1T0 2 
X = C(N) 

Y = C(N + 1) 

IF X < = Y THEN GOTO 380 
C(N) = Y 

C(N + 1) = X 
NEXT N 

NEXT M 


Ordena S, T y X en forma creciente (ver programa de ordenación de números). 


400 
410 
420 


S = C(1) 
T = C (2) 
X = C (3) 


Coloca en S, T y X el valor correctamente ordenado. 


430 
440 
450 
460 
470 
480 
490 
495 
500 
510 
520 
530 
540 
550 


IF B$ (1) = AS (S) THEN GOTO 450 

GOTO 500 

IF B$ (2) = A$ (T) THEN GOTO 470 

GOTO 500 

IF B$ (3) = A$ (X) THEN GOTO 490 

GOTO 500 

PRINT “ACERTASTE” 

GOTO 90 
PRINT “TE EQUIVOCASTE, EL RESULTADO ES:” 
PRINT A$ (5) 

PRINT AS (T) 

PRINT A$ (X) 

GOTO 90 

END 
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Esta parte del programa compara las respuestas con el contenido de la memoria y 
en caso de error de la respuesta correcta. 

El programa podría mejorarse borrando la pantalla (CLS) entre las instrucciones 
490 y 495, estableciendo un tiempo de espera (con PAUSE, WAIT, o instrucción 
similar según el ordenador) y lo mismo entre 530 y 540 (CLS y PAUSE, WAIT). 
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15. Lengua e idiomas 


1. Ejercicio de acentuación 
Presentación del problema 


Este programa da un texto literario de la obra Platero y yo, de Juan Ramón Jimé- 
nez, omitiendo los acentos. El alumno debe escribir las palabras que están mal es- 
critas. 


El texto es el siguiente: 


“ACARICIE A PLATERO Y, COMO PUDE, LO ENGANCHE A LA CARRETILLA, DE- 
LANTE DEL BORRICO MISERABLE, LO OBLIGUE ENTONCES, CON UN CARINOSO 
IMPERIO, Y PLATERO, DE UN TIRON, SACO CARRETILLA Y RUCIO DEL ATOLLA- 
DERO, Y LES SUBIO LA CUESTA». 


Las palabras solución son: 
ACARICIÉ, ENGANCHÉ, OBLIGUÉ, TIRÓN, SACÓ, SUBIÓ 


El programa 


Lo primero que hay que hacer es almacenar el texto en memoria mediante una 
subrutina, para que aparezca en pantalla todas las veces que sea necesario. Las 
palabras solución se almacenan en una lista B$ de dimensión 6, pues hay seis pa- 
labras. 


Se define una variable A$ para que el alumno dé sus respuestas. 


Por último, se hace una comparación entre la respuesta A$ y la solución B$. Si es 
correcta, se pasa a la siguiente solución y, en caso contrario, se obliga a que el 
alumno responda otra vez. Este proceso se repite hasta que se conteste correcta- 
mente todas las palabras, con lo que el programa se termina. 


1 REM ACENTUACION 

10 PRINT “DEBES DAR LAS PALABRAS QUE LLEVEN ACENTO, DE 
UNA EN UNA” 

20 DIM B$(6) 

30 B5(1)="ACARICIE” 

40 B5$(2)="ENGANCHE” 

50 B$(3)="OBLIGUE” 

60 B$(4)="TIRON” 

70 B$(5)="SACO" 

80 B$(6)="SUBIO” 
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En la instrucción 10 aparece lo que el alumno debe hacer durante la ejecución del 
programa. La instrucción 20 da la dimensión de 6 cadenas de caracteres. Por úl- 
timo, de la 30 a la 80 se introducen en la lista B$ las soluciones del problema. 


90 FORI|=1T0 6 
95 PRINT 

100 GOSUB 500 
110 INPUT AS 


La instrucción 90 abre un bucle para las seis respuestas válidas. La 100 envía a la 
subrutina; ésta escribe el texto que se debe acentuar. La 110 sirve para que el 
alumno introduzca su respuesta. 


160 IF A$=B$ (1) THEN GOTO 190 * 
170 PRINT “ERROR; PRUEBA DE NUEVO” 
180 GOTO 110 


La instrucción 160 compara la respuesta del alumno A$ con la solución, B$ (1); si 
no es correcta, pasa el control a la instrucción siguiente (instrucción 170), la cual 
le anuncia que es errónea la respuesta y que pruebe de nuevo; por ello, la ins- 
trucción 180 envía el control del programa a la 110, para que el alumno intro- 
duzca otra respuesta (hasta que no sea correcta no se sale de este bucle). Si la 
respuesta es la correcta, la instrucción 160 envía el control a la 190, 


190 CLS 

200 PRINT “CORRECTO, BUSCA LA PALABRA SIGUIENTE” 
210 NEXT | 

220 PRINT “NO BUSQUES MAS; HAS ENCONTRADO TODAS” 
230 END 


La 190 borra toda la pantalla y la 200 anuncia que la respuesta es correcta y que 
el alumno se prepare para la siguiente. En la 210 se termina el bucle comenzado 
en la instrucción 90. 


El fin del ejercicio se anuncia con la instrucción 220 y el programa termina con la 
instrucción 230 END para que no se ejecuten las instrucciones siguientes, corres- 
pondientes a la subrutina. 


500 PRINT “ACARICIE A PLATERO Y, COMO PUDE, LO ENGANCHE A 
LA CARRETILLA, DELANTE DEL BORRICO MISERABLE, LE OBLI- 
GUE ENTONCES, CON UN CARIÑOSO IMPERIO, Y PLATERO, DE 
UN TIRON, SACO CARRETILLA Y RUCIO DEL ATOLLADERO, Y 
LES SUBIO LA CUESTA”. 

510 RETURN 
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Estas dos últimas instrucciones componen la subrutina, llamada desde la instruc- 
ción 100. La 500 escribe el texto de estudio, que se repetirá seis veces a lo largo 
de la ejecución del programa (una para cada respuesta), y la 510 envía el control 
a la instrucción siguiente: a la 110. 


2. Conjugación de verbos 


Presentación del problema 


Los verbos regulares se dividen en tres grupos, atendiendo a su terminación: en 
AR, en ER y en IR. La conjugación de todos los verbos de un mismo grupo es bas- 
tante sencilla pues sus terminaciones (desinencias) son iguales. Así, por ejemplo, 
para conjugar el presente de indicativo del verbo TARDAR basta con añadir a la 
parte radical TARD las terminaciones 


NOSOTROS 
VOSOTROS 
ELLOS 


El programa 


El programa siguiente da el presente de indicativo, el pretérito imperfecto o co- 
pretérito, el pretérito indefinido o pretérito y el futuro imperfecto o futuro (todos 
estos tiempos forman el modo indicativo simple de un verbo de la primera conju- 
gación, o sea de los terminados en AR). 


Lo primero que hace el programa es pedir el verbo y almacenarlo en A$. A conti- 
nuación calcula la longitud de A$ y extrae una subcadena B$, en la que falta la 
terminación AR. 
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Por último, se imprimen los distintos tiempos del modo indicativo simple. 


1 


REM MODO INDICATIVO SIMPLE DE LA PRIMERA CONJUGACION 


10 INPUT “VERBO”, A$ 
20 L=LEN(A$) 
30 B$=LEFT$ (AS; L-2) 


En la instrucción 10 se introduce el verbo, en la 20 se calcula su longitud y en la 
30 se extrae una subcadena izquierda de longitud L-2. Luego B$ es la parte radi- 
cal del verbo. 


PRINT “PRESENTE DE INPICATIvOs 
PRINT “YO”, B$, “O” 

PRINT “TU”, B$; “AS” 

PRINT “EL”, B$; “A” 

PRINT “NOSOTROS”, B$; “AMOS” 
PRINT “VOSOTROS”, B$; “AIS” 
PRINT “ELLOS”, B5; “AN” 


El bloque 40-100 imprime el presente de indicativo. Por ejemplo, la instrucción 50 
imprime YO, salta a la siguiente zona (debido a la coma) y escribe B$ junto a O, 
pues al estar separados por punto y coma se imprimen seguidos. 


110 
120 
130 
140 
150 
160 
170 


PRINT “PRETERITO IMPERFECTO O CO-PRETERITO” 
PRINT “YO”, B5; “ABA” 

PRINT “TU”, B$; "ABAS” 

PRINT “EL”, B$, "ABA" 

PRINT “NOSOTROS”, B$; “ABAMOS” 

PRINT “VOSOTROS”, B$; “ABAIS” 

PRINT “ELLOS”, B5; “ABAN"” 


Este bloque imprime el pretérito imperfecto o co-pretérito. 


180 
190 
200 
210 
220 
230 
240 


PRINT “PRETERITO INDEFINIDO O PRETERITO” 
PRINT “YO”, BS; “E” 

PRINT “TU”, B$; “ASTE” 

PRINT “EL”, B$; “0” 

PRINT “NOSOTROS”, B5; “AMOS” 

PRINT “VOSOTROS”, B$; “ASTEIS” 

PRINT “ELLOS”, B5; “ARON" 
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Las instrucciones 180-240 imprimen el pretérito indefinido o pretérito. 


250 PRINT “FUTURO IMPERFECTO O FUTURO” 
260 PRINT “YO”, A$; “E” 

270 PRINT “TU”, AS; “AS” 

280 PRINT “EL”, A$; “A” 

290 PRINT “NOSOTROS”, A$; "EMOS" 

300 PRINT “VOSOTROS”, A$; “ElS” 

310 PRINT “ELLOS”, A$; “AN” 


En el bloque 250-310 se toma como parte radical el verbo completo A$, ya que 
la conjugación del futuro incluye en la parte radical la terminación AR. 


3. Identificación de preposiciones inglesas 
Planteamiento del problema 


Se trata de que vayan apareciendo en pantalla una serie de frases en inglés, a las 
que le faltan las preposiciones. El que ejecute el programa deberá dar una prepo- 
sición (respuesta) y la pantalla le comunicará si es correcta o no. Si lo es, pasará a 
la siguiente frase; en caso contrario, volverá a aparecer en la pantalla la misma 
frase. 


Las preposiciones 

Las frases son son 
JOHN DREAMT ... MARY LAST NIGHT OF 
'M LOOKING FORWARD ... MEETING YOU NEXT 

WEEKEND A TO 
MARY IS THINKING ... JOHN AT THE MOMENT OF 
DON'T STAND THERE ... THE DOOR, COME ... 

AND SIT ... AT-IN-DOWN 
YOU CAN COUNT ... ME, | NEVER LET MY 

FRIENDS ... ON-DOWN 
El programa 


El programa es abierto, es decir, sirve para cualquier número N de frases (en este 
caso N = 5) y no únicamente en inglés sino en cualquier otro idioma; basta con 
dar las frases y preposiciones en el idioma elegido. 
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La entrada de los datos, frases y preposiciones se realiza por medio de un bucle 
que tiene la instrucción READ A$, B$, almacenándose en Af las frases y en B$ las 
preposiciones. 


Previamente hay que dimensionar las listas A$ y B$. 


El paso siguiente es definir dos variables Y$ y Z$ que, por medio de su impresión 
en pantalla, informen si la respuesta elegida es correcta o no. 


Por último, el programa ha de contener un bucle en el que se dé una de las 
frases, se introduzca una respuesta y ésta sea comparada con la solución. No se 
pasará a la frase siguiente hasta que no se acierte la respuesta. 


1 REM PREPOSICIONES EN INGLES 
10 INPUT N 
20 DIM A$(N) : DIM B$ (N) 
30 FOR|=1TON 
40 READ A$(1), B$(1) 
50 NEXT | 
60 Y$= "NO ES CORRECTA, REPITE” 
70 Z$ = "CORRECTA, AQUI TIENES OTRA FRASE” 


En la instrucción 10 se introduce el número de frases que hay. La 20 dimensiona 
las N frases y las N preposiciones. De la 30 a la 50, el bucle lee las frases y las 
preposiciones. Por último, en la 60 y 70 se definen dos variables, Y$ y X$, que 
servirán para decir al alumno si ha acertado o no. 


80 FOR!|I=1TON 

90 PRINT “PON LA PREPOSICION QUE FALTA Y SI HAY VARIAS 
DEJA UN ESPACIO ENTRE ELLAS” 

100 PRINT AS(1) 

110 INPUT X$ 


En la instrucción 80 comienza la serie de instrucciones que llevarán al alumno a 
contestar las distintas cuestiones. En la 90 se le dice qué es lo que debe respon- 
der y cómo, en el caso de varias preposiciones. La instrucción 100 le da la frase 
que tiene que completar con la preposición o preposiciones correctas. Por último, 
en la 110 se le pide la respuesta X$. 


160 IF X$ = B$(1) THEN GOTO 190 
170 PRINT Y$ 
180 GOTO 100 
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La instrucción 160 compara la solución B$ con la respuesta X$ dada por el 
alumno; si no es correcta, la 170 le dice que repita, por medio de la variable Y$ 
definida en 60; y la 180 devuelve el control a la instrucción 100, con lo que se da 
la oportunidad al alumno de repetir su respuesta X$. Este proceso no termina 
mientras no dé la respuesta correcta. 


190 CLS 

200 PRINT Z$ 

210 NEXT | 

220 PRINT “¡ENHORABUENA! LLEGASTE A LA META” 


A la instrucción 190 se llega cuando el alumno ha acertado la preposición corres- 
pondiente; esta instrucción borra toda la pantalla. Posteriormente, se da al 
alumno, por medio de la variable Z$ (en la instrucción 200), la información de ha- 
ber acertado con su respuesta. Mediante la instrucción 210 se repite el ciclo con 
otra frase. Con la instrucción 220 se le comunica al alumno el fin de las frases y, 
por tanto, la terminación del programa. 


Este programa sirve para cualquier número de frases, pero si se quiere utilizar 
para el ejemplo presentado, hay que introducir los siguientes datos: 


230 DATA 5 

240 DATA "JOHN DREAMT ... MARY LAST NIGHT”, “OF” 

250 DATA “I'M LOOKING FORWARD ... MEETING YOU NEXT WEE- 
KEND”, “TO” 

260 DATA "MARY IS THINKING ... JOHN AT THE MOMENT”, “OF” 

270 DATA “DON'T STAND THERE ... THE DOOR, COME ... AND SIT ...”, 
“AT IN DOWN” 

280 DATA “YOU CAN COUNT ... ME, | NEVER LET MY FRIENDS ...”, 
“ON DOWN” 

290 END 


4. Repasando la morfología y sintaxis latina 


Representación del problema 


En el texto siguiente: 


“ACEPT ORATIONE EORUM, CAESAR OBSIDES TRADI IMPERAT EOSQUE ADUCCI 
IlUBET. HIS ADDUCTIS, UT IMPERAVERAT, ARBITROS INTER CIVITATES DAT, QUI 
LITEM AESTIMENT” 
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el alumno debe buscar los sustantivos (excluidos los nombres propios) y, una vez 
encontrados todos, ponerlos en el mismo caso, pero en número distinto al que 
tienen en el texto. 


Las soluciones son: 


ORATIONE ORATIONIBUS 
OBSIDES - OBSIDEM 
ARBITROS ARBITRUM 
CIVITATES CIVITATEM 
LITEM LITES 


(CAESAR ES NOMBRE PROPIO) 


El programa 


Se utilizan dos variables B$(I) y C$(I) para introducir en la memoria del microorde- 
nador las palabras solución. 


El texto, incluido en dos variables T$ y U$, se almacena en una subrutina para no 
sobrepasar el límite de 255 caracteres. 


En un bucle se comprueba si cada respuesta es correcta. Cuando se han encon- 
trado todos los sustantivos, se pide que se pongan en distinto número. De nuevo, 
las respuestas se comparan con las palabras solución y en caso de equivocación se 
da el resultado correcto. 


1 REM LATIN 
10 READ N 
20 DIM B$ (N) 
25 DIM C$ (N) 
30 FORI|I=1TON 
40 READ Bf (1) 
50 READ C$ (1) 
60 NEXT | 


En estas instrucciones (en combinación con las instrucciones DATA del final del 
programa) se introducen los sustantivos en distintos números C$ (|). 
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70 
80 


90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 


GOSUB 500 

PRINT “BUSCA TODOS LOS SUSTANTIVOS EN EL TEXTO SI- 
GUIENTE” 

INPUT A$ 

FOR | = 1TON 

IF AS = B$ (1) THEN GOTO 130 

GOTO 160 

PRINT “ACERTASTE” 

C=C+1 

GOTO 190 

NEXT | 

PRINT “TE EQUIVOCASTE, BUSCA DE NUEVO” 
GOTO 90 

IFC < N THEN GOTO 90 

PRINT “YA HAS ENCONTRADO TODOS” 


Una vez puesto el texto en pantalla (instrucción 70 GOSUB 500) se piden los sus- 
tantivos (instrucciones 80 y 90) y en el bucle siguiente (instrucciones 100 a 160) 
se comprueba si la palabra dada es alguna de las almacenadas en la tabla B$ (1); 
si no es, se pide una nueva respuesta, y si es, se suma uno al contador C. 
Cuando C coincide con el número de sustantivos existentes (instrucciones 160 a 
170) se imprime un mensaje. 


210 
220 
230 


240 
250 
260 
270 
280 
290 
300 
310 
320 


FOR K = 1 TO 2000 : NEXT K 

CLS 

PRINT “PON LOS SUSTANTIVOS ANTERIORES EN EL MISMO 
CASO PERO EN DISTINTO NUMERO” 

FOR | = 1TO N 

PRINT B$ (1) 

INPUT D$ 

IF D$ = C$ (1) THEN GOTO 300 

PRINT “TE EQUIVOCASTE, LA RESPUESTA ES:”; C$ (1) 
GOTO 310 

PRINT “DE ACUERDO” 

NEXT | 

STOP 


Después de un tiempo de espera (instrucción 210), se borra la pantalla (instruc- 
ción 220) y se hace una nueva pregunta (230). La respuesta se compara con el 
contenido de la memoria C$ (1). 
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500 
510 
520 


530 
540 


T$ = "ACCEPTA ORATIONE EORUM, CAESAR OBSIDES TRADI 
IMPERAT EOSQUE ADUCCI IUBET” 

U$ = "HISADDUCTIS, UT IMPERAVERAT, ARBITROS INTER 
CIVITATES DAT, QUI LITEM AESTIMENT” 

PRINT T$ 

PRINT U$ 

RETURN 


Con la subrutina (500-540) se almacena el texto en dos variables, T$ y U$, y a 
continuación se imprime 


590 
600 


610 
620 


DATA 5 

DATA ORATIONE, ORATIONIBUS, OBSIDES, OBSIDEM, ARBI- 
TROS, ARBITRUM 

DATA CIVITATES, CIVITATEM, LITEM, LITES 

END 


En las instrucciones DATA se incluyen las distintas soluciones. 
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16. Control de una pequeña 
biblioteca 


Presentación del problema 


El almacenamiento y el tratamiento de información se hace habitualmente con 
dispositivos que utilizan discos o disquetes. Sin embargo, si la cantidad de infor- 
mación no ocupa excesiva memoria, se puede utilizar un simple grabador, repro- 
ductor de casetes. 


El programa 


El programa siguiente sirve para controlar una pequeña biblioteca, para lo cual 
permite llevar a cabo las siguientes acciones: 


— Crear el archivo de libros. 

— Cargar el archivo desde casete. 
— Añadir libros al archivo ya creado. 
— Buscar libros de una materia. 

— Buscar libros de un autor. 


— Buscar el autor de un título. 


Para construir el programa se define la variable L$ (1) mediante la concatenación 
de la variable A$ (1) que incluye un número clave con dos cifras (correspondiente 
a la materia), la variable B$ (I) con un número clave (correspondiente a la subma- 
teria), la variable C$ (1) que incluye el autor de la obra (para el que se reservan 30 
caracteres) y la variable D$ (I) que incluye el título del libro (para el que también 
se reservan 30 caracteres). 


Una vez ordenados los libros introducidos, tomando diversos fragmentos de la va- 


riable L$ (1), con la ayuda de las instrucciones LEFT$ ( ), MID$ ( ) y RIGHTS ( ) se 
puede localizar materia, autor o título según se desee. 
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1 REM BIBLIOTECA 

10 CLEAR 3000 

20 DIM Af$ (50), B$ (50), C$ (50), D$ (50), E$ (50), L$ (50) 
30 DIM C(50), D(50) 

40 CLS 

50 PRINT “BIBLIOTECA” 

60 PRINT : PRINT : PRINT - 

70 PRINT “ESCRIBE:” 

80 PRINT 

90 PRINT ”1 - PARA CREAR ARCHIVO” 
100 PRINT “2 - PARA MODIFICAR ARCHIVO” 
110 PRINT ”3 - PARA CARGAR ARCHIVO” 
120 PRINT “4 - PARA BUSCAR LIBROS DE UNA MATERIA” 
130 PRINT “5 - PARA BUSCAR LIBROS DE UN AUTOR” 
140 PRINT “6 - PARA BUSCAR EL ÁUTOR DE UN TITULO” 
150 PRINT “7 - PARA ACABAR” 
160 INPUT A 
170 IFA<1ORA > 7 THEN GOTO 40 
180 ON A GOSUB 1000, 2000, 3000, 4000, 5000, 6000, 7000 
190 GOTO 40 


La instrucción 10 limpia la parte de memoria correspondiente a las cadenas, lo 
que permite introducir el número de caracteres necesarios. 


El resto de las instrucciones sirve para reservar posiciones de memoria, y presenta 
las opciones disponibles. Cada una de estas opciones corresponde a una subru- 
tina diferente. 


1000 REM ARCHIVO 

1010 CLS 

1020 |=1 

1030 LOCATE 1,5 : INPUT “CLAVE”; AS (1) 

1040 LOCATE 1,7 : INPUT “SUBCLAVE”: B$ (1) 

1050 LOCATE 1,9 : INPUT “AUTOR”; C$ (1) 

1060 LOCATE 1,11 : INPUT “TITULO”; D$ (1) 

1070 LOCATE 1,20 : INPUT “DE ACUERDO (S/N)”; F$ 
1080 IF F$= “N” OR F$ = “n” THEN GOTO 1030 


Con estas instrucciones se introducen los datos correspondientes a cada libro. 
1100 C (1) = LEN (C$(1)) 
1110 D(I) = LEN(DS(1)) 


1120 FOR J = C(I) + 1T0O 30 
1130 C$ (1) = C$ (1) + ”” 
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1140 NEXT J 

1150 FOR K =D (I) + TO 30 
1160 D$ (1) = D$ (1) + ”” 
1170 NEXT K 


Aquí añaden espacios en blanco a C$ (1) y DS (1) hasta que alcancen una longitud 
de 30. 


1180  L$ (1) = AS (1) + BS (1) + C$ (1) + DS (1) 
1190 INPUT “MAS LIBROS”; F$ 

1200 IF F$= “N” OR F$ = “n” THEN GOTO 1300 
1210 |=1+1 

1220 GOTO 1030 


Este bloque de instrucciones concatena las cadenas A$ (1), B$ (1), C$ (1) y DS (1), 
formando la cadena L$ (1). Además permite introducir un nuevo libro, en cuyo 
caso se incrementa la variable | en una unidad y a continuación se reinicia el pro- 
ceso. 


1300 FORJ=1TOI-1 
1310 FORK =J+1T01 
1320 IF L$ (1) < L$ (K) THEN GOTO 1370 


) 
1330 X$ = L$ (J) 
1340  L$ (1) = L$ (K) 
1350 L$ (K) = X$ 
1360 NEXT K 
1370 NEXT J 


E 


pr 


este grupo de instrucciones se ordenan las variables L$ (1). 


1400 CLS , 

1410 PRINT “PREPARAR EL CASETE PARA GRABAR” 
1415 PRINT : PRINT : INPUT “PREPARADO S/N”; F$ 
1420 IF F$= "N” OR F$ = “n” THEN GOTO 1400 
1430 OPEN "CAS : LIBR” FOR OUTPUT AS 4 1 
1440 FORJ=1T01| 

1450 PRINT 4 1, L$ (J) 

1460 NEXT J 

1470 CLOSE + 1 

1480 PRINT "PULSA UNA TECLA” 

1490 A$ = INKEY$ 

1500 IF A$ = ”” THEN GOTO 1490 

1510 RETURN 
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Mediante estas instrucciones, los datos almacenados en L$ (I) se graban en un ca- 
sete. Las instrucciones 1480 y 1510 permiten volver al punto de partida. 


2000 
2010 
2020 


REM MODIFICAR 
CLS 
GOTO 1030 


Volviendo a la instrucción 1030 pueden añadirse libros a un archivo ya creado. 


3000 
3005 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 


REM CARGAR 

CLS , 

PRINT “PREPARAR EL CASETE PARA CARGAR” 
PRINT : PRINT : PRINT : INPUT “PREPARADO S/N”; F$ 
IF. FS = “N” OR F$ = “n” THEN GOTO 3010 
OPEN "CAS : LIBR” FOR INPUT AS 4 1 

| =1 

IF EOF(1) = —1 THEN CLOSE +1 : RETURN 
INPUT 4 1, L$ (1) 

|=1+1 

GOTO 3060 


Se van leyendo libros desde el casete, hasta que en la instrucción 3060 se detecta 
el final del archivo, pasándose entonces al menú principal. 


4000 
4005 
4010 
4020 
4030 


4040 
4050 
4060 
4070 
4080 


REM MATERIA 

CLS 

INPUT “INTRODUZCA LA CLAVE DE UNA MATERIA”; H$ 
FOR J = 1TO |! 

IF H$ = LEFTS (L$(), 2) THEN PRINT LEFTS(L$(), 3) : PRINT 
MID$ (L$(J), 4,30) : PRINT RIGHTS$(L$(J), 30) 

NEXT J 

PRINT “PULSA UNA TECLA PARA VOLVER” 

A$ = INKEY$ 

IF AS = “” THEN 4060 

RETURN 


Se localizan los libros correspondientes a una clave determinada. 
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5000 
5005 
5010 
5020 
5030 
5040 
5050 
5060 
5070 


5080 
5090 
5100 
5110 
5120 


REM AUTOR 

cs 

INPUT "INTRODUCE EL NOMBRE DEL AUTOR”; N$ 
U = LEN (N$) 

FOR X = U + 1T0 30 

N$ = N$ + ”” 

NEXT X 

FOR J = 1TO! 

IF MIDS (L$(1), 4,30) = N$ THEN PRINT LEFT$ (L$(J), 3) : PRINT 
MID$ (L$ (3), 4,30) : PRINT RIGHT$ (L$ (1), 30) 
NEXT J 

PRINT “PULSA UNA TECLA PARA VOLVER” 

A$ = INKEY$ 

IF A$ = ”” THEN GOTO 5100 

RETURN 


Se introduce el nombre de un autor al que una vez determinada su longitud se le 
añaden espacios en blanco hasta completar 30 carácteres. El bucle comprendido 
entre 5060 y 5080 localiza el autor en la cadena L$(1). 


6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 


6090 
6100 
6110 
6120 
6130 


REM TITULO 

CLS 

INPUT “INTRODUZCA TITULO”; T$ 

V = LEN(T$) 

FOR R = V + TO 30 

15 =T5 +" 

NEXT R 

FOR J = 1T0O! 

IF TS = RIGHTS(L$(J), 30) THEN PRINT LEFTS(L$(J), 3) : PRINT 
MID$(L$(J), 4,30) : PRINT RIGHTS$(L$(J), 30) 
NEXT J A 

PRINT “PULSA UNA TECLA” 

A$ = INKEY$ 

IF A$ = “” THEN GOTO 6110 

RETURN 40 


Con este último bloque de instrucciones se introduce un, título, se completa hasta 
30 caracteres y, en caso de localizarse entre los existentes, se imprimen la clave, 
subclave, el autor y el propio título. 


7000 
7010 


REM FIN 
END 
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17. Un programa para escribir 
y otro para dibujar 


1. Programa para escribir 


El programa siguiente permite llevar a cabo estas acciones: 

— Escribir una página de 20 líneas con 38 caracteres cada una. 

— Ver una página completa en la pantalla. 

— Guardar una página en un casete. 

— Reproducir una página que se ha grabado previamente en un casete. 


1 REM ESCRIBIR 
10 CLEAR 3000 
20 DIM D$(22, 40) 
30 CLS 
40 PRINT “ESCRIBE :” 
50 PRINT “1-PARA ESCRIBIR PAGINA” 
60 PRINT “2-PARA VER PAGINA” 
70 PRINT “3-PARA GUARDAR PAGINA EN CASETE” 
80 PRINT “4-PARA REPRODUCIR PAGINA DE CASETE” 
90 PRINT “5-PARA ACABAR” 
100 INPUT A 
110 IFA<1ORA > 5 THEN GOTO 30 
120 ON A GOSUB 1000, 2000, 3000, 4000, 5000 
130 GOTO 30 


Con estas instrucciones se reserva espacio en la memoria y se presentan las op- 
ciones de que dispone el programa. 
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1000 REM ESCRIBIR 

1005 PRINT “ESCRIBA/PARA ACABAR” : FOR K = 1 TO 2000 : NEXT K 
1010 CLS 

1020 l|=1:J]=1 

1030 IFJ<1ORJ= 38 THEN J = 1 

1040 A$ = INKEY$ 

1050 IF A$ = ““THEN GOTO 1040 

1055 IF A$ = ”” THEN CLS : RETURN 

1060 IF ASC(A$) = 13 THEN J = 38; GOTO 1110 

1070 IF ASC(A$) = 8 THEN PRINT CHR$(29); ””; CHR$(29); : J =J — 


1: D$ (1, J) = “" : GOTO 1030 
1080 D$(l,J) = A$ : PRINT DS$(l, J): 
1085 J=J+1 


1090 IF J = 38 THEN GOTO 1110 
1100 GOTO 1040 

1110 PRINT 

1115 1F|= 20 THEN CLS : RETURN 
1120 |=1+1 

1130 GOTO 1030 


En la instrucción 1040 se almacenan en la cadena A$ los caracteres que se van 
pulsando. 


Cuando se pulsa la tecla | RETURN | (código 13) se produce un salto de línea. 
Cuando se pulsa la tecla (DELETE o similar), con código 8, se borra el ca- 
rácter que se acaba de escribir, ya que PRINT CHR$(29) hace que el cursor retro- 
ceda. 


Por otra parte, en la cadena D$(l, J) se van almacenando los distintos caracteres 
pulsados. 


2000 REM VER PAGINA 

2005 PRINT “CUANDO ACABES PULSA UNA TECLA PARA SEGUIR” : 
FOR K = 1 TO 1000 : NEXT K 

2010 CLS 

2015 FOR!|= 1T0 20 

2020 FOR J = 1 TO 38 

2030 PRINT DS$(l, J); 

2040 NEXT J 

2045 PRINT 

2050 NEXT | 

2070 A$ = INKEY$ 

2080 IF A$ = "” THEN GOTO 2070 

2090 CLS 

2100 RETURN 
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Los caracteres almacenados en D$(l, J) se imprimen en la pantalla. Al pulsar una 
tecla se vuelve al menú inicial. 


3000 REM GUARDAR PAGINA 

3010 CLS 

3020 PRINT “PREPARAR EL CASETE PARA GRABAR” 
3030 PRINT : PRINT : INPUT “PREPARADO(S/N)”; F$ 
3040 IF F$ = “N" OR F$ = “n” THEN GOTO 3020 
3050 OPEN “CAS : PAG” FOR OUTPUT AS +1 

3060 FOR |= 1TO 20 

3070 FOR J = 1 TO 38 

3080 PRINT + 1, D$(l, J) 

3090 NEXT J 

3100 NEXT | 

3110 CLOSE + 1 

3120 RETURN 


Los caracteres almacenados en D$(l, J) se almacenan en un casete. 


4000 REM CARGAR 

4010 CLS 

4020 PRINT “PREPARAR EL CASETE PARA PRODUCIR” 
4030 INPUT “PREPARADO(S/N)”; F$ 

4040 IF F$= "N” OR F$ = “n” THEN GOTO 4020 
4050 OPEN "CAS : PAG” FOR INPUT AS 4 1 

4060 FOR != 1 TO 20 

4070 FOR J = 1 TO 38 

4080 INPUT 4 1, D$(l, J) 

4090 NEXT J 

4100 NEXT | 

4110 CLOSE + 1 

4120 RETURN 


Los caracteres almacenados en un casete se introducen en la memoria del mi- 
croordenador en las variables D$(l, J). Para ver la página puede utilizarse a conti- 
nuación la opción 2. 


5000 REM ACABAR 
5010 END 
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2. Programa para dibujar 


El programa siguiente permite realizar estas acciones: 


— Dibujar en la pantalla con la ayuda de un joystick (moviéndolo en ocho direc- 
ciones) o, inicialmente, usando las teclas de movimiento del cursor. Borrar 
cuando sea necesario y almacenar en la memoria del microordenador los 
puntos de un dibujo que hayan resultado satisfactorios. 


— Ver un dibujo en la pantalla cuando los puntos que van a ayudar a trazarlo es- 
tán almacenados en la memoria. 


— Guardar en un casete un dibujo que apareció en la pantalla. 


— Cargar de un casete los puntos necesarios para reproducir un dibujo en la 
pantalla de un microordenador. 


DIM X(100), Y(100) 

CLS 

PRINT “DIBUJAR” 

PRINT “1-PARA DIBUJAR” 

PRINT “2-PARA VER DIBUJO” 

PRINT “3-PARA GUARDAR DIBUJO EN CASETE” 
PRINT “4-PARA CARGAR DIBUJO EN CASETE” 
PRINT “5-PARA ACABAR” 

INPUT A 

IFA < 1 OR A > 5 THEN GOTO 20 

ON A GOSUB 1000, 2000, 3000, 4000, 5000 
GOTO 20 


Con este bloque de instrucciones se reserva espacio en la memoria para almace- 
nar las coordenadas de los puntos que formarán el dibujo. Se presentan después 
las opciones de que dispone el programa. 


1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 


REM DIBUJAR 

CLS 

PRINT “ESCRIBE D : DIBUJAR” 
PRINT “ESCRIBE B : BORRAR” 
PRINT “ESCRIBE V : GUARDAR PUNTO” 
PRINT “ESCRIBE S PARA SALIR” 
PRINT “PULSA UNA TECLA” 

C$ = INKEY$ 

IF C$ = ”” THEN GOTO 1070 
SCREEN 2 
X=128:1=8:L=031=1 
A$ = INKEY$ : A = STICK (0) 
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1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1340 


IF AS = “” AND A = 0 THEN GOTO 1110: 

IF L = O THEN PSET(X, Y) 

IF L = 1 THEN PRESET(X, Y) 

IF L = 2 THEN X(J) = X : Y(J) = Y :J] =J + 1: PSET (X, Y) : L= 
IFA = 1 THEN Y = Y — 
IFA = 2 THEN Y = 
IFA = 3 THEN X 
IFA = 4 THEN X 
IFA = 5 THEN Y 
IFA = 6 THEN Y 
IFA = 7 THEN X 
IFA = 8 THEN X 
IFX > 255 THEN 
IFX < 1 THEN X 
IF Y < 1 THEN Y 
IF Y < 190 THEN 
IF AS = "B” THE 
IF AS = "D” THEN L = 
IFAS = “V” THEN L = 

IF AS = “S” THEN SCREEN O : RETURN 
GOTO 1110 


X=X+1 


E 
É 
+ 
+1:X=X-1 


1 

¡e 

1 

1: Y=Y+1 
dE 

1 

li 


A A 


Y 
X 
X 
+ 
Y 

X 

X 

= 255 
1 

1 


ZÉ <p momo mo" 
pr 


190 

1 

0 
2 


0 


Con la ayuda de las teclas del cursor (STICK(0)) o con el joystick (STICK(1)) o 
(STICK(2)) se desplaza el cursor en una pantalla de alta resolución. 


Si previamente se ha pulsado D, el cursor deja huella de su paso, y si se ha pul- 
sado B, dicha huella queda borrada. 


Cuando se pulsa V, se almacena en las variables X(l) e Y(I) el punto donde se en- 
cuentra el cursor. Pulsando S el proceso finaliza. 


2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 


2080 
2090 
2100 
2110 
2120 


REM VER DIBUJO 

CLS 

PRINT “PULSA S PARA SALIR” 
FOR K = 1 TO 1000 : NEXT K 


SCREEN 2 
CLS 
FOR | = 0TO ] 


IF X(1 + 1) = 0 AND Y(I + 1) = O THEN X(Y + 1) = X(I) : Y(Y 
+ 1) = Y (1) 

LINE(X(1), Y(0) — (X(1 + 1), Y(l + 1) 

NEXT | 

A$ = INKEY$ 

IF A$ = ”S” THEN CLS : CLOSE + 1 : SCREEN O : RETURN 20 
GOTO 2100 
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Los puntos almacenados (después de haber pulsado V) se unen entre sí en la pan- 
talla del microordenador mediante tramos rectos. Este proceso se termina cuando 
se llega a una coordenada nula para X e Y. 


3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3120 
3130 
3140 
3150 
3160 


REM GUARDAR DIBUJO 

CLS 

PRINT “PREPARAR EL CASETE PARA GRABAR” 
PRINT : PRINT : INPUT “PREPARADO (S/N)”; F$ 
IF. F$ = “N” OR F$ = “n” THEN GOTO 3020 
OPEN "CAS : DIB” FOR OUTPUT AS 4 | 

FOR! =0T0J 

PRINT % 1, X(1) 

PRINT 4 1, Y() 

NEXT | 

CLOSE 1 

PRINT “PULSA UNA TECLA” 

A$ = INKEY$ 

IF A$ = “” THEN GOTO 3130 

CLS 

RETURN 


Las coordenadas X(I), Y(I) almacenadas en la memoria del microordenador se gra- 
ban en un casete. 


4000 
4010 
4020 
4030 
4040 
4050 
4060 
4070 
4080 
4090 
4100 
4110 


REM CARGAR DIBUJO 
CLS 
PRINT “PREPARA EL CASETE PARA REPRODUCIR” 
PRINT : PRINT : INPUT “PREPARADO(S/N)”; F$ 
IF. FS = “N” OR F$ = “n” THEN GOTO 4020 
OPEN "CAS : DIB” FOR INPUT AS 4 1 

= 1 
IF EOF(1) = -1 THEN CLOSE + 1 : RETURN 
INPUT 4% 1, X(1) 
INPUT 4 1, Y(1) 
|=1+1 
GOTO 4070 


Desde un casete se llevan a la memoria del microordenador las coordenadas de 
los puntos de un dibujo, y dichas coordenadas se almacenan en las variables X(l) 


e Y(1). 
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Cuando en la instrucción 4070 se detecta el final del archivo, acaba todo el pro- 
ceso. 


5000 REM ACABAR 
5010 END 
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